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Riassunto. — Si calcola lo sparpagliamento integrale di uno sciame nel- 
Varia, mediato sullo sviluppo longitudinale, in base alla teoria della ca- 
scata elettrofotoniea sviluppata, in approssimazione unidimensionale, 
dal gruppo di Padova. Le espressioni ricavate descrivono, in accordo con 
i risultati di ROBERG e NORDHEIM, lo sparpagliamento nella regione delle 
basse e medie energie. Nella regione delle alte energie esse si convertono 
nelle relazioni stabilite da BORSELLINO. | 


Seguendo i criteri e le approssimazioni usate in un precedente lavoro (*), 
abbiamo valutato, alle basse energie, le grandezze integrali che definiscono lo 
sparpagliamento, laterale degli sciami mediato sullo sviluppo longitudinale, 
ammettendo che la componente elettrofotonica in approssimazione verticale 
sia descritta dagli spettri integrali di bassa energia calcolati dal gruppo di 
Padova (2) (2) (*). Le notazioni non esplicitate che qui compaiono corrispon- 
dono pertanto a quelle usate nei citati lavori. 

Detti TIE,, 1350) & T(E,, W,t) gli spettri integrali degli elettroni e dei 
fotoni di energia E < m e, rispettivamente, W< M, lo sparpagliamento inte- 
grale, mediato sullo sviluppo longitudinale dello sciame, risulta definito dalle 
seguenti relazioni che, per N= 0,1, 2, danno rispettivamente il valore medio 
quadratico dell’angolo di deviazione, dell’accoppiamento an golare laterale e della 


(1) G. Poranr e C. VILLI: Nuovo Cimento, 9, 825 (1952). 

(2) N. DALLAPORTA € E. CLEMENTEL: Nuovo Cimento, 8, 235 (1946). 

(3) N. DALLAPORTA € G. Poranı: Nuovo Cimento, 4, 1 (1947). 

(4) E. CLEMENTEL e L. FABBRICHESI: Nuovo Cimento, 6, 399 (1949); 5, 77 (1948). 
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deviazione laterale: 


(OR, — Ay B)[(O™*R* ys, (Be; E, t)dt, componente elettronica 
0 
(1) 
(OR y= UE; W)| OR, (Bo, W,t)dt, componente fotonica 


0 


in cui F(B,, E) e G(E,, W) sono le tracce integrali elettronica e fotonica: 


[ E, | 5 Cd) 20 D(E, 1) p(E, 0, 1) — q(B, 0, 1) 1 
AH By = | no Oa, y ee 
De Dem ao nn e BDACD) 
(2) | 7 i 1 
p E g( W, 1) E 
G(E —. | RE aw CA = 
em = | | | Ra) 


Nella valutazione delle L-trasformate si sono usate le espressioni complete 
di Rossi-Greisen (cfr. (5); formule (2.37)) allo scopo di tenere conto dell’effetto 
dei poli sull’asse reale negativo del piano (A, s), conseguenza del processo mol- 
tiplicativo. 

Inoltre, onde evitare incontrollate approssimazioni ed in vista della neces- 
sità di calcolare derivate sino al secondo ordine rispetto al parametro di tra- 
sformazione À, si & evitato di esprimere la funzione B.(Ey, E) (cfr. (5), pag. 7) 
tramite la funzione ®(H, s) scrivendola, tenuto conto della correzione di polo, 


nella forma: 


6+ ico 


1 F LS (2,0, aye LUE as a 
B;(E,, E, A 1 2e | 2b) A | [A — A,(s)][A — 2(s)] 


X =) ER E(B, A, 3 = 0) ds ? 
€) \No 


6—ico 


—g(W, s)dW ; 


"(m iw y(W) + x(W) 
ze Ao( W, s) 


(5) B. Rossı e K. GREISEN: Rev. Mod. Phys., 13, 241 (1941). 
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4 A SRE 
ea) (pee 9 Fei FS IE 
ee s+1le ai, 


le: ih; ==. 


Applicando il teorema della Faltung e servendosi di elementari proprieta 
della trasformazione di Laplace, le (1) si possono scrivere nella forma: 


| 2 RNS 2 n| F. 1 A On A dE' 
| GE ies = ERE (Ky, E) A | T(E, E ; À. s) Dar LB , E, An r) B® pe 
(3) | fo D 


2 x an’ 
COR" y= EG (Er W) 24 2 a 


| Iw 
"o 


L’integrale sull’energia EH’, calcolato dopo aver esplicitato le L-trasformate degli 
spettri elettronico e fotonico di bassa energia, fornisce al successivo integrale 
in campo complesso un polo nel punto 7 = s 42, Applicando il teorema dei 
residui si ottiene: 


= E, 2 4 BE, 2 
(4a) CD | ) >, 1H) = es fx(E) 5 B<e', 
Paz P 
, E\22 E.\: f Me 
(4b) OR" y= (5) > JOUE =) px(W), WE, 
1 
in cui: 
(5a) TE) = 
5 DE, = (a eo Ay ek? RES 
FC na (L)F (Bo) E) | £! ae E) DAN [! re » ] Yla=o; a eho) 9, 
=e =i) mit as Wyn = = B)f(e +E)P 7; == 12 


(bb) Oe Wee 


hi EH (1)ox B 2 No — Hl a OUT, s+ Vie EN 
ea ejay) () e +W Site ES Me 


di =o ; pa EN LOUE (no — We WANNE on — dr 


| TES E', A, Ss) JA L;(E', W, A, r) EB” hese 


nl 
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Le funzioni che compaiono nelle (5a, b) sono espresse da 


oe re en 
il DE GED Ale EDS A DN)” 
O(s +2)f(H, 90,8 +2) 
Eco Noe 4 Dr eee 
AR emery ce ner 
x Dios EU A Se ee 
ER fete Du rs one 
2(5, aie A)Q(E, À Saale 2) 
0) oe iR NS: Eee Ne si == a oe 
ee 212 en OL 
; 1 D, 
P(E, A, s #2)=(5 +20) n De: en ee, 
4 1 pie 
QUE, 2, 8 +2) = (5 +2) 3.18 
(6) | . 
A Ÿ W all 112 Mt = « 
dE" we À ee en (E, A, 3 re), 
x exp 4 _ ) pv E)| dE" + q(#, À, s +2) ; 
Cs + 2) 


sont) tte) 


"g(W, À, s +2) 3 
AN: 
We 


a WM + xW) ur ie] 
| Pa-o(B, A, 8 +2) = puns s +2) g(W, 8 +2) 11 W2-0 dW 
5 


Le funzioni p(B, A, s+2) e q(E,A, s +2), che tiene conto della ionizza- 
zione, corrispondono alle (7) (8) del lavoro citato (4), nelle quali perd non com- 
pare la correzione di polo, inessenziale nella teoria verticale. Le (5) sono valide 
nella regione delle basse energie (E<e, W<e) in cui le perdite radiative 
sono trascurabili rispetto a quelle per ionizzazione. 

Nella regione delle medie energie, che supporremo limitata fra &! e No, COT- 
rispondenti ai limiti per le sezioni d’urto asintotiche di formazione di coppie 
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e di bremsstrahlung, le precedenti funzioni si modificano in quanto la fun- 
zione g(W, s +2) si muta nella funzione C(s + 2)/W e, come conseguenza 
della sensibile diminuzione del termine correttivo sulla sezione d’urto totale 
coppie + Compton, la funzione —A,(W,s-+2) tende al valore asintotico 
dy — 0,773. Indicando con IW(E) e J(W) i nuovi valori che si ottengono 
per le (5) dopo aver posto y: —1 (k — 1,2, 3,4), ed w, = 1 (k —1,2) e 
sostituito 8 con Æ e, rispettivamente, con W, le (4) divengono nella regione 
delle medie energie: 


E. 2 £ = + FE E, . Man 
(Ta) E22" RX = (= 2, cu i) = (2) (=) ; © <sE<m 
— E, DRE : _ E. 2 a 
(7b) Kol — ( ) AZ OM) ( 5 E<W<n- 
Ww) oo W 


| | DER mine E>r 
Bo. \n Sa 
= ER 
(Sb) (©) Li? =) 3C(1) Dar ; 


Y = 1 Be (À z 18.0 W> e 
Me Fa RE RENE = ae ae? 


in quanto al limite per X e W tendente all’energia 7 si ha I®(E) = JW) —() 
per 41 e le tracce integrali (2) si identificano con le traece calcolate in 
approssimazione « A » (5). 

Nella Tabella I sono confrontati i valori dell’angolo e della deviazione qua- 
dratica, mediati suilo sviluppo longitudinale dello sciame, calcolati secondo le 
formule (4) (7) (8), con quelli ottenuti da ROBERG e NORDHEIM (7). Il valore 
BP = 86 MeV assunto da questi Autori invece di 6 = 98 MeV, com’e stato 


(°) A. BORSELLINO: Nuovo Cimento, 6, 543 (1949); 7, 323 (1950); 7, 638 (1950). 
(7) J. ROBERG e I. W. NORDHEIM: Phys. Rev., 75, 444 (1949). 
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TABELLA I. 
ROBERG e NORDHEIM Formule (4a) (7a) (Sa) 
FG ee WBS Fee en FE 
207° CR} (Oe (Rs 
co (E,/E)? 0,182 (2,/E} 0,214 (E,/E)? 0,1896 (E.[E)* 0,2417 
10 0,16 0,18 
7 0,15 0,17 
5 0,14 0,16 0,14 0,17 
3 0,13 0,140 0,115 0,15 
2 0,11 0,125 0,11 0,13 
175 0,095 0,105 0,10 0,11 
1 (2,/B) 0,085 (#,/B)? 0,085 (E./ß)? 0,08 (£,/B) 0,08 
0,75 0,135 0,120 0,13 0,10 
0,5 0,22 0,165 0,21 0,11 
0,4 0,28 0,190 0,25 0,20 
0,3 0,35 0,235 0,30 0,27 
0,2 0,57 0,31 0,46 0,37 
0,15 0,80 0,48 0,64 0,44 
0,10 1,23 0,46 0,88 0,52 
0,05 2,3 0,64 1,25 | 0,57 
ROBERG e NORDHEIM Formule (4b) (7b) (Sb) 
EJB WJB - == 
K OF dw x > De x OP} € Re 
co (E./W)? 0,060 (#,/W)? 0,38 (L/W)? 0,0586 (L/W)? 0,2508 
10 0,054 0,32 
m 0,050 0,285 
5 0,047 0,260 0,042 0,22 
3 0,040 0,220 0,036 0,20 
2 0,036 0,195 0.038 0,18 
1.5 0,034 0,175 0,030 0,17 
1 (E,/ß)? 0,029 (E,/ß)? 0,155 (E./ß)? 0,025 (#,/B)? 0,125 
0,75 0,043 0,255 0,035 0,23 
0,5 0,062 0,42 0,050 0,37 
0,4 0,076 | 0,53 0,065 0,47 
0,3 0,100 0,72 0,085 0,62 
0,2 0,145 1,02 0,125 0,92 
0,15 0,185 1,30 0,155 1,20 
0,10 0,26 1,80 0,205 1,68 


0,05 0,40 2,70 0,275 2,09 
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assunto nei nostri calcoli, influisce in modo poco rilevante sui valori relativi 
alla componente elettronica e in modo piuttosto sensibile su quelli della com- 
ponente fotonica. Per valori del rapporto E/ß e WIB maggiori di 6,4 le pre- 
cedenti relazioni forniscono i valori asintotici dello sparpagliamento laterale. 
Questo fatto, che è diretta conseguenza delle condizioni imposte nella teoria 
verticale, impedisce di descrivere il comportamento laterale dello sciame nella 
parte bassa delle alte energie. 


Ringraziamo i laureandi B. DEQUAL e L. TAFFARA che hanno eseguito i 
lunghi calcoli numerici. 


SUMMARY 


The integral spread, averaged over the longitudinal extension of an air shower, 
is caleulated according with the unidimensional theory of the electrophotonic cascade 
developed by the Padua group. The expressions thus derived, which give numerical 
results in accordance with ROBERG and NorDHEIM’s ones, describe the lateral spread 
in the low and medium energy region. In the high energy region they give the relations 
stated by BORSELLINO. 
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Über die Ausbreitung sehr langer elektrischer Wellen 
um die Erde und die Signale des Blitzes. 


W. O. SCHUMANN 


Elektrophysikalisches Institut der Technischen Hochschule - München 


(ricevuto il 5 Ottobre 1952) 


Zusammenfassung. — Es wird mit Hilfe der Methode der singulären 
Eigenfunktionen eine Formel für die Ausbreitung langer elektrischer 
Wellen um die Erde mit Berücksichtigung der Ionosphäre aufgestellt. 
Mit Hilfe der gewonnenen Formel wird die Ausbreitung elektromagne- 
tischer Signale einer Blitzentladung durch Fonrierintegration gewonnen, 
indem für den zeitlichen Verlauf einr Blitzentladung einige einfache, 
prinzipiell verschiedene Zeitabhängigkeiten angenommen werden. 


1. — Ausbreitung dieser Wellen um die Erde. 


Schon seit längerer Zeit besteht das Bedürfnis, für die Wellenausbreitung 
eines Dipols auf große Entfernungen eine einfache, leicht übersehbare Formel 
zu haben, die exakt begründet ist. Dies läßt sich in einfacher Form nur für 
lange elektrische Wellen durchführen, deren radiale Verteilung sich durch eine 
Taylorreihe angenähert darstellen läßt. Für diese Wellen ist der Ausbreitungs- 
raum praktisch die Luftschicht zwischen der Erde und der Unterkante der D 
Schicht der Ionosphäre in ca 60-70 km Höhe. Die D Schicht selbst verhält 
sich wie ein Leiter von der Leitfähigkeit x, etwa wie von Süßwasser. Da die 
Dämpfung der Wellen in der D Schicht weitaus größer ist als an der Erdober- 
fläche, kann man das Problem so idealisieren, daß eine Kugel im Radius 
ry =a mit dem Leitwert unendlich, umgeben ist von einer konzentrischen 
Luftschicht und diese wiederum für r = b begrenzt ist durch eine homogene 
Schicht eines Leiters. Schon G. N. WATSON (!) hat in seiner klassischen Arbeit 


(1) G. N. WATSON: Proc. Roy. Soc. London., A 95, 546 (1919); s. a. W. O. ScHu- 
MANN: Zeits. f. ang. Phys., B. 4, Heft 12 (1952) wo diese Methode benützt ist. 
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dieses Problem ganz allgemein behandelt und es durch seine berühmte Trans- 
formation gelöst. Wir wollen aber keinen Gebrauch davon machen, sondern 
eine neuere Methode anwenden, die meines Wissens zuerst von H. G. BOOKER 
und W. WALKINSHAW (?) angewendet wurde und unabhängig davon von A. 
SOMMERFELD #) als Methode der « singulären » Eigenfunktionen beschrieben 
wurde. Diese Methode ist weniger abstrakt und dem Denken des H. F. Techni- 
kers näherliegend. 

Setzen wir die Maxwell’schen Gleichungen für ein sphärisches Koordina- 
tionssystem an (7, 9,y) und nehmen keine Abhängigkeit von m an, so können 
wir für die Feldkomponenten E,, EB, und H,, die uns speziell für den verti- 
kalen Dipol interessieren, mit dem zeitlichen Gesetz exp [jot] schreiben: 


(1) Ek, = zone (sin 9: pe ‚Eu = = a RHONE a 
r Ve sin 4 20 - 20° 


wobei die Funktion /7T die dreidimensionale Wellengleichung 


(0) 


(2) AVAL SOE == ; = no ov el ’ 


zu erfüllen hat. Besitzt das Medium Leitfähigkeit x, so ist für € der Wert 
— jx/m zu Setzen. Die Funktion JJ hat die Grenzbedingungen HF, = 0, bei 
y — a, Erdoberfläche, und E, und Hy stetig bei r = b, Grenze Luft-Jonosphäre 
zu erfüllen. Es gilt also 


wil 
[ r—a@ 5 == (Ne 
3 | 
(3) NM | 1 à \=(4 wll (Ven Ce 
mt Ve rh We à ) SN LEE : 


Es folgt daraus, daß (/7 dr/Z/dr) an der Grenze r = b stetig ist, was für 
später wichtig ist. Befindet sich an der Stelle 9 = 0, r =r, ein senkrecht, 
radial gerichteter Dipol, so ist dessen Potentialfunktion, ohne Störung durch 
die Erde und die Ionosphäre durch 


Jl exp re jkR | 


(4) I, = —— 
Ir, Ve R 


@) H: G. Booker und W. WALKINSHAW: Conference of the Phys. Soc. and of the 
Royal Meteorological Soc. 8th April 1946, Appendix IT (beschrieben in H. BREMMER: 
Terrestrial Radio Waves, (1949), S. 202). 

. 3) A. SOMMERFELD: Partielle Differentialgleichungen (Leipzig, 1947), S. 216. 
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gegeben, wobei J der Strom, ! die Länge des Dipols und À der Abstand des 
Punktes vom Ort des Dipols ist. Die Funktion // muss also nahe am Ort des 
Dipols in die durch (4) vorgeschriebene Singularität übergehen. 

Die Lösung der Differentialgleichung (2) führt auf Kugelfunktionen P,(cos 9) 
und Bessel’sche Funktionen 


IT 771.2 
Er — y+ 12 (0), 0 — ki 
20 


Die Potentialfunktion ist also, abgesehen von einem später zu bestimmenden 
konstanten Faktor im Luftraum gegeben durch 


[&,(kır) sin Gé, (Arr) | -P,(cos 0) = x,(kır) P,(cos 0) , 


— 
On 


und in der Ionosphäre durch 
(6) GE(kr)P,(cos 0), ky = Vaux: exp[— ja/4] . 


Die Konstanten €, und e, werden durch die Grenzbedingungen (3) bestimmt. 
Die Grenzbedingungen bestimmen nun nach der Methode der singulären 
Higenfunktionen die Ordnung » der Funktionen durch die transzendente Glei- 


chung 


/ 2/ / 2 , kb # 2 + n 
(7) 7» (0)9» (Oa) — Hs (9) (Ou) — N En [m(e)1% (Oa) — 43 Ov) m (Ox) |= 0. 


TO 171,2 


Oo = Mb, «a = Ha, me) = e&"*(0) = | > Hr 12(0)- 


WATSON hat in (!) gezeigt, daß der erste Teil der Gleichung (7), der bei 
unendlich gut leitender Ionosphäre (x — co) allein auftritt, unendlich viele, 
rein imaginäre Wurzeln hat, und eine endliche Zahl reeller Wurzeln zwischen 
den Grenzen + k,b. Das zweite Glied, das nur bei leitender Jonosphäre auf- 
tritt, ergibt zusätzliche komplexe Glieder, die zu den oben erwähnten Wurzeln 
hinzuzufügen sind. 

Auf Grund der Tatsache, daß die Luftschichtdicke b—a & 70 km sehr 
klein gegen den Erdradius a = 6350 km ist, versuchen wir für lange Wellen 


(A in der Größenordnung von ~ 70 km) eine Näherungslösung zu erhalten, 
indem wir 7in eine Taylorreihe bis zur 4. Ableitung entwickeln: 


i id A a Ui A > nn A ; 
(3) le + Ae) = nle) + ne) de + n'(@) = +n"(e): a. 157 Obi 
EN RMON np: rn 
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Die hôheren Differentialquotienten bestimmen wir aus der Differential- 
gleichung für 7 


(9) (0) + 


(10) milo) (0) — m(o)m (0) =— 2, 


(11) 1 (On) (a) — Mv (Ov)Mv (Ga) = 


ae .; | vw + 1) en VE ne vv 1), yy + 1) 1 a 
= N il aie Ao- — Ao | o! ! (1 | sl Ae fe 


Für eine vollkommen leitende Ionosphäre, ~ — ©, ist dieser Ausdruck 
nach Gleichung (7) gleich Null zu setzen und ergibt den dann in Frage kom- 
menden Wert von »=»,. Wie man sieht, genügt (% + 1)/e = 1 für alle 
Ao dieser Gleichung, wenn berücksichtigt wird, daß Ao/o <1 ist. Es ist also 
in dieser Näherung (% + 1) = kla?, zu setzen und », ist eine große Zahl, da 


2700 40 000 
ie = —— =— 
4 Akm 


7 

(2 Wellenlänge), selbst für sehr lange Wellen. 
Wendet man dasselbe Verfahren auf das zweite Glied der Gleichung (7) an, 

so entsteht 

vy + nee y(v + 1) Ao® 


2 ~ Ot au 


2 


(12) nor)? (eu) — mg)» (Qa) = 2} 1 — f 


0 
S 


bei dem der Wert der Klammer sehr nahe an 1 liegt, so lange Ao?/3e klein 
gegen 1 ist, d.h. sobald die Wellenlänge 1 > 2xHVH/|3a ist, H= b— a Höhe 
der Luftschicht, a Radius der Erde; für a/H ~ 100 muß A> 0,35 H sein. 

Aus der Gleichung (7) wird also angenähert, wenn man voraussetzt, daß 
sich »(y + 1)/e2 nur sehr wenig von 1 unterscheidet: 


S 


| re LS 
13 Al en een. 
go) 1 0! ky né (lead) 


Da |k,/k,| <1 ist, bedeutet das 2. Glied nur eine kleine Korrektuy von ». 
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Da 

2/ 2 2 

Ny (0) v »-1/2(0) H,,(@) : 
14 See — = SL ee und =, —— =— dar ks) ea 
zn) 0 "'H},17(6) H,.1(0) A ï AE 
folet 

Gy = il) 1 ( v ; j 

15 es Cees pele 
2) 08 eg 


und hieraus mit 


[Bie Vox exp [— ja/4| und Ga 1, 


L 1 a .a Vio 
(16) v | D Oa (1 7 — = | = N) ET oO 32 . 
oA 2HV joux € 


Slt 


Es folgt hieraus eine Dämpfung der Wellen von der Größe 


| Peat ie 

exp (— He | = ; 
die schon Watson abgeleitet hat ((1) S. 1.) und eine Eindringtiefe, Abfall der 
Amplitude auf 1/e, in die Ionosphäre bei x = 10" 1/Qm, H— 70 km, 
@ = 6370 km, für À = 30 km, von ~1/2 km, bei f = 10 Hz, A = 30000 km 
von ~ 15 km. 

Da die Ordnung » eine komplexe und nicht ganze Zahl wird, enthält die 
Kugelfunktion P,(cos 0) bei 0 =z eine singuläre Stelle, wo sie logarithmisch 
unendlich wird. Die Funktionen 7,(kr)P,(cos 0) und &(k,r)P,(cos 9) nach 
Gl. (5) und (6) werden deshalb von SOMMERFELD als « singuläre » Eigenfunk- 
tionen bezeichnet. Physikalisch bedeutet das, daß an der singulären Stelle sich 
eine lokal konzentrierte Strahlungsquelle befindet, und wir werden im Fol- 
genden unsere Potentialfunktion dem singulären Charakter des Hertz’schen 
Dipols nach Gl. (4) für À — 0 anzupassen haben. Für Kugelfunktionen hoher 
Ordnung gilt die Näherungsformel 


/ 2 1 7 
P,(cos 0) = |} cos Ho Lee : OO gah 
zy sind 2 y 


Für komplexes » mit positiv reellem Teil wie hier ist die Kugelfunktion 
für 0 + 0 regular, P,(+ 1) +1 und divergiert für 0 > x als 
sin var 


(Ta) B,(c080) = 


cos 0 — 1 


log (x — 0} ... 


(4) sa. W. O. SCHUMANN: Zeüs. f. Naturfor., Ta, 15 (1952). 
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Nehmen wir unsere Strahlungsquelle bei d=0 an, so haben wir statt 
0 7x0 zu setzen und erhalten als Darstellung 


[> 5) (x 0) 4 


Ist » komplex, so wird mit cos? = (1/2) (exp [je] + exp — el) 


Os 


ere 
P,— COS 0) = Van C 


er! 1 ; 
as Pen) lee - 


ee ily +5 (a — 9) 
2mwjsind DE } 


wenn die Exponentialfunktion mit negativ imaginärem Teil als unbedeutend 
weggelassen wird. 
Die Potentialfunktion // wird nun in der Form 


(19) T(r, 6) = > D,y,(kır) P,(— cos 8) , 

angesetzt, wobei im Luftraum 

(20) pr) = LT) = Er) + E(k?) , 

(s. Gl. (5)) und in der Tonosphare 

(21) (kor) = G&rlkır) 5 

bedeutet. In unserem vereinfachtem Fall reduziert sich die Summe auf ein 
einziges Glied, da sich nur ein Wert von » aus Gl. (7) bzw. Gl. (13) ergibt. 
Die Faktoren D, können bestimmt werden, da die Funktionen y, für verschie- 


dene Indizes » zu einander orthogonal sind. Dies folet aus der Differential- 
gleichung der y Funktion 


2{y 
r d (rp) — [k? 
dr? 


(22) r2— vv +1), y, = 0, 


die sowohl im Luftraum als auch in der Ionosphare gilt. Schreibt man die 
Differentialgleichungen fur die Funktion y, und y, hin, multipliziert die erste 
mit y, und die zweite mit y, und subtrahiert beide Ausdriicke und integriert 
so entsteht 


co 
MY) gy Ou) | 
dr Poil 


a 


co 


(23) Po +1) — u(u+1)] | yyy dr =| rv, 
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Diese Integration ist von + =a bis r = co durchzuführen. Dabei wird 
auch die Grenze 7 = b durchschnitten. An dieser Grenze ist aber nach Gl. (3) 
ry, d(rp,)/dr stetig, und der Ausdruck in Gl. (23) ebenfalls. Er trägt also 
zum Integral nichts bei. Da nun für r = a, d(ry)/dr verschwindet, und für 
y= oo wegen des asymptotialen Verlaufs der Hankelfunktion der ganze Aus- 
druck verschwindet, ist 


co 


| vpn Ole =O, 


0 


co 
Um die Normierungsintegrale | y, dr zu finden, bestimmt man nach Gl. (23) 


a 


co 
co d(ry) d(rpu) | 
7 BEL, : Wr ein di 
N, = lim | Yup, dr = lim | dr dr E 
2% ge [of +1) — u(u +1)|| 


indem man Zähler und Nenner nach « differenziert und dann zur Grenze 
übergeht. Man erhält dann 


Um diesen Ausdruck auszuwerten, ist es besser, von den Funktionen %,. 
zu den Funktionen ®, = oy,(o) überzugehen, d.h. von den Funktionen 


2 | IT 
&,(0) = | on Ja Poe 1/2(0) , 


zu den Funktionen 


re 
TO 
N, = = H, .1)2(0) , 
wie in Gl. (7). Dann wird 

CO 
‚db d | 
A ci CE | 
iN ug |* 
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Die Differentialgleichung dieser Funktion ist wie in (9) 
’ | 1 | 
D" + 1-7 À P—0. 


Um die Rechnung zu vereinfachen, vernachlässigen wir den kleinen Beitrag, 
den die Ionosphäre zu dem Normalisierungsintegral liefert, d.h. wir integrieren 
> = b statt bis r = oo, d.h. nur innerhalb der Luftschicht. 


nur von 7 = a bis 
Entwickeln wir auch hier wieder die Funktion ® 


Ag? Aw? 
Dior) = B(Q.) + D'(o) Ae + Bla) + PUR): TG, 
! I 7 à Ao? 
D'(or) = B'(0a) + B"(oa)* A9 + Pie)» 
TRE, (1 vw + 1) 0 
dv 0° | 0° dv 
Dee g (1 Ads =) 0, 
0° \ 0? 
und 
oo" 2p(v +1) 08 , w+) {y + 1)\ 28° 
a ee ol 
ov 0° 0 dv 0? 0? dv 
und setzen 
(1 x v(v se =) 
Où 
als sehr klein voraus, d.h. wir setzen 
à wD, 2» +1 ‘5 vy + 1 
Pa, er Das pee ya que 
oy Oa 0 
d®, 2 My + 1) 0D, M +1, 
= = 2y + eh a D. 
dv 5 | med tI Ps 0° wy 0? 


= 


und setzen schließlieh ©, = 0 wegen der Grenzbedingung dr/1/dr = 0 an der 


Erdoberfläche, r — a, so erhalten wir 
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Mit diesem Wert von N bestimmen wir nun den Potentialausdruck der 


Gl. (19) 


co CO 

| I(r, 6)y, (kr) dr = D | yi(kr) dr- P,(— cos 0) = D,N,P,(— cos 0) . 
(5) Um schließlich die Anpassung unseres Potentialausdrucks an die Singula- 
rität des Dipols bei 0 = 0, r = r, nach Gl. (4) 


Ji  exp[—jkR 
Il, R ee ] 


Arr © V & 


vorzunehmen, setzen wir in unmittelbarer Nähe von R—0, wo 77 in IL, 


übergehen muß 


1 oo 
athe 5 v cr) dr = 
A à P,( — cos 6) | TQ O)y (kr) j 
a : “exp [— jkR] 
i € j : »,(kr) dr, 
Arr Ve P,(— COS 0) : R y 
0 —0, R=0; oe 


=; — Ir, COS Os r, radiale Koordinate des Dipols. 


m 71a) GUS == HOU) pr (kr) > y, (ro), 


R=nVm +0, exp[—jkR] >1, 


ergibt sich, indem man den Integrationsbereich auf die sehr kleinen I ntervalle 


von ro(1—) bis (1 + 7) beschränkt, mit Hilfe von Gl. (7a): 


Jl 
2 k, To) 


D,N, = — —eyıl A: 
D # PE NN Sin vz 


mit ry =a, wenn der Sender auf der Erdoberfläche steht. Es ergibt sich 
(5) Genaueres s. A. SOMMERFELD: Partielle Dijjerentialgleichungen (Leipzig, 1947) 


S. 218. 
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schließlich das Hertz’sche Potential auf der Erdoberfläche für r = a 
IU eK) ML 
DD EC 3, her P(— cos 
Hua 4awVe I, Snyx ne 


und da 
y ® (kya ) 


JI 1 1 
P(— cos 6) . 


I — 5 = / = u . 
dawV € H(t ii el SIN VT 
\ Q 


2 exp [jyz] , 


Nun ist 
j 
2 


: iV > F 
sin yt = 2 (exp [pr] — exp [— jra]) & — 


wenn » einen negativ imaginären Teil hat. Daher ist mit Gl. 18 


P,(— cos 9) |. ne Re tte. 
——— S— | — xp [—! 2 
sin vz jay sin 0 IE Ne IE 


und damit wird 
Ail! 1 


ik. TA. = © 
4dawVe H 


k, und a sin 0 = D, D Entfernung am Erdumfang gemessen, 


OS 
i ale 9 
| jay sin 0 ep I? + 1/29], 


und mit » & 
EN LE V. 2 1 ee 
Do 2 HET TO. sn 0 pP ) [4 . 


Die radiale Feldstärke auf der Erdoberflache ist nach Gl. (1) 
j 080 wl 
ee 
AN Gases 120 00 


und nach der Differentialgleichung der Kugelfunktionen ist dies 


74 - Il Nuovo Cimento, 
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so daß 


-exp [— jv + 1/2) 


N A 1 Vis ee 


aceD ee sin 4 


Wenn | = 2h gesetzt wird, h Antennenhöhe oberhalb der Erde und 


VE = 1207, D = Dou cs - 


€ 
wird 
Dr | 
D De NG —3 / _ Sch PRA 4 - eX A DI, , 
4) Bu = 61097 Vio | — 5 exp Ei + 129, 


und wenn alle Längen in km gemessen werden 


3x .h Vio 19 
OF I [3 = ) 
(25) Ea — 0,2 | ae ae VD ye exp [— jv + 1/2) 6], 


(hk, H und D in km), 
wobei nach Gl. 16 


: N | D) 1 1/50 D 
exp [— j(v + 1/2) 0] = exp - 105 ony hr . 


Wünscht man die Feldstärke abhängig von der Strahlungsleistung eines Dipols 


über unendlich gut leitender Erde, so ist diese 
Te + h\? 2 
We = 1607: If A 10! > HP. 


Dies gibt mit Gl. 24 


SE ] or ee 1 / 
IT [62] 
26) Er. 191 I VNwWor— V- _ exp |— | = 
20) “ref 6 EWHVD Vo sind FE 2Hc V 274 2 ne 
und wenn H und D in km gemessen werden 
i = 
2 : | re Ar. il a at i o D 
2) 3, = 900 / vw, == Ve ex V1 vim, 
ref 15 Ÿ Niw Hu VD Vol no p 2¢V Qux H /m, 
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oder in Wellenlängen À ausgedrückt 


ul [Ric / 0 ia 1 1 
(28) Er = 300 VN cw I —— exp | — = ——. 
ae se Dim | sin 0 400:V3 Vi H 


mV/m, 


(x in.1/Qm, À in km). 


(s. a. W. O. ScuumMANN: Zeits. f. ang. Phys., Bd. 4, Heft 12 (1952) wo diese Formel auf 
anderem Wege abgeleitet ist). 


2. — Diseussion der Formeln und Ausbreitung einer Blitzentladung. 


a) Frequenzabhängigkeit. — Vergleicht man die Ausbreitungsformeln 
Gl. (27) und (28) mit schon bekannten Formeln, so ist besonders für diese sehr 
langen Wellen charakteristisch das Auftreten der Frequenz bzw. der Wellen- 
länge nicht nur im Dämpfungsfaktor, sondern auch im Amplitudenfaktor. 
Es bedeutet dies eine wesentliche Verbesserung des Empfangs mit wachsender 
Wellenlänge bei gegebener Strahlungsleistung. Der Grund scheint darin zu 
liegen, daß bei diesen großen Wellen auch schon des Nahfeld der Antenne 
durch die Ionosphäre so beeinflußt wird, daß die abgestrahlte Energie mit Zu- 
nehmender Wellenlänge wächst. Vergleicht man diese Ausbreitungsformel mit 
den Formeln von AUSTIN (*), von ESPENSCHTED, ANDERSON und BAILEY (?), 
und von ZINKE ($), so ergibt sich mit x= 10° S/m, H= 70 km, D= 10000km, 
0 — x/2 nach der neuen Formel mit 


À / D 
Ben = 4,3 7 | VX. xp |— 0,0021 mV/m 
te <p a pl VA Lire 
für 
est) 300 3000 30000 km 
F, 
f= G30 0,284 2,14 8,23 mV/m. 
VN, 
ach AUSTIN 
i 300 x | 6 0.0014 D 
= = Ne) Ze EXP I ner m 
af Sapam, ae Wera = 
E NS axe 
Oa Orla: 0,0239 mV/m. 
VN 
(6) Le W. AUSTIN: Proc. Inst. Rad. Engineers, 20, 689 (1932). 
(7) L. ESPENSCHIED, C. N. ANDERSON und A. Baicex: Proc. Inst. Rad. Eng., 14, 
7 (1926 


)- 
(8) O. ZINKE: Frequenz, 1, 16 (1947). 
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Nach ESPENSCHIED, ANDERSON und BAILEY, 


a 
i= Dale 
E A idler À 
—— 0), Olls ay 0,0357 mV/m . 
VN 
A = 30 300 km - 
Nach ZINKE 
15 So 
ip = VN pls -exp [— 0,003 D/A] 15 = 60-V5.— 
VD | sin 0 ‘ VH 


E 
a ET 1,89 1,89 mV/m. 


A =30 300 3.000 30000 km 


Die Ubereinstimmung der neuen Formel ist am besten fiir 2 = 30 km mit 
der Austin’schen Formel. Die Zink’sche Formel gibt fiir 2 = 30 km viel 
größere Werte, stimmt aber mit der neuen Formel für 2 = 1000-2000 km 
etwa überein. Es ist zu berücksichtigen, daß Austins Formel für Wellen- 
längen 2 = 300-25 000 m aufgestellt ist, die Espenschiedformel für À = 5 250 
17500 m, während die neue Formel etwa von A = 30 km an eilt. Die Zink’sche 
Formel ist für À = 2-20 km gedacht. Das wesentliche der neuen Formel ist 
die Steigerung der Empfangsfeldstärke um mehr als das 1000 fache, wenn 
die Wellenlänge von 30 km auf 30000 km steigt, bei gleicher Strahlungs- 
leistung der Antenne. Übereinstimmend mit ZINKE ist, daß die Entfernung D 
nur mit 1/VD in die Formel eingeht, in Gegensatz zu AUSTIN und ESPEN- 
SCHIED, WO 1/D auftritt. Dagegen tritt in der Zink’schen Formel die Höhe 
der Luftschicht H mit 1/VH auf, während sie in der neuen Formel mit 
1/H erscheint. Die Abhängigkeit des Dämpfungsfaktors von der Frequenz 
ist sehr stark durch den Zustand der Ionosphäre bedingt. Da wir der 
Ionosphäre einen konstanten Leitwert zuschreiben, tritt, wie üblich, in 
solchen Fällen 1/VA im Dämpfungsfaktor auf. Das setzt voraus, daß die 
Schwingungsfrequenz wm sehr klein gegen die Stoßfrequenz der Elektronen im 
Plasma ist, die ihrerseits dem Druck proportional ist, also mit wachsender 
Hohe H abnimmt. ZINKE nimmt im Gebiet zwischen 2 = 2 und 20 km bei 
À — 6 km die Stoßfrequenz gleich der aufgeprägten Frequenz an und gelangt 
dadurch zu einem Faktor 1/2 im Dämpfungsfaktor. Während bei Wellen- 
längen größer als 20 km auch bei ihm ein frequenzunabhängiger Leitwert 
zustande kommt, d.h. dort auch 1/V À im Dämpfungsfaktor auftritt. Ver- 
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glichen mit den exakten Berechnungen BREMMERS (*) auf S. 229, Fig. 65, seines 
Buches mit x — 4-10-4 S/m, H= 100 km, À = 10 km, D = 10000 km, ergibt 
sich nach der neuen Formel die Feldstärke von 12,2 uV/m gegen das dort 
bestimmte 2,5 uV/m. In diesem Bereich ist die Formel offenbar nicht mehr 
brauchbar, da A/H= 1/10 ist, während der Gültigkeitsbereich der Formel bis 
etwa A/H = 1/2 geht. 

Anders sind die Frequenzabhängigkeiten, wenn der Strom des Dipols ge- 
geben ist, Gl. 25. Dann ist die radiale Feldstärke 


Re 1 4 0. D 1 
Saga Vo "exp | 77 =a mV/m 
A Hymn Ne sin 4 = H 400V3 V rm IV] 
bzw. 
29) re Ws 7 eg ya 
| ue = 0,2 | — É © ayy ‘him | al ~ EXP - 5 V: 2 | ‘ 
9 Hi Die sin 9 Den x V2 
bzw 


exp |— 


a Ip /2 DE 
12 er HE 


Es ergibt sich hieraus für jede Entfernung D eine günstigste Frequenz 
Mopt bzw. eine günstigste Wellenlänge Ao», für maximales £,, 


die proportional mit dem Quadrat der Distanz wächst. 
So ergibt sich z.B. mit D — mH, H=70 km, und x— 107? S/m, 


m 10 30 40 50 100 
D 700 2100 2 800 3500 7000 km 
Aopt 2 18 32 50 200 km 


foot 150000 16700 9400 6000 1500 Hz 


Für andere Frequenzen gilt 


2 | ges Er à > 

FE, | 0) | / i) rn SH? x 
ea ie eu ’ OO On tp Fo en 

Bae Oopt Oopt D € 


(9) H. BREMMER: Terrestrial Radio Waves (Elsevier, 1949). 
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und man erhält für gegebenen konstanten Antennenstrom für verschiedene 
Frequenzen ein Spektrum nach Fig. 1. 
Das Spektrum fällt für &/oopt <1 rasch auf Null nach Formel E,/ By ax 
=}? V wl@opt und erstreckt sich sehr weit ins hochfrequente Gebiet. Je 
größer die Distanz D ist, um so mehr 
spielen die langen Wellen eine maßgebende 
Rolle im Spektrum. Besteht das Stromsignal 
aus verschiedenen Frequenzen, so wird bei 
kleineren Distanzen D, wo die optimale 
w 


Op Frequenz @opt ziemlich hoch liegt, das Ge- 
ie 1 9 


biet der niederen Frequenzen, der langen 
Wellen, stark benachteiligt sem. Umgekehrt 
ist es bei größeren Distanzen, wo die op- 


Fig. 1. — Ausbreitungsspektrum des 
Antennenstromes. 


timale Frequenz ziemlich tief liegt, und wo 
die niederen Frequenzen dominieren und die hohen relativ schwach vertre- 
ten sind. 


“ Blitzspektrum ” 
und Form des ausgesendeten Signals (Atmospherics). 


Die Zeichen, die eine Blitzentladung aussendet, hängen erstens von dem 
zeitlichen Verlauf der Blitzentladung ab und zweitens von der Verformung, 
die das Zeichen während seiner Wanderung erleidet. Faßt man den Blitz als 
einen erzwungenen Stromverlauf auf, so ist die spektrale Verteilung der Ver- 
formung durch Gl. (24) und (25) nämlich Vp exp [— «/2p1/?], p = jo, ge- 
geben (wenigstens für die Frequenzen für die diese Gl. gilt). Es sollen nun 
einige einfache zeitliche Formen in Bezug au ihren zeitlichen Verlauf besprochen 
werden. 


1. — Dirae Impuls des Stromes. 


Setzt man eine Blitzentladung als einen kurz dauernden rechteckigen Strom- 
impuls von Strom J und der kurzen Dauer 7, vom Moment Jr voraus, in der 
Grenze als einen Dirac-Impuls von diesem Moment, so hat sein Spektrum 
die Größe Jr für alle Frequenzen. Nach Durchlaufen einer Entfernung von 
D= 10000 km tritt das Maximum von E, für eine Frequenz von fopt = 705 Hz, 
Aopt = 425 km auf. 

Nimmt man eine Höhe der Blitzentladung h = 5000 m an, J, = 5-104 A 
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— 10-s, so folgt aus Gl. (29) Ey ax = 0,150 mV/m für À opt. Für 
flfont = 1/100, f = 7,05 Hz ist H = 0,25E, ,.; für fe 100,9 = 10500) Ey 
ist dagegen E, schon auf 1,24 -10=E, ax gefallen. 

Um den zeitlichen Verlauf des Impulses nach Durchlaufen der Strecke D 
zu finden, gehen wir von der spektralen Verteilung nach Gl. (24) und (25) 
aus, die sich mit p = jo in der Form 


D = 


ES m te RG sn 2 ; 
(30) 18). = WENN exp Ten vo, 


X lol 


darstellt. Das Fourierintegral über die Frequenz (!°) ergibt den zeitlichen 
Verlauf: 


IDEE 


60 = ‘= 
2 4H? x’ 


ete a 1 1) 
Gt 6 m=" — | = 1) —, exp (— u )» 2 = 
1 re 


ol? 94/7 \2% dx 


wobei die Zeit vom Beginn des Eintreffens des Signals an gerechnet ist. 
Fig. 2 zeigt den Wert von (BC) ae. 

Aus der Figur sieht man, daB der Höchstwert des Signals mit 1/D7? ab- 
nimmt, und daß es sich proportional zu 


D? verbreitert. Bemerkenswert ist, daß Er or 
h C 
es bel 3 
; 2 Ie NR 
NCAA 2 
durch Null geht, und daß dann eine nie- 1 


drige negative Halbwelle folgt, die sehr 


langsam mit der Zeit ¢ verschwindet. 


01 0,5 


Besteht der Blitz aus mehreren aufeinan- T 
derfolgenden Dirac-Impulsen, so erhält Fig.,2- =Ditau-Empuls. des Stromes 
man für jeden solchen Impuls ein Bild 

nach Art der Fig. 2. Diese Bilder sind entsprechend zeitlich relativ zueinander 


zu verschieben und dann zu überlagern. 


3. — Rechteekiger Impuls endlicher Länge und konstanter Höhe. 


Dieser Verlauf läßt sich darstellen durch die Überlagerung zweier « unit 
steps » Einheitsstufen-Funktionen, entgegengesetzten Vorzeichens (Fig. 3), von 


(2%) G. A. CAMPBELL und R. M. Foster: Fourier Integrals for Practical Applications, 
Bell. Tel. System Monograph B 584 (New York, 1931), S. 93, No. 806, wo alle Regeln 
u. Integrale angegeben sind. 


1132 W. O. SCHUMANN 


denen die eine zur Zeit t = 0 den Strom von 0 auf den Wert J springen läßt, 
und die zweite nach der Zeit T den Rücksprung auf den Strom Null durch- 
fübrt. Die «unit step» Funktion hat das Spek- 
trum 1/p. Dies gibt multipliziert mit dem Spektrum 
die Verteilungsfunktion 


1 xp [ 1/2 p1/2] Ds € 
— exp [— ad? mu = 2 
Vp I ! : 4H2 x 
ih . 
Mit 
Fig. 3. — Rechteckimpuls 
les Stromes. 1 oe 
peat 2 = 17.9.100 Ems under 
Or 
ist 
JDE 11 
i, = . 
H: 4,52-10? ’ 
also 


ihe Jah == 7D km 


1 | 
JD = él D OS D — 10 km 
4,5 | 
1 | me 
D=10H = ire SIR, D -700 km 
AD) 
il | 
D = 100H = a Ss D == 7D) Ladies 
45 


Aus der Spektralfunktion folet mit dem Fourierintegral der zeitliche 
Verlauf 


1 
= exp [— «/4t] . 
Vat 


Das Maximum 


1 À " œ 
Ym = — WUE Te = u —= = 
dl 


EE ve 


erreicht und es ist 


y | CUT À 
M — — Veex — r/2t]. 
Ym t x P l ; ] 
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Den Verlauf der Kurve zeigt Fig. 4. 


t 1 t 
TO ee 02 Ze 0 eo. 100: 
T Ym T Ym 


Wenn nach der Zeit T eine zweite Stufenfunktion gleicher Höhe und ent- 
vegengesetzten Vorzeichens einsetzt, so entsteht die Verteilung eines ursprüng- 
lich rechteckigen Impulses von der Breite 7. 
Fig. 5 gibt eine Skizze des Verlaufs fir den Fall 
P= 7247 und 107. 

Da t = x/2 quadratisch mit dem Abstand 
zunimmt, verschiebt sich das Bild der Fig. 5 


von großen T/t zu immer kleineren Tr während 
seiner Wanderung Fig. 4. — Stufenimpuls des 
A : = ce Stromes. 
Für kleine 7 sieht man den Ubergang zum 
Dirac-Impuls, Fig. 2. Fiir eine Reihe aufeinan- 
derfolgender Rechteckimpulse erhält man durch Uberlagerung mehrerer Bil- 
der nach Fig. 5 die zeitliche Verteilung. 


4. — Exponentiell abklingender Impuls. 


Verläuft der Strom nach einer plötzlich 
bei t= 0 einsetzenden abklingenden Expo- 
nentialfunktion exp [— ft], t > 0, so ist seine 
spektrale Verteilung durch 1/(p+ B) gegeben 
und die spektrale Verteilung des Signals 


Fig. 5. - Rechteckiger Stromimpuls. 


ee I N 
Da NEC durch 
De TOT p! 2exp[— ol? pre] 


p +6 


Die zugehôrige Zeitfunktion ist (CAMPBELL und FOSTER, 1. c., S. 40, No. 219) 


t 
exp [— Pt] | exp [Bt]@(t) at, 
0 
mit 
Gt) = | zZ ar exp — «/4t], 


SE 


(wie in Gl. (31) S. 1131 Funktion des Dirac Impulses Fig. 2). Für sehr kleine 
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Zeiten ft <1 verläuft sie wie die Funktion der »unit step» Funktion 


1 
= exp CAR 
Vat pl 


um aber mit wachsender Zeit rascher abzuklingen, durch Null zu gehen, und 
mit einer negativen Halbwelle zu enden. Für sehr kleine f# nähert sie sich 
der Form der Einheitsstufenfunktion (Fig. 4, S. 1133), für große 5 mehr der 
Form G(t) des Dirac-Impulses. Der Nulldurchgang liegt bei sehr großen P 
nahe bei t = «/2. Mit abnehmenden f, langsamerem Abklingen wächst 7 an, 
bis bei 6 +0, T— oo geht. Es ist also immer eine auf dem positiven An- 
fangsteil folgende, negative langsam abklingende Halbwelle vorhanden. 


5. — Ladungs- bzw. Nahfeldimpulse. 


Anders wird das Bild der Signale, wenn man sich bestimmte Ladungs- 
verläufe, bzw. Feldverläufe im Nahfeld des Dipols als gegeben denkt. Da 
Strom und Ladung nach der Gl. 


à = — dq/dt = — jog = — pq, 


zusammenhängen, wird jetzt die Verteilungsfunktion für die Feldstärke LH, 
nach Gl. (30) S 1131 


p?? exp [— al?pu?]. 


Verglichen mit der früheren Verteilungsfunktion 4/p exp [— 4/x pl], die 
das Spektrum Fig. 1 gab, ist die spektrale Verteilung jetzt ¢3(@/@ p¢)*?- 
-exp | — 3Vo/oopt] wobei @opt = 9(2/x) ist. also 9 mal größer wie in Fig. 1. 
Außerdem fällt die Intensität der Frequenzen für © < @opt Jetzt viel rascher 
ab, als in Fig. 1. Ausgedrückt durch 2 = V a«o/2, erhält man im ersten Fall 
Ylym = 62" exp [— +] und im zweiten Fall ¢3(2/3)? exp [— 2], wobei im ersten 
Fall das Maximum bei z = 1, im zweiten Fall bei z = 3 auftritt. Fig. 6 gibt 
eine vergleichende Skizze beider Verteilungen. 

Da « bei Entfernungen von D = 700-7000 km von (1/4,5) : (10° bis 10°?) s 
geht, liegen die Frequenzen in diesem Bereich bei f & 1,44 : (10° bis 10%)z* s7", 
Speziellfüre = 1 bei f= 1,44 - (10° bis 102)-Hz, bel z= 3 bet f= 1,3: (105-102) Hz. 
Setzt man die gleichen zeitlichen Verläufe, wie bisher für den Strom ange- 
nommen, jetzt fiir die Ladungen (bzw. die Nahfeldstärken des Blitzes) an, so 
sind die entsprechenden spektralen Verteilungen mit p zu multiplizieren. Nach 
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den Regeln der Fourierintegration (bzw. der Laplace-Transformation) heißt 
das, daß die neuen zeitlichen Verläufe als Differentialquotienten nach der Zeit 
der bisherigen Verläufe gewonnen werden. 


6. — Dirae-Einheitsimpuls der Ladung. 

(Entsprechend einer kurz dauernden positiven und einer unmittelbar darauf 
folgenden kurzen negativen Halbwelle 
des Stromes). Der zeitliche Verlauf des 


y 
Signals ist jetzt nach dem Spektrum 17 @ © 


pl exp [ap]; 


durch ee Ve 
2 


1 1 4 3 5 
eae a iy CD [—1/4@]- a 
2 mar? SX? Ties, ie = Ausbreitungsspektrum des 
1 Antennenstromes (1), der Antennen- 
- [1 — 122 + 12%], L = - ladung (2). 


1.) YsYm = exp [— 2] 


13 
2.) YsUn = «| | exp (j— 2) - 
gegeben. 
Die Kurve geht zweimal durch 
Null füræ =(1/2)+ V1/6 = 0,093 und 
0,305 entsprechend Fig. 6, und hat 2 


Maximaund ein Minimum. Fig. 7 stellt 


schematisch den Verlauf von & = yo?” 
Fig. 7. — Dirac-Impuls der Ladung. dar. 


7. — Rechteckiger Ladungsimpuls endlicher Lange mit konstanter Hohe. 


In diesem Fall erhalten wir für die Feldstärke mit der Einheitsstufen- 
funktion die spektrale Verteilung 


vp exp[— va ve], 
wie beim Dirac-Impuls des Stromes, entsprechend Fig. 2, S. 1131. Geht der 


Impuls nach der Zeit T wieder sprungweise auf Null zurück, so sind 2 Kurven 
nach Fig. 8 im zeitlichen Abstand T zu überlagern. 
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Da « proportional D? ist, nimmt 7/x mit wachsender Distanz ab, d.h. der 
Impuls geht allmählich von der Form 3 in 2 und dann 1 ete. über. Für sehr 
eroBe Abstände nähert er sich der Form des Dirac Impulses. Der Rechteck- 


Fig. 8. — Rechteckiger Ladungsimpuls. 
LT 005% 
a NEN € 
SZENE 168 


impuls in der Ladung (bzw. im Nahfeld) entspricht 2 kurzen starken Strom- 
stößen in zeitlichen Abstand 7 in entgegengesetzter Richtung. 


8. — Exponentiell abklingender Ladungsimpuls. 


Dieser Verlauf entspricht einem kurzen starken Stromstoß in einer Richtung 
und einer kleinen lang andauernden negativen Stromwelle, die allmählich 
verschwindet. Dieser Verlauf ist etwa aus einer Kurve abzuleiten, die No- 
RINDER (!!), angibt. Der spektrale Verlauf ist jetzt 


pe? 


DER ve vpl, 


und dies gibt den zeitlichen Verlauf 
: y 
exp [— ft] | exp [Pt]@(t) dt 


0 


ù 
y= 
y yt 


(1) H. NORINDER: Arkiv för Malematik, Astronomi och Fysik, 33a, No. 16 (1946) 
Fig. 25, 8. 35, und Consiglio Nazionale delle Ricerche: Propagation of Atmospherics, 
(1948), S. 168. 
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t 


yf = 98 exp Bel | exp [ft] G(t) dt + G(t) 
0 
G@) = E 1 : we 
D UT | 24/7 tl2 exp [—a/4t] . 


Die Funktion G(t) ist in Fig. 2 (S. 1131) dargestellt. Der Verlauf von y 
ist stark vom Dämpfungsfaktor 5 der Ladungsverteilung abhängig. Für sehr 
große ß nähert sich der Verlauf dem Wert d@/dt, wie er in Fig. 7, 8. 1135, 
dargestellt ist. Für sehr kleine 8 hingegen dem Wert G(t), Fig. 2, S. 1131. 
Man hat im zeitlichen Verlauf bei kleinen ß eine positive Welle, die durch 
Null geht bei t = «/2 und der dann eine schwache negative Welle folgt. Mit 
wachsender Dämpfung wird die positive Welle immer kürzer, die negative 
Welle wird größer und es tritt ein zweiter Nulldurchgang auf, dem dann 
schließlich eine lange, sehr geringe positive Welle folgt. 

Ein Impuls, der bei t = 0 allmählich beginnt und dann wieder abklingt, 
läßt sich immer durch die Differenz zweier Exponentialfunktionen mit ver- 
schiedenen Dämpfungen darstellen. Sein Spektrum wäre dann durch 


js 1 eee 
es A 
I» HA DH AL VIE 


wo Bo > Pi ist, gegeben und sein zeitlicher Verlauf durch die Differenz der 
zeitlichen Verläufe, die den beiden Spektralfunktionen entsprechen. Eine 
genauere Diskussion scheint aber erst notwendig, wenn die Experimente ge- 
nauere Unterlagen über Anstieg und Abklingen dieser Impulse geliefert haben. 
Eine lineare Änderung der Ladung (oder des Nahfeldes) von einem kon- 
stanten Wert zu einem anderen konstanten Wert ist nach der Gl. i — — dg/dé 
identisch mit einem rechteckigen Strompuls und bringt dadurch Zeichen rech 
Fig. 5 hervor. 

Wir haben im vorhergehenden einige ganz einfache zeitliche Verläufe von 
Strom und Ladung, nämlich den Dirac-Impuls, dem rechteckigen Impuls 
endlicher Breite und den exponentiell abklingenden Impuls diskutiert. Man 
sieht, wie verschieden in der Form die entstehenden ausgesendeten Signale 
sind und wie stark sie sich während ihrer Wanderung in ihrer Form ändern. 

In der Hauptsache sieht man, daß Signale, die ursprünglich nur positive 
Ordinaten haben (unidirectional) im Lauf ihrer Wanderung auch negative 
Anteile in den ausgesendeten Signalen erzeugen. Da die grundlegende spektrale 
Funktion der Ausbreitung V p exp [— x«/?-p"?] zunächst nur für lange Wellen: 
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d.h. geringe Frequenzen gilt, und die geringen Frequenzen den zeitlichen Ver- 
lauf nur für größere Zeiten wesentlich bestimmen, dürften die angegebenen 
zeitlichen Verläufe für größere Zeiten richtig sein, während der Verlauf zu 
Beginn des Signals, für kleine t, zweifelhaft ist. Deswegen dürften die ange- 
gebenen zeitlichen Formen für große Wanderungsstrecken D genauer stimmen 
als für kleine Entfernungen. 


RIASSUNTO (*) 


Utilizzando il metodo delle autofunzioni singolari, si ricava una formula per la 
propagazione delle onde elettriche lunghe intorno alla Terra, tenendo conto della iono- 
sfera. Con l’ausilio della formula ottenuta, si ricava, per mezzo di una integrazione di 
Fourier, la propagazione dei segnali elettromagnetici dovuti a una scarica a scintilla, 
ammettendo per il decorso temporale della scarica a scintilla alcune semplici dipen- 
denze temporali, in linea di principio differenti fra di loro. 


(*) Traduzione a cura della Redazione. 
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Application of second quantization methods 
to the classical statistical mechanics. 


M. SCHONBERG 


Centre de Physique Nucléaire de V Université Libre - Bruxelles 


(ricevuto il 6 Ottobre 1952) 


Summary. — Methods similar to those of second quantization are applied 
to the Liouville equation of the classical statistical mechanics. An equa- 
tion similar to the Schrodinger equation of a quantized field is given. 
It is shown that the interpretation rules of the quantal type of the field 
formalism lead to the rules of the classical statistical mechanics, the 
particles being treated as indistinguishable. The application of the Fock 
treatment of second quantization leads to the introduction of wave funct- 
ions in phase space, the probability density of the classical statistical 
mechanics being the square of the absolute value of the wave function 
in phase space. The choice of the sien in the commutation rules of the 
field operators leads to symmetrical or anti-symmetrical wave functions in 
phase space. To each kind of symmetry corresponds a different statistics, 
Bose or Fermi, as in quantum theory. The Boltzmann statistics cor- 
responds to phase space wave functions without symmetry conditions. 
The wave functionals for the states with zero particles of the fields in, 
phase space are analogous to those of the quantized fields of quantum 
theory. 


1. — Introduction. 


The method of second quantization was first introduced by DIRAC 
quantum theory of systems of non interacting Bose particles and later extended 
by JoRDAN and Kreın (?) to the case of interacting Bose particles, and by 


\ 
( 


1) P. A. M. Dirac: Proc. Roy. Soc. A 114, 243 (192 
2) P. Jorpan and O. KLEIN: Zeits. |. Phys., 45, 7 


in the 
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JORDAN (?) and JORDAN and WIGNER (!) to the case of systems of Fermi ‘| 
particles. The second quantization methods are essentially based on two 
properties of the wave equation of à system of similar particles: 

a) The wave equation is linear and contains symmetrically the variables 
of the similar particles. 

b) Only symmetrical (Bose particles) or anti-symmetrical solutions (Fermi 
particles) of the wave equation are required. 

The Schrödinger equation for n particles of mass m 


= ely mf Neh OL | 
(1) ih - = > i ST 3 en fas 


@(x— x’) denoting the potential of the forces arising from the interaction 
between two particles (*) is replaced by a Schrédinger equation for the wave 
functional y, 


RR 
(2) ih = = KHitas 
h? d \2 
(3) TH = ums on | 2 (x) (=. p(x) dx + 
1 if à ! Y ! = ! 
id || y*(x)p*(x') D(x — x')y(x')p(x) dx dx’, 


y(x) and w*(x) being operators defined by the commutation rules: 


(4) [y(x), p*(x')].= d(x — x’) [y(x), p(x’) |4= [p*(x), w*(x)]1= 0. 


The signs and + correspond to the Bose and Fermi statistics, respectively, 
and 
(5) [A, Bs = AB Sai BA . 


The equation (2) for the wave functional 7, is actually more general than the 
equation (1), because it leads to (1) when the operator N°, for the number 
of particles 


(6) NO = | p*(x)p(x) dx, 


has a value » in the state of motion described by y, 


(7) Ne tel tree 
(>) P. JORDAN: Zeits. f. Phys., 44, 473 (1927). 
(4) P. JorpAN and E. Wiener: Zeits. f. Phys., 47, 631 (1928). 
(*) ® depends only an |x—x'|. 
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But there are solutions of (2) which do not satisfy a condition of the type (7). 
Moreover, the nature of the statistics obeyed by the particles is already de- 
termined by equation (2) and the commutation rules (4), whereas in the case 
of (1) it must be imposed as an extra condition for the solution. Equation (2) 
describes the states of motion of a quantized field and the particles are the 
quanta of the field. 

The commutation rules (4) do not involve the Planck constant h and the 
procedure of second quantization can be applied to the Liouville equation of 
the classical statistical mechanics 


dfn 4 ; 1 à 
(8) / = (He Tale = | > = — 5 > Pix xv), In} = — atin + 
u (FG woe 
le opr dp, dx 


pı denoting the momentum of the particle of position vector x; and f(x, … 
Xn} Pry ++) Pn) being the probability density in the phase space of the » particle 
system. We shall prove in section 2 that an equation analogous to (2), and 
leading to (8), can be introduced: 


SOLS 5 
(10) 0 von Ry , 


9 


(11) R =i (yr, p)(P 


2m 


» W(x, p) dxdp + 
1 


i iif 7 sk ! ! t Ü ! \ ! la 
= 5 | p*(x, p)p*(x', p') (Dix — x), p(x’, p )y(x, p)): dx dp dx’ dp’. 
The operators y(x, p) and y*(x, p) are defined by commutation rules ana- 
logous to (4): à 


D 


(12) (w(x, pP), p*(x', p')]z = 0(x— x')ô(p— p') ; 


(13) (p(x, p), y(x', p’)|+ = [ p*(x, P); p*(x', p')]+ = 0. 


The choice of the sign in the commutation rules (12)-(13) determines the sta- 
tistics of the particles, which are here considered as indistinguishable. There 
is an operator N,, for the number of particles 


- 


(14) Nop = | p*(x, p)y(x, p) dx dp , 


having the eigenvalues 0,1, 2,... and equation (8) is a consequence of (10), 
when 


(15) May (n = positive integer) . 


75 - Il Nuovo Cimento. 
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The possibility of introducting the Bose and Fermi statistics (*) in à clas- 
sical treatment is rather surprising, because the probability density f, is always 
positive and must be taken as symmetrical in the variables of the n particles. 
As a matter of fact, we shall see that the theory based on the equation (10) 
leads to the introduction of a function ©, such that 


(16) (aise RS Pose PA NO ce Paes aii: 
The Liouville equation (8) for f, follows from that for ©, 


(17) >: = (Wale Onn) n <i ILinOn ’ 


since the square of the absolute value of a solution of the Liouville equation 
is also a solution of the same equation. ©, is symmetrical in the case of 
Bose particles and antisymmetrical in the case of Fermi particles. Thus we 
are led to a kind of wave function in classical statistical mechanics. 
Equation (17) may be considered as the classical wave equation, the hermitian 
operator L, playing the part of classical hamiltonian operator. 

In general, a condition of the type (15) is not satisfied by y(t) and the 
number of particles will not be determined, so that the solutions of (10) 
describe grand ensembles. We shall see that the average value N of Ways, inal 
the state of motion described by 7, and computed as in quantum mechanics, 
is constant (time independent) 


(18) N = (20), Nop 20) = | x*ONopxl0) de, 


the integration being performed with respect to the variables used to describe 
our special kind of « quantized » field in the phase space of a single particle. 
4 must, of course, be normalized: 


(19) | x*ydu =1. 

The rules of physical interpretation of the formalism based on y are ana- 
logous to those of quantum theory, but they lead to the results given by 
classical statistical mechanics whenever the number of particles is determined. 
This fact shows that a statistical formalism of the classical type can be reduced 
to a statistical formalism of the quantum type. It is possible to introduce 


(*) There are of course several essential modifications introduced by the quant- 
ization. | 
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an operator for any mechanical quantity. The operator for the energy is 


as ; D? 
(20) de = | p*(x, p) . p(x, p)dxdp + 


| w*(x, p)p*(x', p') D(x— x')y(x', p')y(x, p)dxdpdx' dp’, 


f 
bo} Hi 


and does not coincide with the operator X which appears in the equation of 
motion (10). 

The passage from the field type of formalism to the ordinary description 
in phase space can be easily done by using the following expansion of 7(t): 


(21) AURA UE + 
ES Typ al: 
PPT 


The 7,(x,.--; Pı...) are the functionals of a complete orthogonal system defined 


| Ojft; Kr + &n; Pry 3 Pn) Xn(Æ1s vs En Pas u.) Pn) dA ... dt dpi pre 


in section 5. Each term in the expansion (21) is an eigenfunctional of Non; 
the eigenvalue being the corresponding value of n.  ÿn(%,...: Prs +++) de- 
scribes the n-particle state in which the particles are at the points (x, pi) 
in the phase space © of a single particle. The expansion (21) is closely related 
to Fock’s (5) method of second quantization, but it is more conyenient to 
handle than the Fock sequences of functions. 
| Oult; ry es An} Pise Ph) "dx... dx„dp,... dPn ; 

gives the probability of finding only n particles and at the points (x, pi) of 
the phase space of a single particle 2. Equation (21) shows clearly how 7(t) 
describes a non stationary grand ensemble. 

As in quantum mechanics, we can introduce time dependent operators and 
a time independent wave functional by a contact transformation: 
(22) A(t) = exp [iRt]A exp [— iX]. 
The time dependent density operators are particularly interesting: 
(23) Onlb; X15 ce.) in; Pis.) Pr)= 

= y*lt; &, Pı)y*(t; Xe, Pa) + Y*lb; Æns Pa) PES Xny Pn) +» P(t; Xu, Pr) - 


We shall prove in section 8 that 
eine rest Page) (sy ln sie) oak 


== | = D(x; — Xn41)0n11(6; Lys 000) Kny Xni15 Pıs +++) Pry Pr-ı) diet dPnis 4 
À 1 1 


N n+ 


(5) V. Fock: Zeits. f. Phys., 75. 622 (1932). 
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These equations are a generalization of those used by Yvon (*), Born and 
GREEN (?) and Kirkwoop (%) in classical statistical mechanics. The latter 
equations can be obtained from (25) by taking the average values of both 
sides, as will be shown in section 8. 

One of the most interesting elements of the present field formalism is the 
introduction of a kind of vacuum analogous to those of quantum field theories. 
The vacuum is described by a wave functional y, defined by the condition 


(25) | p*(x, p)y(x, p)dxdp y — 0, 


oO 


« being any volume of the phase space of a single particle. By introducing 
suitable variables, 7, becomes the product of an infinite number of ground-state 
wave functions of quantum-like oscillators, in the Bose case. The vacuum 
described by y is, of course, defined in the phase-space of a single particle, 
not in ordinary space as those of quantum field theories. 

The wave formalism in phase space corresponding to the wave functions ©, 
and the wave equation (17) is, in a certain sense, more general than that based 
on the second quantization methods, since it allows us to treat also the 
Boltzmann statistics by considering wave functions ©, which are neither sym- 
metrical nor anti-symmetrical (distinguishable particles). The relation between 
the wave function ©, and the distribution function f, is not given by (16) in 
the case of the Boltzmann statistics, we have instead 


(26) ETES RÉ Dita ioe > MONET Sr Preys Pr) 2 
perm. ic 

because f, dx,...dx,dp,...dp, gives the probability of finding » particles in 
the elements of volume (dx,, dp,),..., (dx,, dp.) of 2, without specification 
of which of the similar distinguishable particles is at any of the elements of 
volume. It will be shown that is it also possible to introduce rules of physical 
interpretation of the quantum type in the phase space wave formalism, which 
lead to the usual rules of the classical statistical mechanics. 


2. — The equation of motion of the states x. 


We shall from now on denote a point of the phase space of a single particle 


by t. The corresponding elements of phase space volume will be denoted 


(6) J. Yvon: La théorie statistique des fluides et l'équation d'état (Paris, 1935). 
(7) M. Born and H. 8. Green: Proc. Roy. Soc., A 188, 10 (1947). 
(8) J. G. Kirkwoop: Journ. of Chem. Phys., 14, 180 (1946). 
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by dr. It is convenient to split X into two parts: 


(27) R = Rkin + R pot » 
(28) Run =i [PH vtr) dr = — À fyx(np: PO) dr 
A L 2m’ = m. dx : 
Ona 0 %(-! 1 1 1 
(29) Rpot = 5 || p*(t)p*(t')(O(x— x"), pr )y(T)): dr dt! = 


ie ams me ID dy(z' 
=, ere | pte) Ee eee ne es 


= y(t) : dr dr’. 
Ts p(T) 

The integrations are extended to the entire phase space of a single particle ©. 
(When no indication of the integration domain is given, the integration is 
taken over that phase space). It is easily seen that Akin and Ayo are her- 
mitian operators: 


(30) Rin = Risin Rpot = Rot Gee TRY 
It results from the commutation rules (12)-(13) that: 


(31). Tye), Raul = (PE wre), Ly, Fun = Elan, PO) 


1 /1 


| [y*(x), Roo) = à | (Bi x"), p(op*(r)}ytr") dr, 
(32) | | 


(p(t), Root] == i | y*(c)(Pe— x), p(r')p(T))2 AT" ; 


(33) | D Pl 2 
Hence 
s ù ; 
(34) [p*(t) y(t), Rin] = — Ÿp: A p*(t)p(t) } , 
(35) [w*(t) p(t), Rootl = 4 | (D(x— x'), v*(r)p*(r')y(r')p(r)): dr, 


and by integrating both sides of (54) and (35) we get: 
(36) lions Rxin| = [NES rot — [None R| = (Ur 


It follows from (36) and (10) that: 


Cl 
(37) v dt (NopX) = Nop Ky = RN opX) . 
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If y(t) is an eigenfunctional of Non 
(38) N opX&) = NY) (n = e— number), 
the same will happen with z(t) for any value of # 

(39) Nop x(t) en), 


because N,,7(t) is the solution of (10) determined by the initial value ny(t,) 
and the same holds for ny(t). The eigenvalue n can only be 0 or a positive 
integer, as in quantum field theory. (We shall prove later this theorem)- 
We shall denote by y,(t) any functional satisfying both (10) and (39). 

We can now obtain a very simple relation between the equation (10) and 
the Liouville equation (8) 


ù 


(10) NL Late (Ur ter Co) 

with 

(41) Yall Thy ens Ta) = lta) p(t) Wenzl) = LT Ya) - 
It follows from (12) and (13) that a 

(42) [y(t'), p*(z)(t)] = d(¢—t') y(t) , 

with 

(43) d(t— t') = d(x — x’)d(p— p'), 

hence 

(44) TL p(Ti), p*( = 2 à = Go) II y(t.) 


It is easily seen, by taking into account (12), (13) and (44), that 


(45) Il UT), PT )p*(T Le y(t’) = O(x — x) = 
= | TI vee, wee dpa 220) «22 
— 9) 3 LTT vtr), pro “EO = 
= php) 3 rn) (ox 2), TT vr) + 
+ oe ve ie 28 


— 0(t — Tt’) > (oe > TT a) O(t — T,) - 
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FU | 
We get from (2 
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28) and (44) 
2 n | n pP: n 
(46) m p(Ti), Rain = i( > > =, H ve) , 
Te ne U rea n 
and from (29) and (45) 
a n 4 n 
(47) he p(T), Ryo — à > ey D(x, a Xe) in) yr) =- 
t=1 IN Stier 7 XL n 
+ à | p*(x) | > d(x— x,), y(t) Il vie) dr, 
« r=1 ae 
hence 


(48) pi vn), R|- iM. TI y( (x) cie 


ei 


=F 0 [ y*(T) (2 P(x — Xr), w(t) i 


t=1 
It follows from (10), (41) and (48) that: 


d 
(49) 1 5 Ynlts Try + 


n 


TT vtr) Rott) = RTL verdzal0 + 


+: (mu, 1 IH p(T 1) et Alt) + a [ vita > P(x — x,), p(T) II u) Aryl) - 
1=1 2 r=1 L=1 n+1 


By integrating poth sides of (44) we get 


(50) TT PT) N op = (Nop ac n) in y(t) ’ 
WE 1=1 
hence 
(51) ny(t) II WT) Zall oy (TN op Xnlt) = 
Et 
= (Noy + +1)y D IL (T1) Anl) > 
(52) Nopp(t ily = — y(t) TL peor 
= l=1 


Since V,, has no negative eigenvalues, it follows from (52) that 


(53) 
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By taking into account the definition of A and (53), we see that the first and 
third terms in the right hand side of (49) are nil. Thus we get equation (40). 
It is still necessary to prove that there are non vanishing functionals of the 


type y,(t; T%,...), this point will become obvious later. 


3. — The unitary operator of the motion of the states y. 
The operator A being hermitian, exp [— U] is unitary. We shall use 
the following notation: 
(54) U(t) = exp [— ÜUR]. 
The Liouville equation (8) can be written as follows 


ore fr 
9% ( 
po) ot 


== Dale ; 


L, being the hermitian operator 


ig CS > of, -) 
Lx ; 


(56) lbp, = 0 > +. LES 


We have obviously 


(57) y(t) = a(t) 7(0), 
and 
(58) fu) = exp [— üL,1f,(0) = Un(t)fn(0) . 


We have already found a relation between equations (8) and (10) in section 2. 


Now we shall go further by showing that Q/(t) can be expressed in terms of 
the U,: 


[y(cy (Cpa) pr) Unlt) ya)... pr) } dr...dr, Pr. 
The P, are the projection operators corresponding to the eigenvalues n of N,,: 


co co 
(60) Nop — > nP, ’ > P,=1 ’ PE = On . 


n=1 n=0 


It follows from the spectral decomposition (60) that 


(61) Masala == WIE 
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hence for any 7: 
(62) Non Ln X) =F ah 


In order to prove (59) it is sufficient to show that the right hand side satisfies 
the two following equations 


oo 1 + À 
(63) JA SE » a | p*(tn) no UA) KG) ooo Gea) Ola coe Chia Ie, == QVM) = Al , 
n=1 a 
Ve y Ir, > as [ pr ).. vn) U nlt) { pity) ... pr.) }dty ... At, P | = 
dt | nl J n 1 n v 1) +. PAT) gATı «+. n nf 


| p* (Tn) 292 y*(t,) U,,(t) { p(t) ECS PT») } dr, O20 dr, 12 | ’ 


(65) i = RU ; U(0) — 1. 


Equation (63) is a consequence of the second equation (60) and the following 
one : 


(66) — | W* (Tr) oe P* (Tr) Ptr) pn) AT: --- Chien Hn > Rn e 
Equation (66) is obviously true for n=1. I it is valid for m = r it will 


also hold for n — r + 1, because Y(T)%Xr+ IS a %,, AS a consequence of (50), 
and we have: 


il s 3 
(67) m mt D*(Tr41) + Y(t) PT) «-- P(Tr41) AT... AT... Za = 
2 i * | En * * ; 
Se | Ole) | a | w* (Tr) oe p*(T1)Y(T) --- w(t) dt, ... At, | Y(Tr11) Yru = 


ite, gn 
— oi | AT:1Y*(Tr41) D(Tr41) Ari = à + NO rie 


It results from (62) that P,y is a %n, provided it is not nil, therefore we have 


(68) | — | ay (Tn) rene y*(t1) p(t) + (Tn) AT «+ ie, pur 


and (63) is equivalent to y IP a 


n=0 


M 
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It results from (48) that: 
(69) [R, p* (Tr) 300 p*(t)] — {Lp { p(T) QU W(Tr)} a ae 


— if(S D(x — x), p*(Tnr) … penn) w(t) dr. 


n+1 
Hence 
1 ie : - 
70) Rf tn) pm) OA (0 { ltr) … pre) fT … An yn = 
er Hat, ; 
| w* (Tn) ... W*(T) On(t) { RY(T1) .… p(t) „ar, we Ta Yn + 


= : = 
Be se | ya). P* (Tr) Ln Unlt){ PT) ++ p(tn) } ty ÂTa Xn + 


LA 


+ Sf yr Tr) pH) Y*( (nf 2 oe, Tal y(t) p(t). vn} AT. AT Ys 


n+1 


since the Z are hermitian. It results from (53) that the first and third terms 
in the right hand side of (70) are nil, therefore: 


al » 
(71) —R | y*(r)… p(t) Un(t) {ylts) … p(r)} du …… den Pa = 
2 N | y* (Tr) aie p*(n)U,(t) (p(t) re y(t) } dr, woe due 0 


Equation (64) follows immediately from (71), since XP, — 9. Thus we have 
proven (59). 

It is interesting to notice that the expansion (59) does not correspond to 
a perturbation method. If we would know the U, for all the values of n we 
could get exactly Q(t). Usually we know only U, and U,. In the trivial 
case of non interacting particles we can determine all the U,. 

In order to proceed further we must investigate the nature of the spectral 
decomposition of V,,. This will be done in the next section. 


4. — The spectral decomposition of N,,- 


We shall now prove that the eigenvalues of the operator N(w) 


(72) N(w) = | p*(t)y(t) dr , 

o 
© being any volume of the phase space 2 of a single particle, are 0 and the 
positive integers. In particular when wm coincides with the entire phase space, 
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the operator N(m) becomes N,,. Let us introduce a complete set of orthogonal 
and normalized functions for the volume w, which will be denoted by u,(T; ©): 


(73) [ui D ed aaa ome 


co 
(74) > us(t; our; w) = Ô(T—T'). 
We shall expand the operator w(t) in terms of the w,: 
co CO: 
(75) w(t) = > au,(t; ©) k(t) = > 0.us(T;0).» 
s=0 s=0 
The commutation rules (12) and (13) can be obtained by taking: 


(76) ds) Gis la. = One’ , [as Oy |= [sg Go| = 0. 


The expression of N(w) in terms of the (a, a*) is: 


(77) N(w) = > US — ¥ W,(o) 4 


s=0 s=0 


It is easily seen that the N, commute: 
(78) [N.(w), Ns(œ)| = 0. 


The determination of the eigenvalues of the N, is particularly simple in the 
Fermi case, because of the symmetry of the commutation rules with respect 
to the a and a*. aa, and aa; must have the same eigenvalues, therefore 
N, and 1— N, have the same eigenvalues, which must be 0 and 1. (The X, 
are projection operators). ; 

In the case of the Bose statistics, the determination of the eigenvalues is 
somewhat more complicated. The procedure is exactly the same as in quantum 
theory (harmonic oscillator). Let n be an eigenvalue of N, and ¢, the cor- 
responding eigenfunction. aif, is an eigenfunction corresponding to the eigen- 
value n + 1, because 


(79) d'a ar) ar (a a) ds Gn = I eer, 
and asf, does not vanish since n must be positive and we have: 
(80) (a2Cn5 an) =: (Cry al) = (% + ING. 5 So) e205 
On the other hand we have 


(81) ata @slr) = aslataln) — Abn = (W—A)aln , 
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so that af, will be an eigenfunction corresponding to the eigenvalue n— 1, 
if it does not vanish. Since 


(82) (On, Msn) a Was RE) = Men, on) ’ 


a.<, will not vanish if n 40. If N,(w) would have a non integer eigenvalue, 
it would also have negative eigenvalues, since it would be possible to get suc- 
cessively non vanishing eigenfunctions corresponding to eigenvalues n— 1, 
n—2,.... Therefore the eigenvalues must be the positive integers and 0, 
because from any positive integer eigenvalue follows the existence of all the _ 
others, since we can form successively non vanishing eigenfunctions correspond- 
ing to all the other positive integers and 0. 

In the case of the Bose statistics, the formalism is closely related to the 
quantum theory of the harmonic oscillators. This is easily seen by splitting 
each a into a hermitian and a anti-hermitian part: 


ct 1 r a 1 g 
(83) thy = VE (Ps — tas) , ie = VE (Bs + is), 
(84) 1% as" | = —— 10,7, os By) = [0 sr] == Nie 
We have: 
= 7 Me 2 7 1 92 2 
(85) N.@), (ra Le No) ED pro 1), 


Thus we have for N,(m) an expression similar to that of the hamiltonian of 
a quantal harmonic oscillator. In the case of m = Q, we have: 


1 
en Nog = 5D (dal). 


EB, 


It follows from the theory of the harmonic oscillator that the 7(x) describing 
a state in which there are no particles is 


| 


CO 
(87) % Ne) = exp — 


s=0 


— (0) 


DE 


s 


Ÿ 


1 2 
5% 


je 


This wave functional is not normalized. Strietly speaking, it cannot be nor- 
malized, because the normalization factor for each oscillator wave function 
is x 1, so that we would have to multiply 7,(a) by a vanishing factor. This 
difficulty is not serious, because in the physical interpretation we shall need 
the probability distribution corresponding to the square of the absolute value 
of the wave functional, exactly as in quantum mechanies. It is well known 
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that a probability distribution in an infinite dimensional product space is 
determined by the probability distributions in the factor spaces (see, for instance, 
P. R. Harumos, Measure Theory, (New York, 1950, chapter VII). In the case 
of the eigenfunctions of the N(w), the factor spaces are one dimensional and 
the corresponding probability densities are given by squares of normalized 
Hermite functions, as known from the quantum theory of the harmonic 
oscillator. 

In the case of the Fermi statistics, the present formalism is closely related 
to the quantum theory of spin. Indeed, the commutation rules (76) are 
satisfied by the Jordan-Wigner operators 


1 Te 1 
(88) a, = 9 ES (0° — io”) A We = € Ego), ae a) 
(89) Ets) == GO Ms of? À 


the o* being the Pauli operator spin vectors for a system of an infinite number 
of electrons. It results from (88) and (89) that: 


a 
(90) dd Aria, Je 


Thus we get a new proof that the eigenvalues of N, are 0 and 1. In order 
to get the wave function of the state in which there are no particles, we shall 
use a representation in which the asa, are diagonal and a wave fanoHonal 
depending on the eigenvalues of the Ne 


(91) seen = TT born, 


s=0 


In this representation the integration with respect to the N is an ordinary 
sum and we get 


(92) NR NN: 


so that the wave function of the vacuum is normalized. In the case of the 
Bose statisties we can also diagonalize the N, and a normalized wave function 
of the vacuum is given by a formula similar to (91): 


(93) Kol N)= 1] den? (Ko), xo (N)) =1 


s=0 


We shall now introduce the %s, 


sh 


n 
* ae * 
(94) VA ie Gs, As, ve. Us, Xo = IT as Lo - 


we = 
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It is easily seen that: ‘| 
n n 
* * 
(95) Nor 11a, = ee 
| ven b= 
Hence 
(96) N ops, ..8, —— Hl RES ; 


so that 7, will be an eigenfunctional of N., describing a n-particle state, 


if it is not nil. Since for any value of n 

(97) (54.08, 9 Las) = (Kor Ws, ee Us, Us, 0 Ms Yo) > 
we have in particular 

(98) (ls Ks) = (Kor 4s,4s, Xo) = (Lor Xo) #0, 

so that: 

(99) X, F0. 


Now we must distinguish the two statistics. In the case of the Fermi sta- 
tistics x,,..,, vanishes obviously when two of the s are equal. If they are all 
different, we have 


- 4 * = * * * 
(100) Us, ve Us, Ws, se Ms, = (a5 a5) (as, Us, _,) ++ (ds as), 
so that 

(101) (Xs, meas Key le) = (fo, Xo ys Su sm for LAT. 


In the case of the Bose statistics %,,...,, is always a non vanishing fun- 
ctional. Indeed, this is true for n = 1, but if it is true for n =r it will also 
be true for n =r + 1, since 


19 * * * 4 , * * } * | 
(102) As Us, r+ Us, Us, ve Ws, ss = Ws, ve As Us, vee as (Gs. Ws, 24) + 
* 
7g (Os. spt + 9s, Spt as Ei Ge) Mate Ue, lente) 
and 
io 

(8 (B) aa ny (B) (B) en 

(103) Ws Ka) = (1 Sr > Os, mt Ki) Gi 
i=ı 


(8) (8) En 
de, i (Ksı-.5,) Xs, .5,) TE 0 F 
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Moreover, we get by successive applications of (103) the following equation 


(104) Cae Vice) = > Ös se, Os, Sp, see 0) 


2 


k 
perm, k 5 % 


8) WEIN = (CB, „(B (8) 
55 (Xo ’ ho ) = n! Oe (Gas de 40 ) , 


with 


(105) Ne D ö 


N! perm. k 


We can write (101) as follows 


EF) yf) as GDF UE) UF 
(106) Cr Kane) = 1! Ons (Yo 140). 
with 
de Ô APRES, SE Ö, , 
| 
1 | On Oe ae: One 1 
rn NOTA 5 \ 
(107) =. De: 
n 5 N. perm. k 2 >“ 
| Ô OR ER: Ô 
Ö,, st Os, So 05, 5% 


The above notations allow us to put together (104) and (106): 


(108) (Kos? Lanta) = W! Os s(%or Xo) - 


It is easily seen that y, , and yy.’ are orthogonal if the values of the 
Secon Seen 7 
set s’ cannot be obtained from those of the set s by permutation. Thereby 
we can write: 


(109) CARTER (Ye ta). 


It is also clear that 


(110) CO sige = 0 (nn). 
The above considerations show that the set of all the functionals (n! 09?)-1/2 - 
lasse with 8, <s,<... < 8, in the Bose case and 1 <s<..<s, in that of 
Fermi, is orthonormal when % is normalized. We shall now prove that this 
set is also complete. The operators N, are a complete set of commutable 
operators in a space S,. Therefore the states that are common eigenstates 
of all the N, form necessarily a complete orthogonal set. We can get these 
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states by taking all the positive integers n and 0, and choosing in all the 
possible different ways n numbers s, (allowing for repetition in the case of 
Bose statistics) and giving to the corresponding N, the eigenvalue 1 and the 
eigenvalue 0 to all the other N,. But thus we get exactly the Hos 82.08 
(GSK Ko Sake 

We shall denote by P,.,, the projection operator corresponding to the 
state 7... The fact that the 7... form a complete orthogonal set in S, 
leads to the equation: 


co 

(111) Ee Den es 
n=1 DESERT ES 

We have obviously: 

(112) a > Pas 


Thus we obtain the complete spectral resolution of No»: 
(113) NO D TD RE RE 


It follows from the definition of a projection operator that 


(114) Pox =(%> Wor (Zor %) = 1; 


and 


~ in * ‘ 
(115) n | Oren SNS = (10 a. A 0 as Xo ; 


therefore we have: 


(116) BORD RER Ae as, Pos, ie We,» 


nn 
In the Fermi case, P, has a very simple expression in terms of the a,a* 
co 


CO 
(117) Poo [aa TAN: 
0 


s s=0 


5. — The basic representation and the Fock method. 


Instead of using a discrete set of functions in the space 2, we shall now 
take the continuous set w(t, T) 


(118) (ur, T)u(t, t) dt = Öö(T—-T), 


(119) [ w(x, Dur, ddr = dr — 7). 
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~ 


The expansions (75) are now replaced by the following 
(120) y(t) = [au Dar,  p(n = [ar dur, Daz, 


and instead of the commutations rules (76) we have: 


(121) faz), (zz = 0(7—7), Ta(r);a(r')]i= [0*(2), a*(7')]1= 0. 


In particular we may use the orthonormal system of the ö(t— 7): 


(122) u(t, 7) = Ô(T— 7), a(t) = y(r), 


so that we can work with the y and y* themselves. In this case N, i 
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Ss Te- 


placed by N(w) (see (72)). The %,,.,, are now replaced by WOKS. (zent) 


(123) Pt ro) = TUT) =e (Ean) Yor: 


It results from (44) that 


(124) pop), LL Pe) = À Oe) TT vr, 
val En bes 
hence 
(125) No m vr NL) II y(t) ’ 
En 1=1 


Nos. cn) beinesche number of points 7, ... 7, inside wm. (For the sake 
of simplicity we shall not consider the case in which some of those points lie 


on the boundary of w). It follows from (125) that 


(126) N() (Tay +. Ta) = N(@; Ty +, Ta) X (Try ++ u)» 
and 
(127) N opX(Tı, ..., Tn) = NXT Tu). 


It is easily seen that 


(128) (Kran TL 29 To) = (Loy YCEn) + PR) PRCT) ++ P*(Tn) Xo) » 


and by using this equation we can show, as in the preceding section, that 


(129) (x(t; Un): AT: ..) Th) = W! OMT, 2 Yo» Xo) 9 


76 - Il Nuovo Cimento. 
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with 
‘ Sn) 1 1 \ ih . 4 
(130) Og (TT) = ET DRAC RUN) Ur) 
+ perm. k 
O(t%— 7) 6(t1 —T2) O(T, — Tp) 
| 6(t%2—7;) Im — 1) Ô(T: — 7%) | 

1631 ae = 4 
( ) F (T, T ) = n! 

dr) SO ae eee OT ir) 


There is also a relation corresponding to (110): 


(132) (are, at) (n Zn). 


The same considerations applied to the discrete case show that the (1, ..., Tn) 
form a complete set of orthogonal functionals. The situation with respect 
to the normalization is somewhat different, because of the continuous set of 
indices. We shall assume that 7, is normalized as before and we shall use 
a definition of the normalization which makes of the 1/V n\( Er «239 En) EOE 
malized functionals. 


We shall introduce projection operators P(t, ..., Tn): 


(133) BNE Gy, 2005 Ca) N= er ces) Eee 


= (Zor P(Tn) ++ P(A) 4) P* (Tr) «+ P* (Tn) Zo - 
It follows from (133) that: 


if ; 
(134) Pre ta) aie (a) nor (Bel (aa) er tes) coo i (CR) 
It is easily seen that: 
(135) Pee Ee a) == On (n # 


(136) PG met ae) Nee 


The analogue of (111) is 


(137) Pes je done, 
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and we have: 
(138) Pee | Pina. du. den 


The equation (138) is an immediate consequence of (137). In order to prove 
(137) it is sufficient to show that the operator in the left hand side applied 


to the y(T1, Ta, ...,7,) leaves them unchanged. This results from the follow- 
ing relation: 


(139) WEP ea, tree: older de gen, 
= | Aie TOC, Cdt dr, (nm) me 


The expansion of a functional 7 in terms of the z(t, ..., 7.) is: 


CO 


il 
(140)  % = do + À - = OCR en Ta ts. du. 
n!. 


n=1 


This expansion is equivalent to the Fock (5) type of formalism, in which any 
state of the field is described by a sequence of functions 0,(t,, ..-, Tn) 
(n = 0,1, 2,..., co). The representation whose basic vectors are the 7(T, ..., Ti) 
will be called the basic representation. The Fock functions ©, are the com- 
ponents of the vector 7 in the basic representation: 


1 Il ne ; 
(141) Dre) Vai (Tien) = i (x. 1a w*(T1) Lo | - 


We get from (140) and (129) the fundamental formulas 


(142) tee | 10,’+$ | | On (ty, vey Tr) |? dr dr}, 
and 
(143) BR) 190, ae > | OE Cig sce OAC cee A due dr, | : 


The matrix elements of an operator A acting on the x, in the basic repre- 
sentation, are 


1 


nm ee en 


VE ER 


D us 
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we shall now introduce the n-particle part A” of an operator À: 
(145) Am— Ab, . 


The only non vanishing matrix elements of A” are: 

+ =) 

a ‚| fa , ms | , IN 
(146) Crest AT TE a ea ERTL 
In the case of an operator which commutes with N.) we have 
(147) [4,P.]=0, 
so that: 

CO 
(148) AND. A > Am, 
n=0 


The density operators o, defined by (23) commute with No» 
(149) ENS On(T1 sees Tn) | = 0 ’ 


as a consequence of (44) and (125). Since 


co il ir | 
BON a> aa, | péter) PT my) AU ag oy AE 
(Gi) ica > | Pire mdr Aer lies Cas 


| Pin. Ca) duridriv(sr) vale 


a y*(t,) ere p*(Tn)Y(Tn) tee p(t) = On(T1 Aa Ta) : 


we have: 


To any additive (*) physical quantity F(T,, ..., 7„), depending symmetri- 
cally on the variables of particles, corresponds an operator Fo, 


al . 
(152) hy = = | TGR AE Onkt gs Ha tee At, & 


*) The value of the corresponding quantity in the case of a system of 7 particles 
I sd \ yi i 


(r > n) is the sum of its values for all the | | different groups of particles. 
n 
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We assume that the quantity # has no meaning for a system with less than 
n-particles. The first term in (20) is the operator corresponding to the kinetic 
energy of a particle and the second term to the potential energy of a pair of 
particles. It follows from (151) that 


(153) I = D (") | F(x, TL (TARE Te) TT ae dr, 
thus: 
(154) Fo, = (,) [Fa Pac tdenedr. (r>n). 


We get from (154), (140) and (129) 


(156) (4 P20) = (] [Bas ml Only ey te) [PAE oe dr. 


As particular cases of (156) we have: 


: (Pi Pr ja 
(157) (X in) | ie ee | lt, Tr) | Ari... der, 
f de .r 35 
(158) CALAIS | > D Plxi—«y)| Old. t)| dt... de. 
mA Y= Üer 


In the derivation of (157) and (158) we took into account the fact that 
l'O, (4, +. Tr) |? is a symmetrical function of the 7, ..., 7, in both statisties. 

We shall now prove that, in the case of a solution z(t) of (10), the O, satisfy 
the Liouville equation (17). It is easily seen that 


(159) Ry (Try 000) Tn) = — UM (Try + Ta) AT se Tn))n » 
by taking into account the relation (48), hence: 


ee) 


(60) Ry(t) =— Y= | Alt; tr Tu)" 
n=1 Vn!: 
SH, (T1 2. En) X (Tay <-> Tn)}y AT ... IT = 
oy > RSS | (A(t; SOSE) Tn)s O, (t; Ty +.) Tn) a X(T ooh) Tn) AT oe dr, = 
n=1 Vn!- 
.dO et a”) 
= =: %o + À ds | er Dre rit )dtie. Are 
n=1 Thane 
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Therefore we have: 


The expansion (140) of the solutions of equation (10), and the equations 
(161) can be derived from (59). Indeed, we have: 


(162) x(t) = Ut) y(0) = 
Le 
= a0) a > | p*(r,) … p*(ta) OL PH) pen) } du ++ dr Pax(0). 


It is obvious that 
(163) Pyy(0) = 9% - wat ©, = constant, 


and as a consequence of (138) and (133) 


J if / / 1, ‘à VA / 
(164) P,Y(0) = ii 110,0: 7, mr eta) OG ca Tae, 
n!- 
with 
~ fé al 74 [2 
(165) 9,0; Ty, - Ta) = Vi (x(T, .., Tn) 4(0)) : 
Vn! 
Since 
(166) (Gals + Ta)ı Kalt ares a) 


= (Yo Y(Tn) +. Y(t) Y* (72) -.- pt) Xo) = NL OM, 7’) , 


we have: 
(167) WT) --- W(Tr) (Ty) … W*(T,) Xo = n! dr, T) Ho « 


It results from (162) and (167) that: 
(168) y(t)= Ox + > B | 0) (nv (Gt OR O38 Tay ST). de 


co il TRE 
= 05% + - = | II. ran TaN dr, "dre 
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This is exactly the expansion of y(t) in terms of the y(t%1,...,7,), with 
O,(t; Ti, ..., Tn) Satisfying (161). We get from (59), (138) and (134) by taking 
into account (167): 


(169) UD —P,+Y [wre ... pr) Po Talt) (p(n) ee pl} An... dr, . 


It results from (168) and (167) that: 


1 


(170) U,,(t)0,,(93 Ti; DOF) Tn) Xo — yp 
vn! 


CE) Er Or 

This formula is a generalization of (40). Indeed, in the case of a n-particle 
state 7, 

(171) a) = An(t) ’ 


and we have: 
(172) p(T) = (Tn) xnlt) = Vn! Un(t)On(05 Tu, …, Ta) Xo 


The quantity in the right hand side being obviously a solution of the Liouville 
equation (17), that in the left hand side must also be a solution of the same 


equation. We get from (172) and (41): 


n(n—1) 


(173) pus ten tr) = Vale 0; mn... tr) Xo, (al, Hd), 


6. — The rules of physical interpretation. 


Let x be normalized: 


(174) (4 0) =|)? +> [ | O, (Ta; +.) dr... tn = 1 : 


n=1- 


We shall assume that x deseribes a state of a grand ensemble and that the 
probabilities are given by the following rule: 

(I) The probability of finding n particles (and only n), in the elements 

of volume dr, of the phase space of a single particle, is |O,(t1, ---, Tr) | dr dle 

This rule is altogether similar to the corresponding one in quantum theory. 

We shall now see that it leads to the usual interpretation rule of the classical 
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statistical mechanics. Let us consider a normalized solution y(t) of equat- 
ion (10). It follows from (161) that the function 


(175) TAGS Tao tale OM: finan tae 


is a positive solution of the Liouville equation (8). Therefore the interpretation 
rule (I) will give the usual statistical interpretation if we take for y(t) the 
solution of (10) determined by the initial values 


(176) 9,0; Ts Tq) = VPs Ty «+5 Fa) 5 


the p,(t1,.….,7,) being the probability distribution functions of the grand 
ensemble in consideration at the time t= 0. The special case of an ensemble 
corresponds to a wave functional of the type y,(t). The choice (176) is not 
the most general, since we may choose arbitrarily the phases of the ©, at the 
time t = 0 without changing the probability distribution. 

All the mechanical quantities which are considered in the classical sta- 
tistical mechanics correspond to operators F,, of the type (152). It is easily 
seen that the average value of the quantity F(t,,...,T,) 18 


(177) F(t,, ..; Tn) aa (X I) ’ 


as a consequence ot (156). In the present formalism we may introduce more 
general operators than those of type (152). It seems natural to assume that 
those more general operators may also correspond to physical quantities, when 
they are hermitian. One such operator is X. It is clear that 


(178) y(t) = (x(t), Ry) = (ru a a x): 


is a kind of measure of the rate of variation of the state of the non stationary 
grand ensemble. The consideration of these more general operators requires 
a more elaborate notion of system, or even grand ensemble, than that of the 
classical statistical mechanics, since it is not enough to give the probabilities 
10,7? but the ©, themselves are necessary. From now on we shall assume 
that the states of the grand ensembles in consideration are described by wave 
functionals y. We can now extend the interpretation rule I: 

(II) The physical quantities are described by hermitian operators A. 
The possible values of A are its eigenvalues A’. The probability of a value A’ 
in the state described by the normalized wave functional y is (y, Py x), Pu 
being the projection operator of A corresponding to the eigenvalue A’. 

(III) The probability of finding the values A7 of a set of commutable ope- 
rators A, is (y, Py, x). P,, being the projection operator on the linear mani- 


APPLICATION OF SECOND QUANTIZATION METHODS ETC. - I 1165 


fold of the eigenfunctionals of the set A, corresponding to the eigenvalues 
A, and x the normalized wave functional of the state in consideration. 

Let us apply the interpretation rule III to the operators N(w,), the w, 
being non overlapping finite volumes of 2 containing the corresponding points T,. 
We shall have to consider also the operator P,, if we want to have the pro- 
bability of finding just particles. Denoting that probability by p,(@,, ..., @1), 
we have: 


(179) PulOdsy <2) O1) = (x, lee dre. 1.) = 
if ; 2 à ; : 3 ‘ 
oe | dr | canner | dr re) | dita | die | dt OT: CA te 


on Mo Opn, Oy (We On 


Thus we get the rule I as a particular case of the general rule III. 
The average value of the operator A in the state 7 is: 


(180) = > ie) = 
A! 

since 

(181) SAP — A7 


The quantity v(t) defined by (178) will be called the state-velocity of the grand 
ensemble. It is the average value of A at time Et: 


Some of the general operators we are considering have also an equivalent 
in the theory of the Liouville equation. Thus X corresponds to the Liouville 
operator L, of equation (56). We get from (160): 


corral 


a) Ry) = ZF | aly 


MU) Lig Obs Tine t,) Als Un - 


Hence 


(184) (x(t), Ry) = > | OS Ta) LOWE; Tay en La) dti des 


= 17 


in particular, in the case of an ensemble with a fixed number of particles n, 
(184) becomes: 


(185) v(t) = JO: Tay ony Ta) Ln Owls ti. Tn) ty +. AE « 
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We may consider ©, as a kind of normalized wave function in the phase space 
of the n particle system, since in this case 


(186) os DNS TS) Pan, u 


and the quantity in the right hand side of (185) as a kind of average value 
OEL 


7. — The above considerations lead us to introduce a set of interpretation 
rules in the phase space of a n particle system, which generalizes the ordinary 
rules of the classical statistical mechanics: 

A) The states of a n particle system are described by wave functions ©, 
which satisfy the Liouville equation (17). 

B) The physical quantities are represented by linear hermitian operators 
æ« acting on the wave functions ©. The possible values of the physical quan- 
tity represented by the operator « are its eigenvalues «’. 

C) The probability of finding the values &; of a set of commutable ope- 
rators œ is | Otis, Tr)%x!On(T1y +++) Tr) AT +. dt», Xr being the projection 
operator on the linear manifold of the eigenvectors of the x, corresponding to 
the eigenvalues «;, in the Hilbert space of the ©,. 

These interpretation rules are equivalent to the ordinary ones of the clas- 
sical statistical mechanics when the only operators in consideration are of 
the form F(z,,...,7,), but no more so when other operators, such as L,, are 
also taken into account. It is easily seen that the average value of the quan- 
tity « in the state ©, is <a) 


(187) KO — | De) O0 (Er ere 
with the normalization (186). In particular we have: 
(188) i (te 


We shall now generalize the method, given in section 5, of relating ope- 
rators in the phase space with operators in the space 2. We shall take as 
the generalization of (152) for an operator «” depending symmetrically and 
essentially on 7 particles, i.e. having no meaning when applied to a © of a 
group of less than r particles, and corresponding to an additive physical 
quantity, the following relation: 


1 
(189) Ayo 


NO) 


| DT)... Pr HT) pr) + AT: +. AT. 


~ 
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By using the integration kernel Æ,4 

r / / à 2 Rn #15 
(190) K (E17 es Cri Try ses Tr) = AN {dt — te OT,  Tu)s 


/ 


and the generalized density operator 0,(T1, ..., Tr; DT) 
(191) lee OMAN CNE CACAO ae 
we may write: 


out 
(192) Ayn =, | WO ia Vien TNO Ciena Les is sor te) [I UT Ur. 
dive: 


We shall also need to consider generalized projection operators P(t), ..., T} 
LA / 
eek): 


(193) Birne A ee wala PA Ra 
ie / al kk FA > ph 7 
(194) Pit... Te} Try ces Tr) = ri p(T) o p*(T)Pop(Tr) «+ Y(t) - 


We have: 


ee n! fi , y 4 
ae > em PA ALES PE ean Teta) 
De r)di dar dt, — 
co 1 
= p* (1) p*(t,) > (oe r!) | p* (Tr41) oe Ÿ (Tr) Pop(Tn) W(Tr41) At 
if / ee) x [24 
LE) vena Nai" ae) > | Dia... „u Ot na ie 
m —= 0 
y(T,) DOS w(T;) re Ÿ (Ti) GI Ÿ (T;)p(T;) (T1) 
Hence 
(196) Or N 
= n! i ; , / A Tees 2 
=») ee DU ee Fale OR u) dan 
re) oats 
so that the generalization of (153) is: 
NT ; , , , 
(197) Aya = > (") | Tas Tia Or) Clio: Dita" TE 


öl der... den Osh heat dee 
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It follows from (197) that: 


CO fn 2 ” 
N n = / / = 7 7 m m 
(198) AymAgs = > ( l | | rue co TON Enno lise Lit oo te) © 


SIR 
nzr,s \ 


m 4 


HAT, T')P(Tiss, Tr} Ty 9220 Tu ltr — Tr41) ve Ô(Tn — Tp) * 


(2 


Asa — Tey) sor Ol Tp = dem, dr. dr OT dd dede tre 


It is easily seen that equation (198) can be written as follows: 


co 


(199) À, À Bts) = » | Ko 0 (Ta ; ... Uns Tag .... T,) ‘ 
n =D" 
Pro tata du. dr dre den 
il = 
(199a) p=z(r+s+ir—s|), @®= ) "ln, ..Ter)d. 


5) 
2 AIR ae teas 


Equation (199) shows that: 


(200) [A su, À pis), =) when [oder p®)] ==. 

It is obvious that the interpretation rule A for symmetrical or anti-sym- 
metrical @’s is a consequence of those of the y formalism. The same will 
happen with the rules B and C, if we associate to each « the r-particle part 
of the corresponding A, and to ©, the wave functional 7,[0,]: 


i! 


(201) 4.9.) = —= 
%1O,] Be 


| Dee tir sr) dre 


By introducing the spectral decomposition of «, « = > 4/7, into (189) we get 
7} 


x 


A eh) —— > Ca ; 
a! 
and the A, , are projection operators, as a consequence of (199): 


(202) AE dune — Oxi An, ñ (Tri Tan = drain!) . 


x! 


The eigenvalues of «” and AY are therefore the same and the corresponding 
projection operators are related according to (189). The probability of finding 
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the value « of A® in the state y,[9,] is: 
(202) (4101, A, 
EN 
= a | O; (Try +++ Tr) (%r(T OU) Tr); Ark (Ca ooeg T,)) z 


OT WUT, le dr. dt, . 


It results from (197) that: 


(204) (Hit en) Au, ti, N Cire cry tan coule 
Hence 
(205) CAICAR 12 7101) = 


if r VA 2 LA 4 te A 
= | OF (Tay 2 Tr) Ka (Ta ces Tri Tay es OT; 2 bp) Os 22200, dt) dt 


so that the probability of finding the value x of A“) in the state 7,[0,] is the 
same as that of finding the value « of « in the state O,. In a similar way 
we can see that the probability of finding the values x, of the commutable 
operators A‘? in the state 4,[9,) coincides with the probability of finding the 
same set of values for the «, in the state ©,. Thus we have proven that the 
interpretation rules A, B and C result from those of the y formalism, in the 
Bose and Fermi cases. 

It is interesting to remark that the operators Rin, and Ayo may be obtained 


(2) 


from the operators L}, and L respectively 
kin pot? 2 


z ® 
206 Le, =i[P- | =- Re. 
a = 2m’ jy m dx ne 
WD à wD D 
207 L®, = 1D(x— x =f - + ; 
( ) pot \ ( ); is x dx I Dx! Dx! ’ 


by the application of the formula (189). 

The wave formalism in phase space can also be applied to the case of the 
Boltzmann statistics by not imposing symmetry conditions on the wave 
functions ©,. 
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8. — The Heisenberg representation. 


In the Heisenberg representation, the wave functional is time independent 
and the operators for the physical quantities vary with time, as in quantum 
mechanics. The passage from the Schrédinger representation to that of Heisen- 
berg is given by equation (23). The equation of motion follows from (23): 


dA 
(208) i = (4, 2] 


The operators commutable with X do not vary with time, they are the integrals 
of motion. X is obviously an integral and No, too. Other general integrals 
are the energy operator X defined by (20) and the operators of the linear and 
angular momentum 


(209) 9 = [v*(nytopar, 
(210) . “= PACE / par. 


Another important integral of the motion is the symmetry operator Yop: 


CAL) Y= Cove + > Wel | On(— Try — Tay cos — Tn) X(T1y es Tn)AT ee ATs - 
n=1VMN-- 


The equation of motion (208) is not valid for the physical quantities that 
vary with time in the Schrodinger representation. For such quantities there 
is a time derivative (dA/dt), in the Schrôdinger representation, and (208) is 
replaced by 

dA . oA vA 


212 et À 
PONS ee nt 


A pe 
= exp [iÆt] A exp —iXt]. 


The general equation for the integrals of motions is: 


(213) Fe Re 0 


In (213) the operators may be taken either in the Schrédinger or the Heisen- 
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berg representation. The component of & 


(213) ee = | y*r)pir)adt, 


are not integrals of the motion, but with & and @ we can form the vector 7 
(214) a a 
m 


whose components are integrals of the motion. 
In the phase space wave theory of section 7, the integrals of motion are 
defined by an equation analogous to (213): 


It is easily seen that in the case in which x, = RAT, .., 0a) (235) ZOes. Over 
into the Liouville equation: 


F 
(216) a Re CHE 0» 


Thus we see that the definition (215) generalizes the ordinary integrals of 
motion of the classical mechanics. Let «” be a r-particle operator which 
corresponds to an additive physical quantity, in the generalized sense of 
section 7, and is an integral of motion for systems with n > vr particles 


: dam 
(217) 0 


+ [a™, L,]= 0, 


re 


x" being defined by (200). It is easily seen that: 


ù 


(218) À anit) = N 


À . 
ot : 
we get from (217), by taking into account (189) and (48) 


3) 
(219) : v ot À ,6 + As, wg — 0: 


This general theorem shows that X and the components of G, ZT and J are 
integrals of motions. 
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No, is an integral of the motion, as shown by (36). We shall now prove 
that all the operators 77, : 


(220) Tl | Oneness, dns 
are integrals of the motion. Indeed, it follows from (48) that: 
(221) [Qr(Try ++) Tr), A] = Ip). pr), A] (tr) + pu) + 

+ p(t). PRT) (Tr) + Y(t), À] = 


Sr \(H,, 0,(T; DOOR) Teale = i | | = Dia an), 9%-+1(Tı, Li) Trin) Atria . 
E n=1 


r+1 

We get from (221), by integration: 

(222) (Gli Tel) = We 

Therefore the 97, are integrals of motion. It follows from (221) that the 
equation of motion of the density operators in the Heisenberg representat- 
ion is: 


ù 


(223) i cas me CEO tag rer 


Sr | > D(x,1— Xn), Ortill; Try +. Ta) dt, - 
N faa! 


We get from (151): 


(224) OES tay voy Tr) = (oi CO ose tr AG) eae 


n! 
=> ay | On? Tap Cree Ta) Op rer. 


In the particular case in which N,,y7 = ny, we have 


x (n) n | ‘ 2 
CP Fol), = nt | ON Te) ar Are (n >r), 
so that: 
Ben ern (n) il. à m) 
(226) Oral(t; Ty» ses 1) — | 0,(t; Ti. m) dr, 
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The al; 71, ..., 1) are the r-particle distribution functions used by Yvon (°), 
Born and GREEN (7) and KIRKwOOD (8). By taking the average values of 
both sides of (223) we get a generalization of the equations of motion used 
by those authors: 


Ned Az re 
or, nor, 


(227) 


n=1 


+ | | > D(x, 44 >= Xn), 9r41(Tı5 2007) rs) dT,41 © 
u r+1 
It results from (151) that 


(228) on) let, 


because the projection operators are commutable, as shown by (135) and (136). 
Therefore do,/d¢ is commutable with all the operators involved in the right 
hand side of (223), so that this equation involves only commutable operators 
and may be treated as a ¢-number equation. Thus we are led again to (227) 
but without averaging. It is possible to replace (226) by an operator equation 
of general validity, since 
(229) | Oil high coop ti) Chis — Wi) cco. We ea ne) = 

= (Wy Fe Dares). 


By successive applications of (229) we get 


Neos 


(230) TG = (Nop — ti DNS tr 2). Nop— DN = oes ryt? 


with 

(231) (Nop — 8)! =L(Nop— 8 +1). 

It is possible to express o, in terms of the o’s of lower orders: 
r-1 

(232) oo the NT let — > dm — r) = 


tr LA Tralee) —T dr a) (et — dr, fe x) = 


= e(t)(o(T;) Zen T,))(o(T;) öltı Ts) O(T. Ta). (ar) SE 


77 - Il Nuovo Cimento. 
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Equation (232) is also valid in the Heisenberg representation 
(233) OOP Gy ce te) =D AG Vis Ce) er Gt) Sa ölt, — Tr )) = 


= olt; met: To) — (m — 2) ee an Tr) 2 Ô(T, — r) 5 


so that: 


rod 
(234) OWS Tip onan Ep) = Opal? iy com (et: Tr) — > dm — 2) = 


= oft m)(elts 7) — in m)let: Ft rm): 
It follows from (234) that: 
(235) al; Try Te) = Olt; )elt; Te) — el; u) O(t1 — 7) , 
(O36) wein 

= 04(t; Ti Ta)O(Te) — Et; TH, T ») {OT SST ENG er) - 


In order to be able to determine the one particle probability density with (227), 


we must make some assumptions about 0% (t; T1 Te), 1.e. about o(t; z)o(t; To) 
In a similar way an assumption on p,(f; 7, 1,)Jo(t; A) is necessary to get the 
two particle probability density from (227). The ordinary assumption of 


kinetic gas theory is the statistical independence which leads to: 


(237) o(t; n)o(t; 7.) = off; t1) Olt; 7) - 


We shall now see that the kinetic gas theory hypothesis of the molecular chaos 
corresponds to a Hartree kind of approximation in the present theory, more 
exactly a Hartree-Fock approximation in which 


1 


(238) ren) = 
3 Ti... ee 


ul; Tr,) ++ Unlt; Tu) 


(239) OMAN 


Un Esti) MACE TS) ia. alt Ta) 
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the « being normalized orthogonal functions: 
(240) ea Du(t; T) AT = 6. 
We get from (225), in the Hartree-Fock approximation: 


(241) ea) = Sluts)? 


I 


> | wilt; m)urlt; 72) + eurlt; Tr)Ur(t; Te) |? 


(243) cr Le Er =— LE. 


The cross terms in (242) correspond to exchange effects. In order to get a 
clearer understanding of the situation let us consider the case of the Boltz- 
mann statistics, in which the Hartree-Fock approximation is replaced by 
Hartree one: 


(244) GG. Tay a) — M Vy Malls Te) > Wall} Tp) - 


In the Dan case equation (227) is still true, provided we replace 
DEREN, on u) by the r-particle probability density f72(t5 Ti, Try +, Tr) 


il 
(245) AVR REA M Ge [ tn Tn) dran Ara 
with: 
(246) lt; tT, +) Tr) = nl; Tu Te) |? = = > ul; te,)|*--|Un(E5 T,) | 
perm. Æ perm. k 


Thus we get 


(47) fr) =D |ults 2, fs tay Te) = XD [malls tr) |* ud; m)" 
1=1 1=1 Vl 


so that 


(248) FC Tay Te) = AP Ca) fa” (re wer); 72) |, 


a 


and for large values of n: 


(249) (O (ts tay Te) BAG Hh GB) 
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The above result shows that, in the case of the Boltzmann statistics, the 
hypothesis of molecular chaos follows from the Hartree approximation. The 
Hartree-Fock approximation allows us to extend the molecular chaos hypo- 
thesis to the case of indistinguishable particles. It is easily seen that 


(250) Se mel ire) (ta) 


— 5 À ut m)hlutés va) |? dautés m) eis]. 


Teer et — ail 
For large values of n we may neglect the last group of terms in (250): 
3 (b; Tas Tay Ts) fa (Cj Ty Tai (C3 Us) aa (ts 5 Ta) A(t; To) fr (4; Ta). 


This expression of fj? is unsatisfactory for the determination of f; by means 
of (227). 


9. — Variational principles. 


The equation of motion (10) can be derived from a variational principle: 
ù 
(251) 5 | (x0), 2} x(t) at = 0. 


There is also a variational principle for the wave equation in phase space: 


t 


ms 
(252) Ont ROE (Garces AU Zeke 157 | en Uo eos Ole, = Ory 
Lo J 


Lo 
It is easily seen that, in the case of the Bose and Fermi statistics, (252) is a 
consequence of (251). 
The variational principles may be used to get approximations, as in quantum 
mechanics. Thus, by taking a ©, of the form (244) and imposing only the 
normalization of the w, 


(253) | ul; Thuilt; tT) dr =1, 
we get from (252) the equations 


sik dUr P° F ; 
(254) Sn Om’ ult; T )) 4 N | (@ (x — x), url; t')ur(t; T)h,uw(t; t') dt’ + 


+ A,(t)u,(t3 7) — u(t; 7) > i drult; t’) | (ee : ur) == 


VAI 


> [uns tz; r')(D(x' — x"), Ur) 9 At" Ke 


a Te a 
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the A,(t) being simply related to the Lagrange multipliers corresponding to the 
conditions (253). Since 


2 
Re 
2m’ 


(255) à nuts n= ( 


Bs | (D(x — x’), | w(t; 7) |? | wv (t5 t')|?), dt! + 2A,(t) | w(t; 7) is 


EE 


we have 


We get from (255), (256) and the first equation (247) 


ù 2 ; 
(257) > jt 0 =(F APG 7) + [ (@~e— x), (PE; HIPC 7) ede’ — 


2m 1 


— 5 [(8&— =), | wilt; nt; 7’) dr’. 


The last group of terms in the right hand side cancels the self-force. Neglect- 
ing this correetion we get an equation involving only fi 


red ù n) ä (n) il Ty g(r ! ! (n) 
ess) > ne n= (men) +( [Bere er men), 


which describes the motion of a particle in an averaged field of force. 
We shall now obtain some conservation principles from (252). The. phase 
space integral in (252) is invariant for space translations. A translation dx 


(infinitesimal) of all the particles transforms ©, into O, + > (00,/dx;) dx; 


LES 

since the integral in (252) remains unchanged, we get 

2 » N N) . n à] 

BI [Ost > = — OA) At dé = | (to) I — Outo) dt... dt, = const. , 
: 1=1 0% . 110% 


by taking into account the Liouville equation. The average value of the 
nr 
operator > d/dx,: is therefore time independent. It is easily seen that > dfox, 


HA 
is an integral: 


(260) [2 m A za), 
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In à similar way we get 


GE a ® ù ù D 
= CRG go ge OEL con Oe, — 10 
PAL ES mi LU 


by giving an infinitesimal rotation to all the position and momentum vectors 


n 
around the «® axis. The operator > (x, /2/dx; + p,/\d/dp;) is also an in- 
eae | 
tegral of the motion: 
n N) d 
261 N — |, L,| = 0. 
tee PAG ox, je Fi sp | 


N 

The operator I d/dx, has a striking similarity with the operator of total 
v=1 

linear momentum in quantum theory. It is also closely related with the 


total momentum of the system in classical theory, since: 


n dA n \ 
262 nen 2 
2 dx, | Dr) 
In a similar way we get: 
(263) > | sie SE ii IN le, =(6., 5 x \pi) 2 
| MIT dpi | N 1=1 Jn 


More generally, to any integral F(t; 1,, ...,7,) we can attach an operator ¢(t): 
(264) F (1) 0, = NO RAT ah ny) Ge) are 
We have 


(265) (FÜ), 110; — (ER On; F) 


n 


—(0,, (er Ho) = (9, 
because F satisfies (216). Therefore ¢ satisfies (215) 


6 : 
(266) i 7 AN 


and is an integral. It results from the above results that L, is the ¥-integral 
attached to the energy H,. The integral L, can be easily obtained from (252) 
by the consideration of an infinitesimal time displacement. ¥ is essentially 
the infinitesimal contact transformation determined by F. 
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10. — Expansions in eigenfunctions Ot bie 


The eigenfunctions of L, are the analogue of the wave functions of free- 
particles in quantum mechanics. We want to find the solutions of the 
equation: 


(267) 1,0 =i(P , 6) = v0. 


2m 
Tt is well known that the general solution of the equation 


Bea oF (Pa 
(268) = (F ; r), 


dt 
can be obtained by taking an arbitrary function of six independent integrals 
of the motion of a free particle. We can take the three components of the 
momentum p, two components of the angular momentum M and the time 
dependent integral 


xp P° 
269 shy enr, 
hace) 2E, 2 


Therefore the general solution of (268) is 


= x 


W denoting an arbitrary function. The general solution of (267) is obviously 
the most general solution of (268) such that 


Ir 
(271) i =F, 
i.e. we have 
M "x 
(272) @ — ©(p, M) exp Lin (pe) 


® denoting an arbitrary function. 


A complete system of orthogonal and normalized eigenfunctions of L, is 
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the following 


Dane i 
(273) O(T; Pos Y) = A O(P — Po) EXP sn ; 2 (p'x) "| , 
97 ; | Po _ 14 
(274) “= lu) =) |; RE 


Po and y being two arbitrary constant vectors. O(r; po, y) corresponds to 
a well defined linear momentum py. Since 


- M, : Rx 
CH Lee = x * (YA w) # ee = = [Uy |" (Mo: ¥)(Uo") ; 
al Up, 20H», 


we may write: 


— 


TO) LONTSEB | 


(277) | O*(x; po, y) Or; ps, v') dr = 


3 | I 4 r N NS r 
7 UO(Po D) 0) (Y = ÿ )Aw + 7 ((y ae ( Ve u.) — O(po = P;) öly AY Mic 
S 0 
so that the O*(t; po, y) are orthogonal and normalized. They are a complete 
set because: 


U 


= De 
(278) | Okt; po, Y) AT; po, Y)Apody = e APR): 


: | exp dy = 0(p— p’) 0(x— x’). 


ma [VA m+ tym) 


We shall now expand y(t): 


(279) p(t) = | a(po Y) Ar; Posy) Apo Ay. 


The commutation rules for the a and a* are: 


(250) [po v}, 4*(Por YI]z = Apr PJÔY— Y). 


We can diagonalize Akin: 
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(281) Akin = | dr | a*(po; Y)O*(T; Po: va(p{, Y)L0(T; Poo y') dp, dy dp; dy! = 


= | a*(p,, Ylalpı, vol? v) dp, dy dp, dy! | O*(t; po, Y)O(T: po, Y)dr= 


Po 
= | a*(po, Y)4(Po; Y) : v|dpsay 5 
= 0 


We shall now pass to the more complicated case of particles moving in 
a spherically symmetrical field of force of potential V(r): 


(282) He ee Vien 


Instead of (268) we have now: 


\ dE 0: 
( ) ot 2m ae) a 


The general solution of (284) is an arbitrary function of six independent inte- 
erals of the motion of a particle in the field of potential V(r). We have the 
energy H,, the three components of the angular momentum M and the two 
integrals /, and 7, obtained in the solution of the equations of motion of the 
particle: 


M eur 
(284) Ih, = > | nz x me = = M | Er P 
2 A2 Dr 
zm > + y? | == V( = u ER 


(p is the azimuthal angle around M. p, denotes the radial component of the 
momentum). 


m J pe 


2mr? 


The integrations must be performed with M and H, constant. The general 
solution of (823) is: 


(236) (HAM: Te. lit), 


Y denoting an arbitrary function. The general form of the eigenfunctions 
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of L, is 
(287) © = O(H,, M, I,) exp (il), 


® being an arbitrary function. There are now complications resulting from 
the non uniform integrals J, and J,, especially in the case of an attractive 
potential. The similarity with the situation found in the application of the 
Wentzel-Brillouin-Kramers method is striking. It is interesting to remark 
that the regions of phase space in which p, is imaginary must not be excluded 
because complex values are allowed. We shall examine in detail the problems 
related to the non uniform integrals in a following paper, in connection with 
the passage to quantum theory. 

I thank Professors GEHENIAU and PRIGOGINE for their interest in this 
work. 


RIASSUNTO (*) 


Si applicano all’equazione di Liouville della meccanica statistica classica metodi 
simili a quelli della seconda quantizzazione. Si da un’equazione simile all’equazione 
di Schrodinger di un campo quantizzato. Si dimostra che le regole d’interpretazione 
del tipo quantico del formalismo del campo conducono, trattando le particelle come 
indistinguibili, alle regole della meccanica statistica classica. L’applicazione del metodo 
di seconda quantizzazione di Fock conduce all’introduzione di funzioni d’onda nello 
spazio delle fasi, la densita di probabilita della meccanica statistica classica essendo 
il quadrato del valore assoluto della funzione d’onda nello spazio delle fasi. La scelta 
del segno nelle regole di commutazione degli operatori del campo conduce a funzioni 
d’onda simmetriche o antisimmetriche nello spazio delle fasi. Ad ogni specie di sim- 
metria corrisponde una differente statistica, di Bose o di Fermi, come nella teoria 
quantistica. La statistica di Boltzmann corrisponde a funzioni d’onda nello spazio delle 
fasi prive di condizioni di simmetria. I funzionali d’onda per gli stati con zero par- 
ticelle dei campi nello spazio delle fasi sono analoghi a quelli dei campi quantizzati della 
teoria quantistica. 


. 
()) Traduzione a cura della Redazione. 
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Il meccanismo della scarica nei contatori ad argon. 


UL. Contr e U. FACCHINI 


Laboratori CISE - Milano 


(ricevuto il 7 Ottobre 1952) 


Riassunto. — Sono studiate le curve di moltiplicazione e le forme degli 
impulsi in un contatore riempito con argon puro e con argon + anidride 
carbonica a pressioni da 150 a 1000 mm Hg. La forma degli impulsi in 
argon nella alta zona proporzionale © spiegata assumendo una produzione 
di fotoni nella valanga di Townsend, abbastanza energetic! da estrarre 
elettroni dal catodo. Questi fotoelettroni sono prodotti in un tempo ft, 
dell’ordine di pochi ys a 150 mm Hg e più piccolo alle pressioni più alte. 
Lo studio delle proprietà dei livelli eccitati dell'atomo di argon permette 
di attribuire i fotoni in questione ai primi quattro livelli. Questi fotoni 
sono prodotti in collisioni che portano alla distruzione di questi livelli 
e non sono fotoni di risonanza. I valori di 7,, sono in accordo con. quelli 
ottenuti estrapolando i risultati di Molnar, alle pressioni da noi usate. 
Il processo fotoelettrico & praticamente il responsabile dell’accensione della 
corrente di corona che si ha a moltiplieazioni dell’ordine di 200. Nello 
studio delle miscele A-CO, si osserva un aumento del valore delle mol- 
tiplicazioni alla soglia di divergenza del processo fotoelettrico; questo 
aumento è lineare con la concentrazione di CO,. Un comportamento di que- 
sto tipo si pud spiegare assumendo lo spegnimento dei livelli eccitati dell’ar- 
gon per urto di seconda specie contro le molecole di CO,; dai risultati speri- 
mentali & ottenuto per la sezione d’urto di questo processo il valore di 
~ 600: 10-16 em? che è da ritenersi pertanto solo un ordine di grandezza. 
Ulteriori studi concernono l’effetto di carica spaziale che permette la 
formazione di una breve zona Geiger alla soglia di divergenza del pro- 
cesso fotoelettrico quando la concentrazione di CO, è >5-10-+. Nessun 
altro processo é stato rilevato oltre quello fotoelettrico nelle miscele A-CO, 
fino a concentrazioni di CO, di 5-10-%. Sono descritti infine i compor- 
tamenti del contatore riempito con argon + 1% CO,; argon + 10% CO, 
e con argon contaminato con vapori di mercurio a pressione parziale 
corrispondente alla temperatura ambiente. 


1181 L. COLLI € U. FACCHINI 


1. — Introduzione. 


Lo scopo di questo lavoro & indagare sui processi elementari interessati 
nel funzionamento dei contatori Geiger e sul meccanismo della formazione 
dell’impulso nella zona Geiger, in contatori contenenti argon di elevata purezza 
e miscele argon-anidride carbonica in concentrazioni da 10-4 a 1071, 

Molte ricerche sono state condotte per costruire una teoria che renda conto 
del funzionamento dei contatori e delle loro proprieta (1) (2). Queste ricerche 
hanno messo in chiaro i seguenti punti fondamentali: 


1) La formazione dell’impulso & determinata da un processo primario 
di ionizzazione a valanga nella zona vicine al filo (valanga di Townsend) e da 
processi di ionizzazione secondaria nel gas e al catodo per mezzo dei quali la 
scarica si spande lungo tutto il filo. 


2) Nel caso di miscele di argon e di un vapore organico le esperienze di 
STEVER (*) e di altri (?) (5) hanno mostrato che la propagazione della scarica 
avviene in una ristretta zona lungo il filo del contatore. 

Questo dimostra che il processo responsabile dello spandersi produce una 
ionizzazione secondaria nel gas, vicino al filo. 


3) Nel caso di riempimenti con solo argon, idrogeno, ecc., le stesse espe- 
rienze di STEVER (?) hanno messo in evidenza che lo spandersi della scaric: 
avviene in tutto il contatore, non localizzato sul filo. 

GREINER (°) ha dimostrato che il processo agente nell’idrogeno & dovuto 

fotoni generati nella valanga di Townsend. 

KORFF e ROSE (7) e NAWIJN e coll. (8) hanno mostrato che nella formazione 
dell’impulso era interessato il catodo: l’interpretazione dei loro risultati é stata 
fatta per mezzo del processo fotoelettrico sul catodo. CURRAN @ STROTHERS (°) 
hanno studiato l’emissione fotoelettrica del catodo causata da fotoni prodotti 
nella valanga e la riduzione di questa emissione dovuta ad aggiunta di vapori 
organic. 


(*) S. C. CURRAN e J. D. CRAGGS: Counting Tubes (London, 1949). 

(@) D. H. WırLkınson: Lonization Chambers and Counters (Cambridge, 1950). 

(5) H. G. STEVER: Phys. Rev., 61, 38 (1942); 59, 765 (1941). 

(2) J. Hine J. VY. DunvortH: Mature, 158, 833 (1946). 

(5) E. ALDER, E. BALDINGER, P. Huser e F. METZGER: Helv. Phys. Acta, 20, 


E. GREINER: Zeits. Phys., 81, 543 (1933). 


(°) 

(7) S. Korrr e M. E. Rose: Phys. Rev., 59, 850 (1941). 

(5) A. NAWIJN: Physica, 9, 481 (1942); A. Nawign e J. DE Jone: Physica, 10, 
513 (1943); A. Nawisn e D. MULDER: Physica, 10, 531 (1943). 

(2) S. C. Curran e J. E. STROTHERS: Proc. Camb. Phil. Soc., 35, 654 (1939). 
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L’interpretazione di questi risultati presenta la difficoltà che in generale 
i gas usati avevano un grado di purezza non superiore al 99%, non sufficiente 
per ottenere riproducibilità in misure di questo tipo. 


4) Contemporaneamente allo spandersi lungo il filo della scarica si ha 

la formazione di una carica positiva che costituisce una guaina intorno al filo. 

L’effetto della carica spaziale puö essere diverso in dipendenza dal mecca- 

nismo che determina lo spandersi della carica e quindi dalla miscela di gas 
usata. 


5) Secondo la teoria di MONTGOMERY e MONTGOMERY (1°), e di KORFF 
€ PRESENT (!!) nei contatori à attivo un altro processo secondario di produzione 
di elettroni, il processo di estrazione dal catodo per urto di ioni che, avvenendo 
dopo un tempo grande in confronto alla propagazione della scarica, ha per 
effetto di causare impulsi secondari suecessivi e quindi l’accensione della cor- 
rente di corona. 

Appare evidente che una teoria quantitativa del comportamento di un 
contatore richiede la conoscenza dei vari processi secondari di produzione di 
elettroni relativi alla miscela di gas di riempimento e al catodo del contatore. 

Lo studio dei processi elementari interessati nelle scariche nei gas è state 
condotto estesamente da parecchi autori; una rassegna di questi studi è data 
da PENNING e DRUYVESTEYN (1?) e da LOEB (1%). 

Questi lavori mettono in luce la forte dipendenza di tali processi dalle con- 
dizioni di purezza, pressione e temperatura del gas. 


2. — Metodo di misura. 


Il metodo di misura usato utilizza l’osservazione degli impulsi dati da par- 
ticelle ionizzanti in un usuale contatore cilindrico. 

A tensioni inferiori ad un certo valore Vp (soglia di moltiplicazione pro- 
porzionale) il contatore funziona come una camera di ionizzazione ad elettrodi 
cilindrici. 

A tensioni superiori a Vp si ha la zona di moltiplicazione proporzionale: 
in questa zona un elettrone primario generato nel contatore produce una 
valanga di N elettroni sul filo (valanga di Townsend). N definisce il coeffi- 
ciente di moltiplicazione ed è funzione crescente della tensione applicata (?). 

Ad una certa tensione il contatore entra nella zona di scarica a corona; 


(2) €. G. Monrcommry e D. D. MoNTGomerY: Phys. Rev., 57, 1034 (1940). 

(4) S. A. Korrr e R. D. PRESENT: Phys. Rev., 65, 274 (1944). 

(2) M. J. DRUYVESTEYN e F. M. PENNING: Rev. Mod. Phys., 12, 87 (1940). 

(3) L. Loge: Fundamental Processes of Electrical Discharge in Gases (New York, 
9 
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in opportune condizioni questa zona è preceduta dalla zona dei conti 0 zona 
Geiger. 

Una parte della misura consiste nel rilevare le curve che danno la mol- 
tiplicazione N in funzione della tensione. 

La ionizzazione primaria € ottenuta con particelle « emesse da uno strato 
di polonio collocato nel contatore e con singoli elettroni veloci, prodotti da 
una sorgente di radiazioni y (nel seguito indicati come impulsi y). 

Data l’elevata ionizzazione totale le particelle « si sono usate per rilevare 
le curve di moltiplicazione nella zona «camera di ionizzazione » e nella prima 
parte della zona proporzionale fino a moltiplicazioni di circa 50. 

A moltiplicazioni pit elevate, fino a circa 10 000, a causa della distorsione 
del coefficiente di moltiplicazione dovuta alla troppa ionizzazione data dalle «, 
si sono usati gli impulsi y, i quali corrispondono ad una ionizzazione primaria 
di cirea 50 keV e quindi realizzano condizioni di moltiplicazione lineare sul 
filo. 

Un’altra parte della misura è fatta studiando la forma degli impulsi di 
particelle « e y nelle varie zone e con diverse costanti di tempo nella catena 
elettronica, in modo da ottenere impulsi pit o meno differenziati. 

Dalle forme cosi studiate si pud dedurre lo svolgimento nel tempo della 
scarica e da cid informazioni sulla natura dei processi in essa interessati. 


8. — Apparecchiatura sperimentale. 


Le misure sono state condotte con diverse apparecchiature consistenti in 
uno © pit contatori inseriti su un impiauto adatto alla loro vuotatura e al 
loro riempimevto con le miscele di gas. 

L’impianto cui riferiamo la maggior parte dei dati ottenuti & illustrato 
nella fig. 1. 

Questo è costituito da un contatore (fig. 2) inserito in un circuito metal- 
lico dove il gas vien fatto circolare a termosifone e che serve a provvedere una 
rapida purificazione dell’argon e un rapido miscelamento dei gas in studio. 

Il diametro del contatore, di ottone, è di 60 mm, la lunghezza utile del 
filo di tungsteno è di 200 mm e il suo diametro di 0,2 mm. Il contatore é 
isolato elettricamente dal resto dell’impianto per mezzo di due tubi di vetro. 

Il contatore pud essere scaldato fino a 200 °C per la degassazione. Il tubo 
del contatore € stato degassato prima del montaggio a 600 °C per qualche ora. 

Gli altri elementi del circuito sono: 

1) il forno a circolazione contenente lega di calcio e magnesio al 10% 
che, scaldato a 400-500 °C serve alla purificazione dell’argon (14); 


(4) L. Corzr e U. Faccuini: Rev. Sc. Instr., 28, 39° (1952); Nuovo Oimento, 8, 
892 (1951). 
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2) la camera in cui viene effettuata la misura della velocita di migra- 
zione degli elettroni nel gas. 


alla pompa 


Fig. 1. — Appareechiatura sperimentale: 1) contatore; 2) forno di purificazione; 3) ca- 

mera di ionizzazione; 4) serbatoi; 5) colonna riscaldata; 6) trappola; 7) manometro; 

8) saldature vetro metallo; 9) bombole di argon e di anidride carbonica; 10) vacuometri; 
11) rubinetti. 


Questa misura dä il controllo della purezza dell’argon con una sensibilita 
dell’ordine di qualche 10-4 per impurezze molecolari (7); 
3) la trappola che puö essere raffreddata con aria liquida allo scopo 
di fermare eventuali impurezze condensabili; 
4) il sistema di serbatoi per la miscelazione ; 
5) la colonna riscaldata che permette la circolazione anche escludendo 
il forno. 
Questa operazione è necessaria per mescolare rapidamente la miscela quando 
contiene un gas assorbibile dal forno; 
6) il manometro metallico del tipo Bourdon per pressioni da 10 à 
2000 mm Hg. 
Le guarnizioni per le tenute sono tutte metalliche e interamente metallici 
sono anche i rubinetti; & evitata la presenza di grassi e di vapori di mercurio. 
Il eireuito descritto viene vuotato con una pompa a diffusione ad olio di 
silicone. I] vuoto viene misurato con un vacuometro di Penning. 
L’insieme viene vuotato per parecchi giorni a 10° mm Hg, degasando il 
contatore e il forno. 
Si à usato argon di diversa provenienza € di purezze iniziali 98 e 99,9%. 
In ogni caso si sono ottenuti risultati riproducibili. Per le miscele si è usata 
anidride carbonica di purezza 97%. 
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L’apparecchiatura elettronica con cui si studiano gli impulsi é illustrata 
nella fig. 3. 

Il contatore @ contenuto in una 
scatola metallica che fa da schermo elet- 
trostatico. Il filo del contatore è connesso 
attraverso una convenzionale catena di 
amplificatori ad un oscillografo Dumont 
tipo 248 la cui traccia pud essere co- 
mandata a mezzo di un sincroscopio. 

Il tempo di salita della catena e di 
0,6 us. Il alcuni casi € stato utile usare 
una catena pitt rapida (0,2 us). 

Le misure di ampiezza degli impulsi 
sono fatte usando una costante di tempo 


di ingresso molto maggiore dei tempi di 
salita degli impulsi (RC circa 200 ps). 
Fra il preamplificatore e l’amplificatore 
viene in questo caso inserita una linea 
di riflessione la quale ha lPufficio di ta- 
gliare la discesa dell’impulso dopo un 
tempo variabile da 1,2 a 6us. In tal 


modo si misura solo la prima salita del- 
Vimpulso che corrisponde alla collezione 
degli elettroni della valanga di Townsend 


Fig. 2. — Contatore: 1) Catodo diot- © Si esclude una componente dell’impulso 
tone; 2) saldatura ad argento; 3) ko-  ritardata nel tempo e dovuta, per esem- 
var; 4) vetro; 5) filo di tungsteno.  pio, ad elettroni prodotti per effetti se- 

condari. Questo e importante per la 
misura della moltiplicazione sul filo quando il contatore funziona neil’ultima 
parte della zona proporzionale dove questi processi secondari sono rilevanti. 


Fig. 3. — Apparecchiatura elettronica: 1) contatore; 2) schermo elettrostatico; 3) pre- 
amplificatore; 4) amplificatore; 5) oscillografo; 6) alta tensione; 7) filtro. 
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La forma dell’impulso si studia osservandone lo sviluppo nel tempo con 
il sincroscopio ed usando una costante di tempo d’ingresso variabile da 1 a 
200 us che provvede un’opportuna differenziazione degli impulsi. 

La forma dell’impulso puo essere osservata per una durata variabile da 
qualche us a 30 ms usando diversi tempi di sparpagliamento della traceia 
dell’oscillografo. 


4. — Misure sull’argon. 


Le curve di moltiplicazione N(V) a varie pressioni di argon sono mostrate 
nella fig. 4. 

Durante le misure l’argon è fatto circolare nel forno di purificazione. Le 
misure in tali condizioni sono riproducibili nel tempo. Non si sono ottenute 


600 1'000 1'500 2'000 2'500 3000 V 


Fig. 4. — Curve di moltiplicazione per argon: I) p = 150 mm Hg; II) p = 500 mm He; 
III) p = 1000 mm Hg. 


variazioni mettendo dell’aria liquida intorno alla trappola e lasciandovela 
alcune ore. | 

Gli impulsi « e y osservati al sincroscopio hanno la forma usuale (fig. 5) 
fino a moltiplicazioni di circa 50, 

A partire da queste moltiplicazioni ogni impulso, è seguito da una succes-' 
sione convergente di impulsi distanziati di un tempo 7, l’uno dall’altro ed! 
avente ragione crescente con N. i 

Questi impulsi si mettono in evidenza con una costante di tempo dell’or- 
dine di 1 us, cioè inferiore alla loro separazione, che & dell’ ordiue di 10 us 
per tutte le pressioni usate. 
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Fig. 5. - Impulso di particella « nella zona 
proporzionale a bassa moltiplicazione. Se- 


enali di taratura = 10 us ROUTE: 


e 
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Il tempo di salita degli impulsi suc- 
cessivi di una stessa successione si al- 
lunga con l’aumentare del numero d’or- 
dine dell’impulso; di conseguenza di- 
minuisce la nettezza della separazione 
fra gli impulsi. 

Questo effetto 
grande quanto più la pressione @ ri- 
dotta. A pressioni di 150 mm Hg non 


diventa tanto pit 


si vede altro che il secondo impulso 
seguito da una discesa continua. 
Raggiunto un valore della moltipli- 


cazione. dell’ordine di 200 nel conta- 


tore si innesca la corona stabile, dopo una brevissima zona instabile ampia 


solo qualche volt. 


Fig. 6. — Impulsi di particelle & in argon 
alla soglia dicorona. 

N = 210; Se- 

RC = 1 us. 


Pressione = 1000mm He; 


enali di taratura = 1 us; 


Fig. 7.— Impulsi di particelle + in argon alla 
soglia di corona. L’allungamento del tempo 
di salitasi osserva negli ultimi fotoimpulsi. 
Pressione = 500 mm Hg; N = 250; Se- 


enali di taratura = 1 us; RC = 1 us. 


Fig. 8. — Impulsi di particella « in argon 
alla soglia di corona. 

Pressione = 150 mm He; N = 180; Se- 

JRO! = IN AES 


gnali di taratura = 1 us; 


Fig. 9. — Impulso di particella « in argon 


alla soglia di corona. 
Pressione = 150 mm Hg; N = 180; Du- 
rata della traccia = 100 us; RC = 10 us. 
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Le forme degli impulsi alla soglia della corona sono rappresentate nelle 
fige. 6, 7, 8 e 9. La serie degli impulsi successivi resta sempre conver- 
gente anche nel caso di impulsi y. 

Risultati del tutto simili sono stati ottenuti con un altro impianto costi- 
tuito da 4 contatori in parallelo con catodi di rame o di ottone e un complesso 
di circolazione e di purificazione del tipo di quello gia descritto. Soltanto i 
valori di N in soglia potevano differire di un fattore 2. 


5. — Discussione dei risultati sull’argon. 


I risultati ottenuti con l’argon ed in particolare la forma degli impulsi a 
moltiplicazioni elevate si possono interpretare per mezzo dei fenomeni che 
governano la produzione e il decadimento degli atomi eccitati nella valanga 
di Townsend. 


A) I fotoimpulsi. — Come & noto, un elettrone liberato nel contatore, 
arrivato in vicinanza del filo, dove il campo è elevato, pud ionizzare per urto 
eli atomi di argon, dando luogo alla valanga di Townsend. Un elettrone liberato 
nel contatore produce un impulso contenente N elettroni, la cui ampiezza è 
proporzionale ad N. 

Consideriamo gli atomi eccitati prodotti nella valanga; sia x il rapporto 
fra il numero di atomi eccitati ed il numero di ionizzazioni prodotte nella 
valanga. II numero di atomi eccitati per elettrone iniziale © Nz. 

Si pud assumere che i livelli eccitati abbiano una vita media che generi- 
camente indichiamo con 7, e che nei processi interessanti il loro decadimento 
nascano dei fotoni abbastanza energici da estrarre elettroni per effetto foto- 
elettrico dal catodo di ottone (soglia fotoelettrica 4-5 eV). 

Indichiamo con ß il numero di fotoni di energia sufficiente per il. processo 
fotoelettrico che in media vengono emessi per atomo eccitato, con G la pro- 
babilita che questi fotoni arrivino al catodo senza essere distrutti; sia inoltre 
K l’angolo solido utile e y il rendimento fotoelettrico medio del catodo per i 
fotoni interessati. 

in un tempo dell’ordine di 7, gli NafG@ fotoni producone NapGKy foto- 
elettroni al catodo. 

I fotoelettroni arrivano al filo dopo un tempo 7, (tempo di transito nel 
contatore) e si moltiplicano per N dando un nuovo impulso di ampiezza 
(NopGKy)N ad una distanza circa T, + T;, dal primo. Questo processo Si 
ripete e trascurando la variazione di N dovuta all’accumularsi di ioni sul filo, 
produce una successione di impulsi di ragione (Naß@Ky) crescente con la 
tensione. 

Alla tensione per cui 


(1) NopGkKy =1 
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si ha, in questa approssimazione, la soglia di divergenza del processo; chia- 
miamo N,, moltiplicazione di soglia, il valore di N che soddisfa la (1), e in- 
dichiamo con V, la tensione corrispondente. 

In questo modo si possono interpretäre le successioni di impulsi descritte 
nel paragrafo 4, il loro comportamento con la tensione e l’accendersi della 
corona. 

I singoli fotoimpulsi della successione possono osservarsi ben separati solo 
se il tempo 7, & molto piccolo rispetto al tempo 7,. 

Se T, non è trascurabile rispetto a 7, si ottiene ovviamente un allunga- 
mento dei tempi di salita dei fotoimpulsi successivi ed una peggiore separa- 
zione fra di essi. Se 7’, poi è dell’ordine di grandezza di 7, non c’è nessuna 
separazione fra i fotoimpulsi. 

Dall’esame delie successioni di impulsi nell’argon si pud osservare che 7, 
non & trascurabile rispetto a 7,. 

Si nota infatti un allungamento dei tempi di salita per i successivi impulsi 
di una stessa serie (fig. 7); questo allungamento diventa pit rilevante al dimi- 
nuire della pressione finchè alle pressioni più basse la separazione dei fotoimpulsi 
non & più realizzata (fig. 8). 

Questi fatti indicano che il tempo 7’, cresce con il diminuire della pressione. 

E interessante ottenere dall’osservazione degli impulsi una stima dei valori 
possibili per 7’, . 

Per far cio si pud misurare 7, e dedurre 7, dalla relazione approssimata 
TT, = Te. 

I valori di 7, sono ottenuti misurando la distanza fra i massimi del primo 
e del secondo impulso della successione, usando una RC di 1 us circa. 

les puod essere calcolato con una certa precisione conoscendo i valori delle 
velocità di migrazione W- degli elettroni (1) (5). Abbiamo per maggior sicu- 
rezza rilevati i valori di W- con la camera illustrata in paragrafo 3, contem- 
poraneamente alle misure di 7%. 

Per valori di #/p (E campo elettrico, p pressione) più elevati di 1,5 volt 
em! mm Hg-! abbiamo usato 1 valori di NIELSEN (1). 

I valori di T,, T, e n— T, = T,, sono raccolti nella tabella I. 

I valori di 7, rappresentano solo l’ordine di grandezza del tempo di per- 
manenza dei livelli eceitati nel gas e nel caso delle pressioni più elevate sono 
soltanto dei limiti superiori. 

L’assunzione di un tempo 7, di decadimento dell’ordine di grandezza 
indicata nella tabella rende conto della forma degli impulsi e come si vedrä 
in seguito fornisce una semplice interpretazione dei risultati nelle miscele 


(15) R. A. NIELSEN: Phys. Rev., 50, 950 (1936). 
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Una causa di allungamento dei tempi di salita negli impulsi successivi 
potrebbe essere dovuta ad un’azione del tempo di sviluppo della valanga elet- 
tronica ed anche ad una diffusione geometrica degli impulsi lungo il contatore. 
Si puö ritenere che questi effetti siano trascurabili poichè in miscele di A + CO, 
(vedi n. 7) si possono osservare parecchie decine di fotoimpulsi ben separati 
e nessun allungamento sensibile dei tempi di salita. 

In queste condizioni azioni geometriche o di valanga avvengono in modo 
analogo a quanto accade nell’A puro e non portano un effetto sensibile. 


B) Il decadimento dei livelli eceitati. — Ci proponiamo ora di analizzare 
più in dettaglio il fenomeno di decadimento dei livelli eccitati nella valanga, 
Consideriamo a questo scopo lo spettro dei livelli dell’atomo di argon che 

è riportato nella tabella II (15) (7). 

Dei primi quattro livelli eceitati il ®P, e 1’1P, sono livelli di risonanza e 
il *P, & ?P, sono livelli metastabili. Tutti questi livelli, come vedremo dopo, 
hanno una certa durata nel gas (15). 

Dei livelli superiori solo pochi possono decadere direttamente allo stato 
fondamentale emettendo fotoni di 14-15 eV; la più parte decade ad uno dei 
primi quattro livelli con fotoni di energia inferiore a 4 eV. La vita media per 
queste transizioni si ritiene generalmente dell’ordine di 10-7-10-°s molto infe- 
riore alla durata dei primi quattro livelli. 

Notiamo che i fotoni emessi in queste transizioni non hanno energia suffi- 
ciente per produrre estrazione fotoelettrica dal catodo e quindi li escludiamo 
dalla discussione. 

Si ha che per qualsiasi distribuzione iniziale dell’eccitazione, dopo un tempo 
corrispondente a qualche transizione radiante, sono presenti nel gas quasi 
esclusivamente livelli di risonanza e metastabili. 

Possiamo limitarci a considerare il comportamento di questi soli livelli. 

I livelli *P, e 1P, possono decadere allo stato fondamentale emettendo 1 
fotoni di risonanza di 11,62-11,82 eV. 

I fotoni di risonanza come è noto subiscono parecchie volte la cattura e 
la riemissione da parte degli atomi non eccitati del gas. 

HOLSTEIN (1°) ha sviluppato una teoria di questo processo e caleola che 
il fattore di imprigionamento 1/g per un cilindro di qualche em di diametro 
(numero di volte che in media un fotone & catturato e riemesso prima di rag- 
giungere le pareti del cilindro) puö essere dell’ordine di qualche centinaio. 


(6) W. Dr Groot e F. M. Pennine: Handb. d. Phys., 28/1 (1933). 

(17) C. E. Moore: Atomic Energy Levels, Circular NBS n. 4467 (1949). 

(3) A. C. G. MiTcHELL e M. W. ZEMANSKY: Resonance Radiations and Excited 
Atoms (New York, 1934). 

(9) ©, HozsteiN: Phys. Rev., 72, 1212 (1947); 83, 1159 (1951). 
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Si pud ritenere secondo questa teoria che la durata di imprigionamento 
sia dell’ordine di qualche 10-5-10-1$, molto più grande della vita media propria 
dei livelli risonanti. 

HOLSTEIN ottiene inoltre che il fattore di imprigionamento € indipendente 
dalla pressione in una zona di pressione in cui la forma della riga di emissione 
sia determinata dall’allargamento di pressione. Non si conoscono misure di 
1/g nel caso dell’argon. 

I livelli metastabili hanno una vita media propria molto lunga; il processo 
che porta al loro decadimento è dovuto in generale agli urti che subiscono 
con gli atomi non eccitati del gas. 

In recenti misure MOLNAR (?°) (2!) e PHELPS (??) hanno studiato la durata 
di questi livelli a pressioni da 1 a 10 mm Hg. Questi autori ottengono vite 
medie dell’ordine di qualche ms e mostrano che la probabilità per secondo 
di distruzione nel gas di questi livelli cresce cirea con il quadrato della pressione. 
Questa distruzione puo essere attribuita secondo MOLNAR ad una sorta di 
urti che probabilmente conduce alla formazione di molecole eccitate instabili 
fra atomi metastabili ed atomi neutri. 

Processi di questo tipo sono stati invocati anche da altri autori (*) per 
Spiegare parecchi fenomeni connessi con le scariche nei gas per quanto le prove 
chiare della loro esistenza siano state date per ora solo nel caso del Hg (74) e 
dell’He (?5). 

In questo processo & ragionevole pensare che possano prodursi dei fotoni 
abbastanza energici, i quali essendo fotoni di risonanza dell’atomo di argon 
arrivano al catodo senza essere assorbiti. 

Una teoria adeguata del decadimento dei livelli di risonanza e metastabili 
non è stata finora sviluppata ed i dati sperimentali sono piuttosto scarsi; il 
decadimento € inoltre compheato dal fatto che possono esistere dei processi di 
urto in cui, a causa della piccola differenza di energia che c’é fra i livelli, del 
Vordine dell’energia cinetica termica degli atomi, si possono avere delle tran- 
sizioni dall’uno all’altro di essi. 

Possiamo perd fare le seguenti osservazioni: alle elevate pressioni il deca- 
dimento di tutti i primi livelli pud essere considerato globalmente tenendo 
in conto le transizioni per collisione. 

Si hanno le due possibilita: 

1) i livelli decadono attraverso l’emissione di fotoni di risonanza che 


(2°) J. P. Mornar: Phys. Rev., 83, 933 e 940 (1951). 
(21) J. P. MOLNAR: Comunicazione privata. 

(2) A. V. Puæzps: Phys. Rev., 82, 567 (1951). 

(23) E. W. Pike: Phys. Rev., 53, 921 (1938). 

(24) A. Q. McCouUBREY: Phys. Rev., 84, 107 (1952). 
(5) R. MEYEROTT: Phys. Rev., 70, 670 (1946). 
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raggiungono il catodo con una probabilita circa indipendente dalla pressione 
che é rappresentata da g/t* essendo 7* la vita media propria dei livelli risonanti; 

2) per mezzo del processo messo in luce da MOLNAR con probabilita 
crescente con il quadrato della pressione. 

RP ragionevole ritenere che quest’ultimo processo divenga predominante 
alle elevate pressioni e quindi il solo a regolare il decadimento di tutti i livelli. 

Ne segue che a queste pressioni i fotoni di risonanza non raggiungono pra- 
ticamente il eatodo e l’effetto fotoelettrico è provocato dai fotoni generati 
nelle collisioni del tipo descritto da MOLNAR. 

In questa descrizione f rappresenta il numero medio di fotoni « energici » 
emessi in questo processo, y il loro rendimento medio al catodo; @ nell’argon 
si pu ritenere 1. 

Abbiamo, in queste considerazioni, estrapolato i risultati di MOLNAR otte- 
nendo i valori (T,,), elencati nella tab. I che sono in accordo con i nostri ordini 
di grandezza. 


TABBLLA I. — Valori del tempo di decadimento Ty, in argon. 
Pressione Ve | I fi LT, =D le, (Ty) 
(mm Hg) volt (us) (us) (us) (us) 

150 | 1 200 | 10,8 se 1 6,3 +1 | 4,5 +2 4,1 
300 1 640 10,7 EU PA UNS 620 1.1 
500 2 070 10,2 +1 DEA 26 0,4 
1000 3 000 10,5 + 1 8,8 +1 Ze Weil 0,1 


Si pud dire che i nostri risultati danno una conferma dell’esistenza di un 
processo che fornisce una rapida distruzione dei livelli in discorso nell’argon 
ad elevate pressioni. 

Una discussione più dettagliata intorno ai processi di decadimento si trova 
in appendice. 


C) Soglia di corona. — Consideriamo Vaspetto degli impulsi in soglia di 

corona. 

Se la corona è dovuta al solo processo fotoelettrico, la soglia di corona € 
raggiunta al valore di moltiplicazione N, definito dalla (1). 

La ragione della successione di fotoimpulsi in questa ipotesi diventa a ter- 
mini costanti. | 

Esaminando agli inizi della corona la forma degli impulsi « e y con RC 
di 1 us si vede che la suecessione dei fotoimpulsi ha andamento convergente. 

I piccoli valori di RC usati, inferiori al tempo di salita degli impulsi, produ- 
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cono una riduzione del valore della ragione della successione dovuta al fatto 
che i tempi di salita vanno allungandosi e quindi gli impulsi susseguenti sono 
più diminuiti in ampiezza. 

Tenendo conto in modo approssimativo di questo fatto si pud valutare 
che la ragione della successione in soglia non Sia inferiore a 0,7-0,9. 

Non si puo escludere pertanto che all’accensione della corona contribui- 
scano anche altri processi che potrebbero non apparire nella struttura degli 
impulsi perché non osservabili nelle condizioni sperimentali usate; si pud perd 
dire che questi processi non portano un contributo rilevante all’accensione 
della corona alle pressioni da noi usate (?f). 


6. — Miscele A+ CO,. Spegnimento degli atomi eccitati. 


Le miscele argon + anidride carbonica sono state studiate a pressioni da 
150 a 1000 mm Hg e con concentrazioni di CO, da 10-4 a 10-'. Descriviamo 
qui i risultati ottenuti con concentrazioni inferiori a 5-107’. 


1% 


a et 2 
1'000 1500 2000 2500 3000 Volt 3'500 


He 


Fig. 10. — Curve di moltiplicazione per miscele argon-anidride carbonica: 
I) p= 150 mm Hg; II) p = 500 mm Hg; III) p = 1000 mm Hg 
© argon; @ argon + 3:10 CO,; O argon + 17-1074 CO,. 


(5) Notiamo che recentemente LAUER (Journ. of Appl. Phys., 23, 300 (1952)) ha 
messo in evidenza la presenza di un processo di estrazione ione-catodo nell’argon usando 
catodi di nichel ben degassati. 
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A) Risultati sperimentali. — La fig. 10 rappresenta alcune curve di mol- 
tiplicazione. 
Si nota che nel tratto iniziale le curve a pari pressione si sovrappongono 
Puna all’altra entro gli errori di misura che si possono ritenere del 20%. 


Gli impulsi « e y ad elevati valori della moltiplicazione sono costituiti da 
una successione di impulsi del tipo descritto nel caso dell’argon (figg. 11 e 12). 


Fig. 11. — Impulso y in argon+5-10~4 CO, Fig, 12.- Impulso yin argon+ 17-107? CQ, 
poco sotto la soglia Geiger. poco sotto la soglia Geiger. 
Pressione = 500 mm Hg; Segnali di ta- Pressione = 150 mm Hg; Segnali di ta- 
ratura = 1 us; RC = I us. ratura = I us; eC = 1 us: 


La ragione della successione € al solito crescente con il crescere della mol- 
tiplicazione fino a raggiungere la divergenza. 

Per queste miscele abbiamo osservato un valore di tale ragione uguale 
a 1 ed anche più grande di 1 senza avere una corona stabile. 

Una discussione dettagliata del comportamento del contatore in soglia di 
divergenza è presentata nel paragrafo 7; ci interessa qui studiare il valore N, 
della moltiplicazione in soglia di divergenza, intendendosi per soglia di diver- 
genza quel valore V, della tensione a cui la ragione della successione & 1. 

I valori ottenuti per N, sono rappresentati in fig. 15. 

Si vede che N, in funzione della concentrazione di CO, e a pressioni di 
150-500-1000 mm Hg cresce circa linearmente con X e che la pendenza delle 
rette N,(X) è più grande per i bassi valori della pressione. 


B) Discussione. — Dal fatto che a parità di tensione applicata non vari 
praticamente il valore della moltiplicazione con la percentuale di CO, come 
si vede per esempio nelle curve di fig. 10, si puö dedurre che la struttura della 
valanga di Townsend non è praticamente influenzata da piecole percentuali 
di CO,. Ne segue che «, numero di atomi eccitati generati nella valanga per 
elettrone, si pud ritenere indipendente dalla concentrazione X della COs. 

Le successioni di impulsi che si osservano in tutte le miscele e pressioni 
descritte in questo paragrafo possono interpretarsi come nel caso dell’argon 
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puro, come successioni di fotoimpulsi: la ragione della successione @ data da 
N(V)«fGK(1 — P)y. 

P indica la probabilita che i livelli eccitati di risonanza e metastabili si 
distruggano per urti con le molecole di CO,, di tipo diverso da quelli sull’argon, 
senza dar luogo a fotoni efficienti per l’estrazione fotoelettrica. 

Il valore di soglia N, è determinato da 


(2) NaPKG(1 — Py =1., 


Le quantita considerate possono dipendere in generale dalla pressione € 
dalla concentrazione di CO,. 

Non si sa con precisione quale sia la dipendenza di y, rendimento fotoelet- 
trico del catodo, dalla concentrazione X, pertanto riteniamo che non sia for- 
temente variabile con essa; x, 6, K sono da ritenersi costanti. 


TABELLA IT. — Livelli eccitati del’ atomo di argon. 
Energia 
Termini di eccitazione Tipo del livello 
eV | 
| 7 
SH) 0 | Stato fondamentale 
= ee | | 
| | 
32 11,545 | Metastabile | 
8P, 11,620 | Risonante | 
sp, 11,720 | Metastabile | 
1p, 11,825 | Risonante | 
IS 124 10) 12,7 — 13,3 
| 1,3P | -14,0-14,9 | 
| LP; DF 166-143 7 
[ 15,68 LAS | 
| 15,86 | onizzazione | 


La forte variazione di N, con X puo essere attribuita o ad una riduzione 
della probabilità G che hanno i fotoni di raggiungere il catodo, o ad un au- 
mento con X della probabilità P che i livelli eccitati di argon si distruggano 
per urto senza dar lucgo a fotoni duri. 

G potrebbe essere ridotto a causa dell’assorbimento dei fotoni da parte 
della CO,. Si vede facilmente che questo processo deve dare per G e quindi 
per N, una dipendenza esponenziale della quantita di CO, contenuta nel 
contatore. 
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I risultati sperimentali non sono in accorde con questa conclusione, in 
quanto N dipende linearmente da X e la dipendenza dalla pressione totale 
è in senso opposto a quella prevista. Questo porta ad escludere l’assorbi- 
mento dei fotoni come processo importante per la riduzione di N, Si puo 
dire che @ — 1 anche in presenza di CO,. 

Per la variazione di P vanno considerati gli urti di seconda specie che gli 
atomi eccitati di argon fanno sulla CO, durante la loro permanenza nel gas. 

Processi di questo tipo sono stati studiati da parecchi autori (7%). In gene- 
rale si ha che l’energia di eccitazione dell’atomo eccitato & ceduta alla mole- 
cola urtata che rimane dissociata o eccitata © ionizzata. 

Nel presente caso la reazione avviene probabilmente secondo lo schema (*?) 


COs A = OA 


Tl CO* rimane eccitato per 6 eV. 

L’energia della reazione é circa 11,5 eV ed & molto vicina a quella dei livelli 
eccitati dell’argon che consideriamo. Questo pud favorire a causa della risonanza 
la probabilita del processo (?*). 

In prima approssimazione consideriamo che gli atomi eccitati si trovino 
in un solo livello e sia N* il numero di questi atomi presenti nel contatore. 
L’equazione di decadimento si pud scrivere: 


N* & 
(3) dN*= — T. dt N *nco,ov dt ; 


fo 


-- (N*/T,) rappresenta il decadimento esponenziale che tiene conto del pro- 
cesso di distruzione in argon 


— N*Nco, ov rappresenta la distruzione dei livelli per urto su CO, 


Noo, = MX = numero di molecole di CO, per cm? 


NA — numero di atomi di argon per cm? 
o __ sezione d’urto del processo di diseccitazione sulla CO, 
® — velocità media. 


(27) H. J. Henning: Ann. der Phys., 18, 529 (1932). 
(28) N. F. Mort e H. 8. Massey: The Theory of Atomic Collisions (Oxford, 1949), 
Cap. XII. 
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La frazione di atomi eccitati 
che viene distrutta in argon pro- 
ducendo i fotoni utili e data evi- 
dentemente da 


(4) HP u 


e si ha infine dalla (2) 


1 ei 
5 N= (1 ny OGL e 
(5) ; apy | svoTl’;,X) 
Nelle approssimazioni fatte, la 
formula (5) dä per N, un anda- 
mento lineare con X. Tenuto conto 


degli errori di misura, questa for- 

mula rappresenta soddisfacente- 
mente i risultati sperimentali. 

Dall’equazione (5) e dai risultati 

oS pot UE LS spossond, esscro de dott 

Fig. 13. — Valori di N, in funzione di X,con- valori di 7,6 alle varie pressioni, 

centrazione di CO, : wpressione—150 mm Hg; essendo noti gli altri parametri. 


® pressione = 500 mm Hg; —e- pressio- Questi valori sono rappresentati 
ne — 1000 mm Hg 


nella tab. III, e sono decrescenti 
proporzionalmente al quadrato della 
pressione. Questo sembra indicare che anche nel nostro argon la distruzione 
dei livelli in discorso sia regolata dallo stesso processo messo in luce 
da MOLNAR. 


TAaBerrA III. — Valori di oT,, e di o dedolli dalle ewrve di fig. 13. 
Pressione noo, NV oT,, (2 
(mm Hg) ö M Git Sat cm? s em? | 
> | 
150 33 820 0,18 : 1024 187,8: 10721 | 458: 10716 | 
500 14 507 0,61 : 1021 23,8 : 10721 595-1076 
1 000 7720 | 1,22. 10° 63311082: 633: 10716 


Utilizzando i valori di (7, )„ di tab. I, ottenuti estrapolando i valori di 
MOLNAR, si ottengono i valori della o riportati nell’ultima colonna della tab. III. 
Informazioni sul diverso comportamento dei livelli risonanti e metastabili 
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dovuto a diversi valori delle sezioni d’urto del processo di distruzione e alla 
probabilita di transizione dall’uno all’altro di questi livelli potrebbero cercarsi 
studiando in maggior dettaglio le curve di N,(X) specialmente nel primo tratto 
dove si realizza lo spegnimento del grosso degli atomi eccitati. 

I limiti della precisione della nostra misura e le scarse informazioni che si 
hanno a disposizione per quel che concerne questi processi, non ci consentono 
di perfezionare i nostri risultati. Cosi pure rimane incerto il significato preciso 
della o data nella tab. III. 


D) Purezza dell’argon. — Si vede dalla fig. 13 che il valore di N, è sen- 
sibile a concentrazioni di CO, dell’ordine di 10-* e meno. 
Da questo si pud dedurre che l’argon usato per le miscele contiene meno 
di 10-4 di CO, o impurezze che abbiano lo stesso effetto di spegnimento dei 
livelli eccitati. | 


7. — Miscele A + CO,. Formazione della zona Geiger. 


1) RISULTATI SPERIMENTALI. 


Vogliamo considerare in questo paragrafo più in dettaglio la forma degli 
impulsi ottenuti con le miscele A + CO, per concentrazioni di CO, < 510-3 
ed in particolare l’aspetto degli impulsi alla soglia di divergenza del processo 
fotoelettronico. 


A) Forma degli impulsi. — Gli impulsi « e y ad elevati valori della mel- 
tiplicazione sono costituiti da successioni di impulsi fotoelettronici di ragione 
crescente con la tensione. 

Dall’esame di queste successioni si ottiene: 

1) la distanza T fra l’uno e l’altro degli impulsi successivi decresce con 
il crescere della percentuale di CO,. I valori ottenuti per T sono illustrati 
nella tab. IV per pressioni di 500 e 1000 mm Hg; 

2) la nettezza con cui si separano i fotoimpulsi successivi & tanto pit 
*grande quanto più grande é la percentuale di CO,. Con 5-10™ di CO, à pres- 
sioni di 500 e 1000 mm Hg si distingue nettamente qualche centinaio di foto- 
impulsi a seguito di un impulso y primario. 

A pressioni di 150 mm Hg si distinguono nettamente i foteimpulsi con 
concentrazioni di CO, più grandi di 10°. 


B) Soglia di corona: miscele con concentrazioni di CO, < 51074. — In queste 
miscele raggiunto il valore della moltiplicazione di soglia di divergenza, 
(N= N,), si ha bruscamente la zona di corona instabile. 
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TABELLA IV. — Valori del tempo T, fra à foloimpulsi in A e in miseele A + CO,. 
Pressione Concentrazione T 
(mm He) di CO, (LS) 

1 000 0 10,5 
) 10-4 8,9 
» 3.104 8,1 
» 9-10-4 6,3 
» Po Or 4,4 
500 Ö 10,2 
» 107% 8,3 

) oll Ome 7 

) 9-10-4 6 
) 17-10-4 4,1 


Gli impulsi y alla soglia sono costituiti da una successione di foteimpulsi 
distinti molto lunga che ne comprende qualche centinaio. La successione ha 
dapprima andamento convergente con ragione prossima ad 1 e procede poi 
a termini pressapoco costanti. 

Ogni impulso si pud considerare come un burst di corona il cui spegni- 
mento & di natura statistica. 

Pochi volt più sopra si ha la corona stabile. 


C) Formazione della zona Geiger: miscele con concentrazioni di CO, fra 


Fig. 14. — Impulso y nella zona Geiger. Fig. 15. — Impulsi y sovrapposti nella 
Ar + ~ 107% CO,. zona Geiger. A + 17-1074 CO,. 
Pressione = 760 mm Hg; Durata della Pressione = 1000 mm Hs; Distanza fra 
traccia — 50)us; AO —= Vous: i fotoimpulsi = 4,4 us; RC = 1 us. 


5-10-4 e 5:10. — A concentrazioni di CO, maggiori di 5-10”? si raggiungono 
valori di N > N,; senza che si accenda la corona instabile. 
Gli impulsi y sono del tipo di fig. 14. 
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La successione ha un andamento dapprima divergente e poi convergente 


fino a spegnersi. La struttura di questi gruppi di fotoimpulsi e indipendente 
dalla ionizzazione iniziale (fig. 15). Per questo il gruppo di fotoimpulsi si pud 


Fig. 16. — Dettaglio di bursts nella zona Fig. 17. — Dettaglio di burst nella zona 

di corona instabile. A + ~ 107% CO,. di corona instabile. A + ~ 1073 CO,. 

Pressione = 760 mm Hg; Durata della Pressione = 760 mm Hg; Durata della 
imaccia == LOO uss “OC == IN us. traccia = 100 us; RC =1 us. 


chiamare impulso Geiger e la zona di tensione, in cui si ha questo comporta- 
mento, zona Geiger. 

In queste miscele la soglia di divergenza V,, definita nel paragrafo 6, corri- 
sponde alla soglia della zona Geiger. 

L’ampiezza di questa zona è di qualche decina di volt. Dopo di essa si ha 
una zona di corona instabile formata di bursts spenti statisticamente. Le fige. 16 
e 17 rappresentano particolari di bursts di corona in cui la struttura a impulsi 
fotoelettronici è messa in evidenza con RC di 1 us. 

L’ampiezza dei singoli fotoimpulsi oscilla periodicamente. 

Alla fine della zona instabile, ampia pochi volt, si accende la corona stabile, 
con andamento oscillante. 

In comportamento di questo tipo @ stato osservato anche con parecchi 
contatori di diametro da 2 a 15 em a pressioni da 200 a 1500 mm Hg e con 
miscele di A 10 di CO, (29) (figg:-18, 19, 20 e=2%). 

E ante notare che nelle zone Geiger e nelle zone di corona instabile, 
mentre gli impulsi y hanno l’aspetto descritto, le particelle « danno luogo ad 
impulsi convergenti e con ragione dell’ordine di 0,5 (fig. 22). 


(2) L. Cort, U. Faccninı e E. Garri: Phys. Rev., 80, 92 (1950); L. Coruı, U. 
FaccHInı e E. Garrr: Phys. Rev., 84, 606. (1951). 
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Fig. 18. — Impulso y nella zona Geiger. Fig. 19. — Burst all’inizio della zona di 
A + m 1073 CO,. corona instabile. A + ~ 1073 CO,. Dia- 

a à » D we = € . ” 12 

Diametro del contatore — 90 mm; Pres- metro del contatore — 20 mm; Pres 
sione 760 mm Hg; Tempo fra i foto- sione = 760 mm Hg; Segnali di tara- 
impulsi = 12 us. . tura = 100 us; RO = 10=us. I gruppi 


di impulsi fotoelettronici appaiono come 
singoli impulsi per via dell’elevata RC. 


\ ARENA 


Fig. 20. — Burst alla fine della zona di Fig. 21. — Bursts nella zona di corona 
corona instabile. A + ~ 1073 CO,. Dia- instabile. A + ~ 107 CO,. 
metro del contatore = 20 mm;  Pres- Diametro del contatore= 130 mm; Pres- 
sione 760 mm He; Segnali di tara- sione — 760 mm Hg; Durata della traceia 
tura = 100 us; RO 10 us. I gruppi = 30 us; RO = 200 us. 

di impulsi fotoelettronici appaiono come 


singol impulsi per via dell’elevata RC. 


2. — DISCUSSIONE DEI RISULTATI SPERIMENTALI. 
A) Distanza dei fotoimpulsi. — Le successioni di impulsi descritte sono 


da ritenersi, come nel caso dell’argon, dovute al processo fotoelettronico sul 
catodo. 


La distanza di tempo 7 fra gli impulsi fotoelettronici successivi € la somma 
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di due tempi T,e T,; T, tempo di decadimento dei livelli eccitati nella miscela 
T, tempo di transito dei fotoelettroni nella miscela. 

Da quanto visto nel paragrafo 6 si pud dedurre che 7, è fortemente ridotto © 
per effetto degli urti di spegnimento 
dei livelli eccitati dell’A su CO, (*); 
T, è pure ridotto per effetto delle CO, 
in quanto questo gas molecolare au- 
menta la velocità di migrazione degli 
elettroni in campo elettrico (1). 

Percentuali di 0,1% di CO, rad- 
doppiano cirea il valore di tale velo- 
cità (1). 

Si pud notare osservando il forte 


effetto di spegnimento mostrato dalle Fig. 22. — Impulso di particella « nella 
curve in fig. 13 che 7, & ridotto pit zona Geiger. 
fortemente che T,. 
Questo è in accordo col fatto che 
la nettezza di separazione dei fotoimpulsi migliora aumentando la percen- 


tuale di CO,. 


Effetto della carica spaziale. — La forma degli impulsi all’inizio della 
corona alle varie concentrazioni e la formazione della zona Geiger descritta 
in 1, C) si spiegano facilmente considerando Veffetto della. carica positiva 
generata sul file. Consideriamo, per semplicita, alla soglia di divergenza, l’effetto 
di riduzione del campo elettrico che & prodotto da un impulso iniziale generato 
da un arrivo uniforme su tutto il filo di q elettroni pr imari per cm. Nascono 
nella valanga di Townsend Go: N,; ioni per em. Questi ioni si allontanano dal 
filo con una velocità Wt = kE essendo k la mobilita ed # il campo elettrico. 

Luibbasspenss di potenziale prodotto da una guaina di ioni di. densità 
qo: Ns collocati a distanza 7* dal filo & rappresentata dalla formula (3) (BE 


a? 


%x2 ” 


(6) we y,-4N:in 


(3°) T, per miscele di A + CO, è data, con le notazioni usate, da 


m 

ji T's, 

Le = Ts 
N, «ßyK 


l 


(31) B. Rossı e H. STAUB: Ionization Chambers and Counters (New York, 1949). 
(32) La formula (6) v vale quando la RC del contatore & piccola, il che € verificato 


nel nostro caso. 


79 - Il Nuovo Cimento. 
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In questa formula si considera il catodo a potenziale 0, il filo a potenziale V, 
in assenza di guaina. I] campo vicino al filo in presenza di guaina corrisponde 
al potenziale V’ del filo; a à il raggio del contatore. 

La moltiplicazione per i fotoelettroni secondari che arrivano sul filo dopo 
il tempo 7 & determinata dal valore di campo che trovano nella zona di va- 
langa. 

Assumiamo in prima approssimazione la mobilita degli ioni positivi di argon 
indipendente dal campo E e dalla percentuale di CO,. Si ha facilmente che 
la distanza r percorsa nel tempo 7 & data dalla legge r*? = yT + b? con 
4 = 2kV/In (a/b) essendo b il raggio del filo e V il potenziale agente sugli 
ioni, che è dell’ordine di Vy. 

La formula (6) diventa 


= 
= 


a” 
PE Vea) us 
À Go2 k n ıT re be 


Si puod concludere che l’effetto di carica spaziale come descritto dalla for- 
mula riduce V, e quindi la ragione della successione di fotoimpulsi tanto più 
è grande N, ed è piccolo 7. 

Un calcolo della formazione degli impulsi che tenga conto di tutta la suc- 
cessione e dello spandersi lungo il filo degli impulsi € molto complicato e non 
è stato da noi abbordato. 

Pertanto, poiché con l’aumento delle percentuali di CO, vengono ridotti 
i valori di 7 ed aumentati i valori di N,, si pud ritenere che l’effetto di ridu- 
zione della ragione per via della carica spaziale à pit forte a percentuali più 
elevate di CO,. 

In particolare a percentuali di CO, > 0,05% (T << 7 us, Ns > 2000) questo 
effetto impedisce il formarsi della corona al valore della divergenza del pro- 
cesso fotoelettrico e produce l’autospegnimento degli impulsi aventi ragione 
iniziale >1. 

Si puö descrivere come segue la forma degli impulsi di fig. 14. La molti- 
plicazione iniziale ha un valore N > N, e la successione procede divergente. 
La carica spaziale formatasi in un primo gruppo di fotoimpulsi riduce il campo 
in modo che N sia ridotto minore di N,; questo succede intorno al fotoimpulso 
di valore massimo del gruppo. Di poi N è tenuto < N, fino al completo spegni- 
mento della successione fotoelettrica. 

La carica € sparsa con questo meccanismo lungo tutto il filo del conta- 
tore. La formazione della zona Geiger è chiaramente dovuta all’effetto di ca- 
rica spaziale, in accordo con quanto è stato assunto da vari autori; da nostri 
risultati si ha perd che questo accade solo con concentrazioni di CO, supe- 
riori a 5*10-4. 

L’andamento dei bursts rilevati nelle figg. 16-17 è spiegato con la forma- 
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zione e l’allontanamento degli ioni positivi, che determina un’oscillazione del 
campo sul filo intorno al valore di soglia. 

Si vede chiaramente con un ragionamento del tipo di quello fatto sopra 
che i valori di N prossimi ad N, si hanno nei massimi e nei minimi delle oscil- 
lazioni. 

La carica spaziale non riesce piü in queste condizioni a spegnere comple- 
tamente gli impulsi fotoelettronici. 


C) Zona Geiger: contatori lenti. — Vogliamo considerare un poco pit in 
dettaglio le proprieta della zona Geiger e della zona di corona instabile. Come 
si & detto queste zone sono molto ristrette e non hanno utilita pratica ai fini 
de] conteggio. 

Notiamo pertanto che sono ottenute con valori di RC da 1 a 200 us e senza 
cireuito di spegnimento esterno. E evidente che usando elevati valori di RC 
0 eireuito di spegnimento esterno si possono ottenere ampie zone Geiger poiche 
con questi artifici si pud impedire per un certo campo di tensioni il riaccen- 
dersi della corona. 

Il comportamento dei nostri contatori & identificabile con quello dei cosi- 
detti slow counters. 


D) Altri processi. — Come si & visto nei paragrafi precedenti, la forma- 
zione della zona Geiger e della corona si possono completamente spiegare a 
mezzo del processo fotoelettronico al catodo. La netta separazione dei foto- 
impulsi e la loro correlazione fino a tempi di qualche ms indica che altri pro- 
cessi giocano un ruolo trascurabile rispetto a questo processo. 

In particolare non è importante, alle pressioni e percentuali usate, un pro- 
cesso di fotoionizzazione nel gas. 

Non si è avuta inoltre evidenza del processo di estrazione ione-catodo, 
invocato da aleuni autori per spiegare la mancanza di estese zone Geiger nel 
contatori ad argon. Questo processo dovrebbe produrre elettroni dopo un 
tempo dell’ordine di 2-8 ms dalla valanga iniziale, corrispondente al tempo 
di transito degli ioni nel contatore nelle varie condizioni. 

Per quelle miscele (5:10-*< X < 5:10”) dove si hanno la zona Geiger e 
la zona di corona instabile si possono trarre alcune conclusioni più precise 
riguardo questo processo. 

a) La fine della zona Geiger è come si & visto determinata dal mecca- 
nismo fotoelettronico. 


b) Gli impulsi fotoelettronici nei bursts di corona instabile sono ben 
correlati fino a tempi corrispondenti al tempo di transito degli ioni 7+; non 
c’è arrivo importante di elettroni al tempo 7°. 


c) I massimi nei bursts di corona instabile capitano a tempi inferiori 
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a T* cioe non si possono attribuire questi massimi ad un arrivo di elettroni 
prodotti al catodo dagli ioni. 


d) Si osservano nella zona Geiger impulsi autospegnentisi ben separati 
e non correlati fra loro, ciascuno dei quali dura qualche centinaio di us. 
La carica totale contenuta in uno di questi impulsi é dell’ordine di 108 
cariche elettroniche. Si deduce da cid che la probabilità del processo ione- 
catodo & inferiore a 107$. 


E) Particelle «. - Il fatto che nella zona Geiger e di corona instabile gli 
impulsi dati da particelle « abbiano aspetto convergente € facilmente spiegato 
invocando la riduzione di moltiplicazione che si ha ad elevate moltiplicazioni, 
per la densa ionizzazione primaria di queste particelle. 

Questo effetto è dovuto alla riduzione del campo sul filo ad opera degli 
ioni positivi formatisi nella valanga che riducono lo sviluppo della valanga 
stessa. 

Si ha allora N, < N, con N, coefficiente di moltiplicazione iniziale per par- 
ticelle x ed N, coefficiente di moltiplicazione iniziale per impulsi y. Nella breve 
zona Geiger si puö assumere N,> NV, ed N,< N,; mentre gli impulsi y hanno 
inizio divergente le particelle « danno impulsi convergenti. 


8. — Miscele A + CO,. Contatori rapidi. 


Non abbiamo condotto uno studio sistematico del comportamento del con- 
tatore con percentuali di CO, maggiori di 0,5%. Abbiamo soltanto tratto aleune 
osservazioni che riteniamo utile descrivere. 

Con pressioni di 500 e 1000 mm Hg e percentuale di 1% di CO, si nota 
ad elevati valori delle moltiplicazioni l’apparire degli impulsi fotoelettronici. 

Il tempo fra i fotoimpulsi € dell’ordine di 1,5-2 us e al valore di molti- 
plicazione N, 20000-40000 si ha la soglia di divergenza del processo che da 
anche la soglia della zona Geiger. Questa zona si estende per qualche centinaio 
di volt. Gli impulsi Geiger hanno nel tratto iniziale della zona l’aspetto di 
fig. 14 e con l’aumento della tensione aumentano di ampiezza. 

Aumentando la tensione appaiono inoltre meno netti i fotoimpulsi sino a 
che l’impulso Geiger è del tipo di fig. 23 e quindi di fig. 24. 

A pressioni di 150 mm Hg la zona Geiger comincia a valori di N, 80 000 
e gli impulsi gia in soglia hanno un aspetto del tipo di fig. 23. 

Con percentuali di CO, del 10% gli impulsi Geiger hanno subito in soglia 
Vaspetto continuo e rapido illustrato dalla fig. 24 a tutte le pressioni usate e 
non c’e traccia di impulsi fotoelettronici. 
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Per questo motivo è difficile misurare à questa percentuale il valore della 
moltiplicazione di soglia. 

Le zone Geiger e le forme di impulsi ora descritte sono del tipo di quelle 
ottenute comunemente con miscele di argon e vapori organici. 

L’andamento rapido e continuo dell’impulso è indice che nello spandersi 


Fig. 23. - Impulso Geiger. A+10~ CO,. Fig. 24. — Impulsi Geiger sovrapposti alla 
Pressione = 500 mm Hg; Segnale di ta- fine della zona Geiger. A + 10-? CO,. 


ratura = 1 us; RO =1 us. Lo span- 

dersi della scarica & dovuto parzialmente 

all’effetto fotocatodico e parzialmente ad 
una produzione di elettroni nel gas. 


x 


Pressione = 500 mm Hg; Segnali di ta- 
ratura = 1 us; RO =1 us. 


della carica & interessato in modo prevalente un processo di liberazione di 
elettroni nel gas. 

Questo risultato & in accordo con la conoscenza attuale che si ha intorno 
al meccanismo dello spandersi della scarica nei contatori rapidi (*). 

Più precisamente si pud dire che a percentuali di CO, dell’ordine dell’1% 
lo spandersi della scarica avviene in parte per il processo fotocatodico e in parte 
per il processo nel gas. A percentuali maggiori invece il processo nel gas € pre- 
valente. 

La rapiditä del processo di catodo e del processo nel gas e gli elevati valori 
di moltiplicazione fanno si che l’effetto di spegnimento degli impulsi ad opera 
della carica spaziale sia più efficiente e permetta la formazione di zone Geiger 
abbastanza lunghe. 

(33) Notiamo che processi di questo tipo possono anche essere provocati dagli atomi 
molto eceitati di argon i quali possono liberare elettroni in urti contro atomi di argon, 
per esempio per mezzo della reazione A* + A> (A,)* > A,’ + e studiata recente- 
mente da HornBECK e MOLNAR (Phys. Rev., 84, 615 e 621 (1951)), oppure per mezzo 
di urti ionizzati contro atomi o molecole di un altro gas. 

Lo spandersi della scarica in queste condizioni pud essere favorito dalla diffusione 
di questi atomi eccitati in una ristretta zona intorno al filo per mezzo dei fotoni duri 
da essi emessi e che si possono considerare fotoni di risonanza. 
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9. — Miscele A + He. 


Descriviamo in questo paragrafo alcuni risultati che abbiamo ottenuto 
contaminando con vapore di mercurio l’argon contenuto nel contatore. La 
pressione parziale del mercurio corrisponde alla sua tensione di vapore ed &, 
alla temperatura ambiente, dell’ordine di 10° mm He, 


10° - = — 


igen 


500 1000 1500 2.000 2500 Volt V 
Fig. 25. — Curve di moltiplicazione per argon contaminato con vapori di mereurio 


alla tensione di vapore a 20°C; [j pressione argon— 200 mm Hg 
x pressione argon = 600 mm Hs; © pressione argon— 900 mm He. 


L’introduzione del mercurio é stata fatta in un impianto analogo a quello 
di fig. 1 comprendente due contatori di caratteristiche uguali a quello di fig. 2. 
La fig. 25 rappresenta le curve di 
moltiplicazione relative a tre pressioni. 
Alla tensione V, indicata nella figura 
si accende bruscamente la corona sta- 
bile. Come si vede i corrispondenti va- 
lori di N, sono di molto inferiori à 
quelli per argon non contaminato. 
A pressioni di 200-500 mm Hg non 
vi € traccia di fotoimpulsi nella forma 
degli impulsi rilevati in soglia di co- 
Fig. 26. — Impulsi di particelle x in ar- rona e osservati con valori di RC 
gon + vapore di mereurio alla soglia di di 1us; a pressioni di 1000 mm Hg 
corona. Pressione = 1000 mm Hg; Se- Au ER x 
gnali di taratura — 1 us; RO=1 us, la forma degli impulsi è rappresen- 
La lenta discesa è dovuta a produzione tata nella fig. 26; la gobba della discesa 
di elettroni nel gas. si riconosce essere il contributo dell’ef- 
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fetto fotocatodico; la discesa continua ha andamento molto più lento di quel 
che comporta la costante di tempo RC usata. 

Abbiamo studiato anche il comportamento di miscele A + 0,1% di CO, 
in presenza del vapore di mercurio ottenendo in questo caso lo stesso compor- 
tamento che si osserva in assenza del vapore di mercurio. 

Questi risultati si interpretano facilmente nello schema discusso nei para- 
erafi 5 e 6 e tenendo conto dell’effetto Penning: ionizzazione per urto del Hg 
da parte di atomi eccitati di argon ("). (Ricordiamo che Venergia di ionizza- 
zione del Hg (10,4 eV) & inferiore all’energia di eecitazione degli atomi eccitati 
di argon negli stati risonanti e metastabili). 

Gli atomi eccitati di À permangono un certo tempo nel gas diffondendo 
intorno al filo oltre che per usuale diffusione gassosa per mezzo dei fotoni di 
risonanza. Questi atomi hanno la possibilità di ionizzare per urto il Hg e l’elet- 
trone prodotto serve a riaccendere la scarica. 

Il nostro risultato in A +4 Hg indica che per l’accensione della corona alle 
basse pressioni questo processo è più importante di quello fotocatodico. Il 
contributo dell’effetto Penning diventa invece dello stesso ordine di quello 
fotocatodico a pressione di 1000 mm Hg, dove anche la concentrazione di 
mercurio & più piccola. é 

La forma degli impulsi in soglia di corona è in accordo con questa descri- 
zione del meccanismo di accensione della corona in quanto l’andamento con- 
tinuo della discesa è indice di un processo di produzione di elettroni nel gas. 

Riguardo al comportamento della miscela A-+ 0,1% CO, + Hg si pud 
osservare che l’effetto di spegnimento della CO, riduce fortemente la vita media 
degli atomi eccitati di argon in modo da rendere trascurabile il numero degli 
urti che tali atomi fanno contro gli atomi di Hg. 

À interessante notare Velevata criticità del comportamento dei contatori 
con argon puro a percentuali piccolissime di vapori 0 gas come il mercurio che 
sono ionizzabili per urto da parte degli atomi eccitati di argon. 


10. — Conclusione. 


In quelle miscele in cui l’effetto fotoelettronico & il processo determinante 
della scarica, si pud, dall’osservazione degli impulsi, distinguere chiaramente 
la parte di essi dovuta alle valanghe di Townsend e quella dovuta al processo 
fotocatodico. Conseguentemente & possibile studiare i processi in questione 
in modo abbastanza preciso. 

Dai nostri risultati si ottiene che il processo fotocatodico è responsabile 
della formazione delle zone Geiger e di corona in miscele fino a concentrazioni 
di CO, di 5-10-? e ad ogni pressione usata. 
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La stessa accuratezza non pud essere ottenuta con quelle miscele in cui 
è presente una produzione di elettroni nel gas. 

Processi di questo tipo sono molto interessanti in quanto producono uno 
spandersi della scarica nel contatore molto rapido. 

I fotoni agenti nel processo fotocatodico sono da attribuirsi ai primi 
livelli eccitati dell’argon. Questi fotoni sono prodotti in collisioni che portano 
alla distruzione di questi livelli e non sono fotoni di risonanza. L’ordine di 
erandezza trovato per la vita media di questi livelli é in accordo con i valori 
che si ottengono estrapolando alle nostre pressioni i risultati di MOLNAR. 

Misure dirette di questa quantita in tutto il campo di pressioni sarebbero 
molto interessanti per definire la natura del processo di distruzione dei livelli 
nell’argon. 


Desideriamo ringraziare il prof. G. BoLLA ed il prof. B. FERRETTI per 
l’assistenza dataci nel corso del presente lavoro, il dr. F. M. PENNING per le 
utili discussioni e i consigli, i dr. MOLNAR e HORNBECK per la comunicazione 
dei loro risultati e l’ing. E. GATTI per l’aiuto datoci nei riguardi dell’apparec- 
chiatura elettronica. Desideriamo infine ringraziare il prof. L. B. LOEB per 
Vinteressante critica e gli utili consigli datici. 


APPENDICE 


1. — Processi di decadimento dei livelli risonanti e metastabili nell’argon. 


A) Imprigionamento dei fotoni di risonanza. — La radiazione di risonanza 
prima di raggiungere il catodo subisce il processo di assorbimento e di riemis- 
sione da parte degli atomi non eccitati del gas. Un calcolo di questa diffusione 
è stato recentemente sviluppato da HOLSTEIN (1°). 

Questo autore ha messo in evidenza come l’assorbimento sia strettamente 
legato alla forma della riga di emissione: questi calcoli valgono in condizioni 
sperimentali tali che la forma della riga di emissione sia uguale a quella della 
riga di assorbimento. 

A pressione superiore a 1 mm Hg la forma della riga è determinata dagli 
urti che l’atomo eccitato fa con gli atomi non eccitati durante la sua perma- 
nenza nel gas. Questi urti hanno l’effetto di allargare la riga. Per rappresen- 
tare la forma di una riga viene assunta generalmente una legge di dispersione 
data da: 


K, dv 


‘its U NT Vere enn TTR 
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K(v) dy = intensità relativa emessa nell’intervallo di frequenza compreso fra 


y e y + dy; 
AS 92 Yo , 
je art rw a 
+m i Yo 
dove A, —lunghezza d’onda al centro della riga, 
Yo — frequenza al centro della riga, 
n — numero degli atomi per cm’, 


1 Jo = pesi Statistici del livello eceitato e del livello fondamentale, 
Yo — larghezza propria della riga, 


yp —larghezza della riga alla pressione 9; 


y, eresce linearmente con p; a seguito di cid K, è indipendente dalla pressione. 

HOLSTEIN ottiene, per una riga di forma rappresentata dalla (1’), che il 
decadimento degli atomi eceitati in un cilindro infinito, avviene prineipal- 
mente secondo una legge esponenziale data da 


’ 


INES (6) N exp E at t 
= 


dove N*(t) — numero degli atomi eccitati al tempo 1, 


No: — numero degli atomi eccitati al tempo { = 0, 


7* — 1/y, — vita media propria del livello, 


1/g —numero di volte che un fotone & in media catturato e 
riemesso prima di raggiungere la superficie del cilindro. 


Nella condizione di validità della (1') si ha per g 
g = 1,115 (xK,R)*”*, 


essendo R il raggio del cilindro. Si ha che g e indipendente dalla pressione 
poiché ne & indipendente K,. 

Ci sono vari motivi per ritenere che la legge (1’) non sia valida a pressioni 
elevate dell’ordine di 1 atm e più e non rappresenti bene la forma della riga 
alle sue ali (cioè per »—v abbastanza grandi). 

Rimandiamo a (1%) per una discussione più dettagliata limitandoci ad osser- 
dare che uno di questi motivi & dovuto al fatto che nella deduzione della for- 
mula si considera che la durata dell’urto (dimensioni atomiche/velocita rela- 
tiva ~ 10-118) sia trascurabile rispetto al tempo fra due urti; il che non à 
valido a pressioni elevate (per esempio ad 1 atm il tempo fra due urti à 
~ 6-10! 8). 

Non si conoscono con precisione i limiti di validita di queste leggi e nes- 
suna misura ci risulta essere stata fatta sull’allargamento di pressione delle 
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riche di risonanza dell’argon. Una verifica della teoria & stata ottenuta da 
HOLSTEIN e collaboratori nel caso del mercurio (%). 


B) Transizioni per wrto fra à livelli. — Gli atomi eccitati in uno dei 
livelli 2P,, *P,, °P,, 1P, possono transire da uno all’altro dei livelli contigui 
per urto con atomi non eccitati. 

Per esempio fra i livelli °P, e °P, si ha la reazione 


INR Se AN ess INE) a IN 


L’energia di reazione è di 0,075 eV; nel caso del passaggio *P, > °P, é 
fornita dall’energia cinetica termica degli atomi. 
Sia W,, il numero di transizioni per secondo °P, — °P}, 
W,, il numero di transizioni inverse per secondo ?P, — ®P,, 
E, ed E, le energie dei due livelli, 
per il principio dei bilanci dettagliati, considerando il sistema in equilibrio 
termodinamico, si ha 


——" = exp [—(B,— B,)/kT] « 


21 


=| = 


In modo analogo vanno considerati gli altri tipi di transizioni fra i livelli. 
I valori delle sezioni d’urto relative a queste transizioni non sono noti nel 
caso dell’argon. 


C) Altri processi nel gas. — Secondo alcuni autori alla distruzione dei livelli 
eccitati contribuisce anche un processo di formazione di molecola fra un atomo 
eccitato ed un atomo normale di argon. Questa molecola si forma in uno stato 
eccitato e pud attraverso i suoi livelli diseccitarsi fino alla dissociazione, non 
essendo stabili le molecole di gas nobile. 

I fotoni emessi in questo tipo di diseccitazione, non essendo fotoni di riso- 
nanza dell’atomo di argon, non sono assorbiti dal gas. 

Questo tipo di processo è stato invocato da PIKE (*) per spiegare alcune 
questioni connesse con i fenomeni di scarica nel neon. 

Successivamente MEYEROTT (?) ha mostrato come a questo processo & da 
attribuirsi la distruzione dei livelli metastabili dell’elio. 

Evidenza dell’esistenza di un processo di distruzione nel gas per gli atomi 
metastabili e stata ottenuta nel caso dell’argon da MOLNAR e PHELPS (?%?2), 
MOLNAR (°!) osserva una durata di questi livelli decrescente con il crescere 
della pressione secondo la legge 


(2') A= 40p + 10p?. 


A probabilità di distruzione per secondo, p pressione in mm Hg. 

Il termine p? che diventa preponderante alle pressioni di qualche em Hg 
puo essere attribuito ad una sorta di processo di decadimento, probabilmente 
formazione di molecole eccitate, collegato con urti tripli nel gas. 


(34) D. ALPERT, A. O. McCougrey e T. HozsreiN: Phys. Rev., 76, 1257 (1949). 
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D) Distruzione alle pareti. — Alle pressioni dell’ordine dei mm Hg gli atomi 
eccitati possono diffondere verso le pareti del tubo contenitore e quivi disecci- 


tarsi per urto 


Questo processo produce un decadimento dei livelli in prima approssima- 
zione esponenziale e di vita media A?/D essendo A la lunghezza di diffusione 
del recipiente (%) e D la costante di diffusione degli atomi, che & proporzio- 


nale a 1/p. 


2. — Calcolo del decadimento. 


Consideriamo per semplieitä solo i livelli *P, e ®P, e scriviamo le equa- 
zioni del decadimento. 


| 


(31) | 


essendo 
N,(t) 


N u(t) 
aN, aN, 


dt’ dt 
(g/t*)N p 
WN e 


Wah » 


iD DD 
Aan a x 


AN» ’ A aly M 


AN, Gar ee ERBEN ® 
en = mn Nr — WıNr + WaNy 13 Na —A,N», 
ay y 


sh bu 2 
dt == WN an + W Nr = PE I ap AN y 5 


il numero di atomi nello stato *P,, risonante 


il numero di atomi nello stato *P, metastabile 

le corrispondenti variazioni per unita di tempo, 

rappresenta il decadimento con imprigionamento del livello 
SP, come descritto in 1, A), 


numero di atomi eceitati nello stato *P, che per secondo pas- 
sano allo stato *P, per urto contro atomi non eccitati di argon, 


numero di atomi *P, che per secondo passano a °P, per urto 
contro atomi non eccitati di argon, 


rappresentano la diffusione alle pareti e la successiva distru- 
zione per urto, 


rappresentano un processo di distruzione nel gas essendo A, 
ed A, le corrispondenti probabilità per secondo. 


Le equazioni (3’) si risolvono facilmente e danno 


(4°) 


| Np = Nz, [ % exp [ét] +.(1 — %o) EXP [ot] | ; 
| My = Vu[Bo exp [ét] + (1 — fr) exp [Et] ; 


(35) Nel caso di recipiente cilindrico 1/A?= (2,4/R)2 + (a/L)? con R raggio e L lun- 
ghezza del cilindro. 
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con N,,e N,,, numeri iniziali di atomi eccitati nei livelli risonante e meta- 
stabile e &, e & dati da 


7) 


n +m ' d(mn — WW. 
(5) Ér2—= D A 1 ae /1 an 


essendo 


m=4r7T + W2 Ae DA = Wa Ay DA 
% € Po Sono legati dalle equazioni 


(6) Nr, = Wa 5 1 &% Nn, “ei Wo 
BoN y em . 


Le equazioni (4') danno per il decadimento dei due livelli una legge a due 
esponenziali di cui quella con decremento &, rappresenta l’andamento asin- 
totico. 

Si pud vedere come uno studio accurato dell’andamento delle popolazioni 
dei livelli permette di ottenere informazioni sui diversi processi di decadi- 
mento. 

I risultati sperimentali sul decadimento dei livelli dell’argon sono troppo 
scarsi perchè si possa da essi trarre con precisione il peso dei vari processi 
interessati. 

Si pud pure vedere con facilità che per ottenere una durata dei li- 
velli decrescente con la pressione con legge quadratica, come osservata da 
MOLNAR, bisogna necessariamente fare intervenire almeno uno dei termini 
Ar, Ay che tengono conto in modo opportuno di una distruzione per urto 
nel gas. 


SUMMARY 


The multiplication curves and the pulse shapes have been studied in a counter 
filled with pure argon and argon plus carbon dioxide mixtures at pressures ranging from 
150 to 1000 mm Hg. The pulse shapes in argon in the high proportional zone and 
up to voltages corresponding to the corona threshold are explained by assuming a pro- 
duction in the Townsend avalanche of photons hard enough to extract electrons from 
the cathode. These photoelectrons are found to be produced during a time T,, of the 
order of a few us when at 150 mm Hg pressure and lower when at higher pressures. 
The analysis of the features of excited levels of argon suggests the four lower excited 
levels to be responsible of the photoelectric effect. The T,, values deduced from our 
measurements are in accordance with those obtained by extrapolating the Molnar values 
of the first metastable levels at the pressures used. This photoelectric process has 
been found to be the main responsible of the sustaining of corona discharge, which 
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starts at multiplication values of about 200. If carbon dioxide concentrations varying 
from 5-10-4 to 5-10-3are added to argon, an increase of N,, the threshold multiplication, 
is obtained, which is linear versus the carbon dioxide percentage. This fact is readily 
explained by assuming the quenching of argon excited levels by second kind collision 
on carbon diexide molecules. From the experimental results a value of the quenching 
cross section may be deduced which is about 600-1077 cm? in a first approximation. 
Further studies concern the space charge effect, at the photoelectric process divergence 
threshold, which allows the formation of a very short Geiger zone, at carbon dioxide 
concentrations over 5-10-4. No processes other than the photoelectric one have been 
found to be efficient in counters filled with argon plus carbon dioxide mixtures up to 
concentration of 5:10-%. At higher percentages the photoelectric process is conside- 
rably reduced; there it has been found to be effective in the discharge building-up a 
process by which electrons are released in the gas body. The counter behaves as a 
regular fast counter. The effect of mercury vapor in the discharge in pure argon is 
also described, which is accounted for by the well known Penning process of mercury 
ionization by collision on excited argon atoms. This process causes the discharge to 
start at extremely low multiplication values. 
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On the Algorithm of Dirac Spurs. 


E. R. CAIANIELLO (*) 


Istituto di Fisica dell’ Università - Torino 


S. FUBINI (**) 


Istituto di Fisica del Politeenico - Torino 


(ricevuto 111 Ottobre 1952) 


Summary. — The spur of the product of any number of Dirac matrices 
is shown to be square root of a determinant. It can be evaluated directly, 
however, by means of rules as simple as those of determinant theory. 
When the number of matrices occurring in the product is >12, further 
considerations show that many terms of the expansion sum up to give 
zero total contribution. A final formula is given which takes all such 
reductions into account. The Appendix contains rules for summing 
over pairs of y which bracket products of Dirac matrices. 


In most calculations of processes involving particles obeying the Dirac 
theory, e.g. of transition probabilities in decay processes or of cross sections 
in electrodynamics or meson theories, it is necessary to evaluate spurs of pro- 
ducts of Dirac matrices. Although this is quite straightforward in principle, 
when the expressions become a little involved a considerable amount of com- 
putational work may be required, and with it arises the need for a careful 
check of the results obtained, to insure that no trivial mistakes were incurred 
into. It is also rather hard to visualize a priori, with the usual computational 
technique, what kinds of symmetries may exist among the different terms 
arising from a higher order calculation, say in electrodynamics. 

In an attempt to improve on this situation, we show here that the eva- 
luation of the spur of the product of any number of Dirac matrices can be 


(*) Now at the Institute for Theoretical Physics (CERN) University of Copenhagen 
Denmark (on leave from the University of Rome, Italy). 
(**) Now at the Istituto di Fisica dell Universita - Torino. 


ON THE ALGORITHM OF DIRAC SPURS 1219 


systematized into an algorithm closely related to, and as simple as, that of 
determinant theory. More precisely, formula (12) shows that the square of 
any such spur is a skew-symmetric determinant of even order; the spur itself, 
however, can be derived from a scheme alike to that of determinant theory 
(sect. 3). 

The expanded spur is very conveniently re-expressed by (17). This for- 
mula permits to realize that remarkable simplifications arise when the order 2n 
of the spur (i.e. the number of its factors) is > 12, which become the more 
effective the larger is n. 

In actual problems, e.g. in electrodynamies, one has to sum over pairs 
of y“ which bracket products of Dirac matrices (1). Although the present 
technique makes it rather easy to sum during and after the expansion of the 
spur, it is desirable to te able to perform the summation before. This question 
is considered in the Appendix. 

The present work will be followed by another in which this technique is 
applied to the problem of summing into one compact expression all the con- 
tributions to a given process from the n' order graphs of the Feynman-Dyson 
theory. 


1. — Let 
(1) pb al Ty = Pi; (u = 1,2,5, 4) 


(à) 


. L is (i Dre DE 
be a vector matrix formed from the scalar p, , the vector pi ; the unit matrix 


and the Dirac matrices (at = ict): 


y | 


yyy” + VV = Lu : 


We are interested in evaluating the spur of a product PP, Pr. TO reduce 
P,P,... PA to a form more suitable to our purposes, we rewrite it as (?) 

[ 5 Port y* ù WE or 
(2) Phase ysePse Piper PE oe 4 ; 

{ 5 Pr 5% 5 ye 


on introducing the dual representation of the y-matrices 


(3) IT; == y? es yyy ys 2 RE a dy ye : (PY PIE a II" = 20, ; 


() R. P. Feynman: Phys. Rev., 76, 769 (1949). 
@) G. Prrrau: Journ. de Phys., 12, 565 (1951). 
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it is seen from (2) that any such product can be written as à products of an 
2 n aS Fix (A) Le 
even number of matrices Q, = Iq)" (u =1, 2, ..., 5) 


Q:Q: CCE en ’ 


where ig! =p) (@ =1,..,4); q =(-1)p. Introduce now the abbre- 
viated notation 


(4) (123... 2n)= %4 Sp (O9 --- Om) « 
From 
9,9: + OQ, = 29,4, = 2(rs), (u =1,..., 5) 


one has 
(5) (123...2n) = > (—1)*(k)(2 3..., k—1, k +1, ..., 2n), 
and, finally, by recursion: 


(6) (123...2n) = > (— 1)? (tyt2)(tst4) ee (don-aton) 5 


where >” is defined by the prescription that the sum must be taken only on 
the permutations of indices satisfying the limitations 


— 
I 
— 


Ua ce tena < ten | ty << tg << ana 3 


p is an integer, even or odd with the parity of the permutation 7,2,2, ... don. 
In the following we say, for short, that the parity of a permutation of indices 
is p, Q,...; relations among such parities are always understood to hold to 
within an arbitrary additive even integer. 

Formula (6), which is the trivial extension to our five dimensional formal- 
ism of a similar expression given by DUFFIN (*), is the starting point for the 
considerations of sect. 3. The whole matter is rendered more intuitive, how- 
ever, by a few algebraic remarks which form the object of the next section. 


2. — Intreduce the 2n alternating variables &, &,..., &,, i.e. variables 
obeying the commutation relations 


(8) seers fs) WU 


(3) R. J. Durrin: Phys. Rev., 77, 683 (1950). 
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Define the algebraic form 


(9) : CO > (hk) En, Ee = > AnrEnËr - 


h<k 


In (9), and throughout this section, the (hk) = (kh) can be any numbers; the 
an are defined by the relations 


MEI for bk 
| 


(10) Apr = (hk)nıx » nk = | 0 TOL % 


Then, 


1 
— (es Anxénér)" - 


(11) or = n!{12...2n)ÉË... Fon = 5 (2, 
((11), when the last term is expanded, yields immediately formula (14) of 
sect. 3). Write now 


Sal 
w= MI améné » 
where all the variables &, &' anticommute as in (8). Then: 


1 


fel Ay EN eS fa fags = 
5152". 5 = Gre ( One UrsEnSusrSe) == 


— (n!)} Aëé vee Sees ... a , 


with A = det. (ax), which shows that 
(12) HR) A 


The fact that a skew-symmetric determinant of even order is the square 
of the form (1 2,...,2”) defined by (5), (6), is quite well known; such form 
are named Pfaffians and their theory is about one hundred years old. In it, 
as we just said, the numbers (hk) are quite arbitrary. (In our case, however, 
they are strongly restricted by the requirement (hk) = dar, whose effect- 
iveness will appear in the next section). 

It goes without saying that (12) supplies an additional method for eval- 
uating Dirac spurs; it states, indeed, their fundamental algebraic property. 
We do not insist on this point here and just confine ourselves to pointing out 


80 - Il Nuovo Cimento. 
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a rule which permits to derive all the terms of (6) frem a scheme as simple 
as that of determinant theory. Form the array 


DD ES CE EAN ee | (1, 2n) 
(02) CONTRER (2, 2n) 
(SDR (3, 2n) 


This symbol can be regarded as representing the Pfaffian (6) if the following 
algorithm is connected to it: 

1) each term (hk) belongs to two lines: line h, containing all terms with 
an index h, and line k, containing all terms with an index k; 

2) crossing out lines h and k reduces (13) to a new Pfaffian: the minor 
of (hk); 

3) (18) is expanded by adding the products of all terms of a fixed line 
by their minors, each product taken with the sign (—1)***, 

Further properties of Pfaffians will be found in the mathematical litera- 

ture (‘). Reference to this algorithm will help in visualizing the results of 
the next section. 


3. — Formula (6) contains >’, defined by the prescription (7). This limi- 
tation can be removed by introducing the quantities a,, defined by the re- 
lations (10). One readily sees then that: 


1 
(14) (1, 2... 20) = —, D (1) gs, ir 


(2n) !! 34 2n-1 2n 
where we remind that (2n)!! = 2n(2n—2)...2=2'! and (2n—1)!! = 
= (2n)!/2”n!. Instead of (7), one can evidently impose other prescriptions 
of the same sort. We use this freedom to give a privileged role to 2x out of 
the 2n vectors q. Let 


Kk<k<..<h,; M<Mm<..<Ns; («+ Bp =n) 


be a permutation R of the indices, with parity o, and let >” denote sum- 
mation over only those permutations satisfying the limitations 


ty ns tg ao lan Se da ae > Von 


r 


(?) See e.g. T. Muir: Contribution to the history of determinants, 5 vols. (the last 
published London, 1930). 
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where i —>j means that à precedes j in R. (14) can then be written as: 


15 Don) S(O NS? (ae ne: : 
da) ( 2) Ge > a) Wizig Vigig Vinny ton 

where p’ means parity of i,i,... 4,, relative to R (ie. p = 0 + p': parity com- 
puted by taking R as giving the fundamental ordering). Consider now the 
permutation S (r +i—a; t+B—n): 


D hei ..< hes: Wig <li N 


where k’, k’ are a permutation of the indices k, and h’, h’ of the indices h. 
Let o be the parity of S relative to R. « being a fixed number, >" splits into 
& +1 partial sums, numbered r= 0,1, 2,..., %, each containing all the pro- 
ducts in which 2r of the q, combine among themselves, 2t of the g, combine 
with 2t of the q, and 2s of the q, combine among themselves. We can accord- 
ingly re-write (15) as 


2r 


(16) Bin) = (IS SEL) > ce 1) an, a, aol 


CA A D AOC ON Nr On ROC Gn, oh; 
a *2t 10072 2s-1 '2 


Here, p’=o + 7"5 ku. k,, is a permutation Ot bone NL PORC a dee 
Cy 


notes summation over all combinations ©, which can be formed for given pr. 

| = : AR arıtır BY af AE 1Q arity 

But p'— À +u +y, where À is the parity of ky y ki, Ki fe 38 the parity 
"an "an (2 H . + / 2 / NT DE 

of ky hi kiki, … k,,h., and » is the parity of h;, hy hy. (relative to the 


order in which the same indices appear in R). Also, from (15) or (16) 


/ . 
DEN NO ee — (EE, … hs); 
2 2r 1 2r 
> (— 1 "an; hi, see Oy’. hi = (hth, vee h},) , 
1 °2 2s-1 2s 


and, denoting with yw’ the parity of Wis N hy. alone (202—t~wt=a—r 
is the number of trauspositions required to separate the k” from the A): 


a Ian: + Ai = > (1) *# Urin; + Ugh; = 
qty Me 2t 


2t 
ttt m, 
sh onen. hs N x 
nf ae x! 
EN ee 
iz ( 1) £ | ( ) 1 ph" MW a 
E LÉMLAORO ose jf 
[4 / 
Weg thy DOC Weg sho | 
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the last being merely à convenient abbreviation to denote the determinant. 
We obtain thus 


tai „ 
œ Ka altes 


20) AED DICHTER re) (Who. Wi 
=0 6 


28) I 


Sat hee) 
and finally: 
: j & kik; Ass k,, = 
(17),- L22n SEND G1) ie) ne 
Be) C, h, h, ons Voy 


where 7 is the parity of the permutation 
Wi he ce The ace hop he Rey ie Wes a ee Là 


aken again with respect to the fundamental ordering 1, 2,..., 2n. 

The formula just obtained rests only on the commutation properties of 
the Dirac matrices: it is therefore valid regardless of the number of dimen- 
sions d of the space spanned by the vectors qg, provided these result from a 
theory in which the same commutation relations hold. (One migth envisage, 
for instance, a generalized Dirac theory with 8-component spinors functions 
of six coordinates, the two additional ones being charge and mass. Or one 
might wish to apply (17) to non-relativistic calculations, with 2-component 
spinors — in this case, however, it is not possible to define a 5-component 
vector q with the mass term as fifth component). 

(17) is completely equivalent to (6); it presents, however, remarkable ad- 
vantages over the latter. Consider, for instance, a standard electrodynamic 
calculation. Simplifications arising from energy-momentum conservation can 
be introduced in (17) before performing its explicit expansion. Furthermore, 
one may take as g. the electron momentum-mass vectors and as q, the vectors 
with all but one components equal to zero arising from the y" of the inter- 
actions with the photon field: this would make a subsequent summation over 
photon polarization quite easy. 

When n>6 (— d +1: we consider only the case d = 5) the formalism 
introduced here permits simplifications not attainable with the customary 
techniques, the greater the larger is n. These result if we take as permutation S 
the fundamental permutation 1,2,..., 2n. It is convenient now to denote 
with k,, ks,..., k,, a permutation of the numbers from 1 to 2x such that 
hy < hig <<... < hors korrı<..<k,, and with h,,,,,..-, Men à permutation of the 
remaining indices such that h,,.ı<.. <a) Myy or 4 +e < Non. The de: 
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terminants appearing in (17) become with this choice simply 


ls.) 00e pr ser; 
Kor+ Korte ee soie k, | ; : = 1) - . con a) 
RAR + oh, 


2x+1 DALDZ so Ye 


(ER Oils od on ¢ (Bees We 


But in a space with d dimensions there are at most d linearly independent 
vectors: it follows that, in our case, all the determinants appearing in (17) of 
order six and greater vanish identically. Only the terms with r=—«, a—4, 
x—2 survive. When n is large, formula (17) must be applied reiterately; 
the choice of « which secures the maximum reduction is & = 3, at all stages 
(it can be shown that the best choice, for d-dimensional gq, is (d + 3)/2 > 
+a>(d +1)/2. The proof of this statement is left to the reader). The final 
result is therefore: 


(18) 2.2) = 12345 6) (7 89... 2n)— 


This [ies 
= »S (—1)"(% ko Ks hs) (Ii Ray --- Mon) 4 
C2 h he 


Pa 


‘Tis ka Ks k 
+3 D |, ; ) ne hen), 


© hz hg No hao 


We are grateful to Prof. G. WATAGHIN for his kind interest in our work. 


APPENDIX 


Let 
RO Oreo es (Rn) = 4 Sp(QiQz ++ Qn) - 


On introducing the 16 matrices of the Dirac algebra (assumed here to be 
Hermitean): 


DUT Pe=itywaryy Fe Sys Tia ays r=} 


one has 
(Al) R, = > (RTS ; 
and also 


(A2) À, = OO SI URL 
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It follows from (A1) that 


(A3) (RR) = I (RITTER). 

This formula permits to reduce the calculation of a spur of given order to 
lower order problems; it can be shown to lead again to (17). From (Al) one 
finds immediately: 


(A4) x yeRyy! = 4(R,) — 2 D (Ray?)y? + 2 D (Rröyey)iyey® — 4(Rry*)y* , 
Q e 
etc. 

Relevant simplifications occur when all qf? = 0. It is just as well then 
to consider the matrices (1) directly: we are in the case usually considered. 
Call, to avoid confusion, X, = yep, (u = 1, 2, 3, 4) such a matrix, and let 
Y, = X,X,... X,. The following identities follow immediately from those re- 
ported above: 


(A5) ee 2 any + (Yony?)Y ; 

(A6) VE — 2 (Yorııy")y" + 2 (YorıiyeyP)iyey , 

(A7) Ir Yayt = 4( Yon) — 4(Yony®)y® 

(A8) dy" Yoyıy! =— 2 > (Yortıye)y® + 2 > (YarniyeyP)iyey?) = — Cha 5 


(A3) results simplified accordingly: in applying it we can take either 


Prior Day (Fan) =) 


RIASSUNTO 


Si dimostra che la traceia del prodotto di un numero qualsiasi di matriei di Dirac 
è uguale alla radice quadrata di un determinante. Essa perd pud essere calcolata diret- 
tamente per mezzo di regole, tanto semplici quanto quelle della teoria dei determi- 
nanti. Se il numero delle matrici che formano il prodotto & >12, molti termini dello 
sviluppo si elidono. Si ottiene una formula finale che tiene eonto di tutte queste ridu- 
zioni. L’appendice contiene delle regole per sommare su coppie di y“ fra le quali si 
trovano dei prodotti di matrici di Dirac. 


NODH TECNICHE 


Sul funzionamento di un particolare circuito a coincidenze. 


F. DEMICHELIS e R. MALVANO 


Istituto di Fisica del Politeenico - Torino 


(ricevuto il 12 Ottobre 1952) 


Riassunto. — Si discute sulla variazione del tempo di risoluzione di un 
particolare cireuito a coincidenze in funzione della energia dei raggi y 
incidenti su contatori a scintillazione. Si segnala come il dispositivo per- 
metta la misurazione dell’energia media dei raggi y. 


In aleune ricerche sui fenomeni di correlazione angolare y-Yy @)esine 
fatto uso di un dispositivo a coincidenze il cui schema & riportato in fig. 1. 
Nella misurazione sperimentale del tempo di risoluzione t di questo disposi- 
tivo a coincidenze si è osservato il comportamento che desideriamo qui 
segnalare. 

T due impulsi di tensione generati nei fotomoltiplicatori F, ed F, (931A) 
da una scintillazione prodottasi rispettivamente nei cristalli C, e C, (Nal[Tl]) 
vanno anzitutto a due amplificatori A, ed A, (fattore di amplificazione ~ 40) 
e da questi giungono alle eriglie di due univibratori asimmetrici costituiti da 
due EFF51. Questi agiscono da « formatori » degli impulsi destinati ad operare, 
mediante un sistema di derivazione capacitiva, sulle due griglie di un doppio 
triodo 6SN7, a uscita catodica, che è destinato a segnalare l’eventuale coinci- 
denza di una coppia di impulsi. Il numeratore N registra il numero di coinci- 
denze rivelate dal dispositivo; I, ed I, registrano il numero di impulsi ad unità 
di tempo, provenienti dai due formatori. 

Il dispositivo a coincidenze trae il suo valore dalla conoscenza del suo 
tempo di risoluzione r, definito come la metä dell’intervallo di tempo entro 
cui devono arrivare i due impulsi perche Vapparecchio li segnali come coin- 
cidenti. 

I metodi consueti per la determinazione sperimentale di 7 in un circuito 
a coincidenze sono i seguenti: 

1) utilizzazione di impulsi artificiali provenienti da un unico generatore 
elettronico di impulsi ed inviati alle due entrate del circuito a coincidenze 
con ritardo di tempo noto uno rispetto all’altro; 


() F. Demicueris e R. MarvaNo: Nuovo Cimento, 9, 1106 (1952). 
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2) utilizzazione di due sorgenti radioattive S, ed S, costituite da uno 
stesso nuclide e poste come in fig. 1, disponendo un opportuno schermo X 
in modo che ciascuno dei due rivelatori riceva solo le radiazioni provenienti 
dalla sorgente più vicina. Il tempo di risoluzione 7 risulta allora: 


en 
NN 


dove N, indica il numero di coincidenze ad unità di tempo registrate al nume- 


ratore N, ed N, e N, sono i numeri di impulsi ad unità di tempo segnalati 
rispettivamente da I, e I, (2). 

Per quanto il metodo 1) sia più rapido e più semplice, il metodo 2) ha il 
vantaggio di lavorare in condizioni il più possibilmente analoghe a quelle della 
ricerca sperimentale che dovra far uso del dispositivo a coincidenze. 

Per questa ragione è stato preferito nella nostra ricerca il metodo 2). 


(?) P. E. Cavanacu: Cap. 60 in O. R. Frisch: Progress in Nuclear Physies 


(London, 1950). 
(?) F. METZGER e M. Deurscn: Phys. Rev., 74, 1640 (1948). 
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Nello studio per la sua migliore messa a punto è risultato il fatto, che 
finora non ci sembra ancor segnalato, sebbene sia fonte di interessanti consi- 
derazioni che qui seguono, che: nel circuito di fig. 1 il tempo di risoluzione T, 
se perfezionando la tecnica scen- 
de a valori minort di 1077 8, 
risulta funzione lineare dell’ ener- 
gia delle radiazioni incidenti. 

Le misurazioni sono state 
effettuate ripetutamente utiliz- 
zando volta a volta un diverso 
nuclide radioattivo. Si ebbero 
cosi tre serie di misurazioni con 
coppie di sorgenti costituite ri- 
spettivamente da: 331, "Ra in 
equilibrio con i suoi sottopro- 
dotti, 37Co. 

Le misurazioni con il “Ra 
sono state eseguite in tre con- 
dizioni diverse: Fig. 2. 

a) senza assorbitore; 
b) con spessore di 9,5 em di Pb (5,65 g/em?) tra sorgente e ricevitore; 
e) con spessore di 1 em di Pb (11,3 g/em?). 

Sicché si poteva disporre in totale di cinque determinazioni. 

I nuclidi da noi utilizzati come sorgenti non emettono raggi y monocro- 
matici. Per rappresentare graficamente i valori del tempo di risoluzione T 
in funzione dell’energia E dei raggi y incidenti, abbiamo preso come variabile 
quella che potrebbe dirsi energia media dei raggi y 0, meglio, il « valore 
baricentrico » dello spettro y 
del nuclide in esame; esso à 
indicato con > n-hy/> n nella 
Tabella ed è espresso in MeV 
h | volt (I) come ascissa di fig. 2. Qui so- 

1 no riportati i 5 punti sperimen- 

tali di cui potevamo. disporre, 

e Vintervallo di energie era 

== tra 0,4e 1,3 MeV. Si nota che, 

3 entro gli errori sperimentali, 

A contenuti entro qualche per- 

j cento, i 5 punti stanno su una 

retta. 

i u Inoltre, misurando il tempo 
900 1000 volt ¢ c an Éd 

di risoluzione con diverse ten- 

sioni ai fotomoltiplicatori e uti- 

lizzando sempre la stessa sor- 

gente, $Co, si ottengono i dati di fig. 3, nella quale sono riportati i valori del 

tempo: di risoluzione t in funzione della tensione V ai fotomoltiplicatori (cur- 

ra I) e i valori dell’altezza h degli impulsi, il cui ammontare più probabile © 


v= 840 volt 


Sy 
eg 


MeV 
se 


in 


t |x10°”s (D 


Pig. 3. 


(?) R. D. Evans: Nucleonics 1, 2, 32 (Ottobre 1947). 
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stato misurato con un oscillografo a doppio raggio, previa amplificazione li- 
neare, in funzione della tensione V ai fotomoltiplicatori (curva IT). 

Per un certo intervallo dei valori della tensione ai fotomoltiplicatori, la 
fig. 3 mostra la linearita tra il tempo di risoluzione e la tensione e tra Valtezza 
degli impulsi e la tensione stessa. 

I risultati di fig. 2 e, essendo noto che l’altezza degli impulsi € proporzio- 
nale all’energia dei raggi y incidenti, i risultati di fig. 3 sono in accordo con 
il fatto che il tempo di risoluzione 7 è, entro l’intervallo da noi studiato, fun- 
zione lineare dell’energia della radiazione incidente. 

Il comportamento del cireuito a coincidenze in funzione dell’energia dei 
raggi y incidenti & facilmente spiegabile tenendo conto dell’influenza che la 
tensione di innesco ha sul tempo di salita di un univibratore (°). Infatti tale 
tempo di salita & determinato essenzialmente dalla larghezza di banda degli 
amplificatori di cui & costituito il circuito e dal guadagno nell’anello rigene- 
rativo. Perd questo vale solo quando la tensione di innesco è appena suffi- 
ciente ad iniziare l’azione rigenerativa; quando invece la tensione di innesco 
è cosi violenta da dominare completamente il passaggio da uno stato all’altro, 
il tempo di salita puö risultare molto più breve. 

Osserviamo che il eurioso fenomeno accennato pud pensarsi un inconve- 
niente per talune utilizzazioni del circuito a coincidenze, nelle quali sarebbe 
comodo disporre di un appareechio a tempo di risoluzione costante. 

Non escludiamo che il fenomeno da noi segnalato sia sfuggito ad altri rela- 
tori, sia perchè non notato, sia perchè non fu spinto il tempo di risoluzione 
fino ai maggiori valori raggiungibili con il dispositivo utilizzato. 

Ma, accanto a questo lato che potremmo dire negativo della segnalazione 
da noi compiuta, riteniamo di dover riconoscere anche un lato positivo. 

Il dispositivo permette una rapida misura, nell’intervallo e con l’appros- 
simazione gid accennati, dell’energia media dei raggi y emessi da un nuclide 
radioattivo. 

L’energia media cosi determinata potrebbe essere utile nello studio delle 
specie radioattive aventi uno spettro y complesso, e una segnalazione del dispo- 
sitivo che facilmente porti a cid ci & sembrato, percid di un qualche interesse. 


Ringraziamo vivamente il prof. E. PERUCCA per il costante interessamento 
e l’aiuto fornitoci. 


(5) W. C. ELMORE e M. SanDs: Electronics (New York, 1949). 


SUMMARY 


We discuss on the variation of the resolving time of a particular coincidence circuit 
versus the energy of y rays, detected by scintillation counters. We point out how this 
apparatus could be used for the measurement of the average energy of y rays. 
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Kernphotospaltung am B!’ und B' dureh y-Strahlen bis 31 MeV 


R. MÜLLER und P. STOLL 


Physikalisches Institut der Eidg. Technischen Hochschule - Zürich 


(ricevuto il 17 Ottobre 1952) 


Die Photoeffekte an leichten Kernen verdienen ein besonderes Interesse, weil 
die theoretische Deutung auf Grund der Compound-Kern-Vorstellung nicht (durch- 
wegs durchgeführt werden kann. Im speziellen Fall der untersuchten 3-Teilchen 
Spaltungen am B!! und B! handelt es sich um Reaktionen, die über zwei Stufen 
verlaufen. In beiden Fällen kann gezeigt werden, daß vorwiegend zuerst ein ange- 
regtes Bes entsteht, das dann unter Alphaemission weiterzerfällt. Es besteht also 
die zusätzliche Möglichkeit, die verschiedenen Zustände im Be$ zu untersuchen. 
Das erste Teilchen, ein Triton oder ein Deuton, kann wegen der guten Diskriminierung 
der Teilchen in der photographischen Platte mit Sicherheit erkannt werden. Im 
Gegensatz zu der Kohlenstoffspaltung C!2(y, 3x)(1?) können hier die Anregungsenergien 
(P*-Werte) in Be® direkt ausgerechnet werden. Folgende Reaktionen wurden mit 
Hilfe von Kodak-NT1-Platten mit Boreinlagerungen untersucht: 


1: BU + hp = H$ + Het + Het QU I3 Mev; 
DR Blo + hy > H? + Het + Het Q = 5,89 MeV. 


Zur Auslösung der Reaktionen ist das Bremsstrahlungsspektrum eines 31 MeV 
Betatrons verwendet worden. 

Die bisher ausgeführten Untersuchungen (!) der Reaktionen 1 und 2 wurden 
mit dem Spektrum der Lithiumgammastrahlung (Z,,: 17,6; 14,8; 12,3 MeV) durch- 
geführt. Die Wirkungsquerschnitte o wurden für die verschiedener Linien mit 
denjenigen der Kohlenstoffspaltung C!2(y, 3x) (?) (2) verglichen. Über den Reaktions- 
mechanismus kann gesagt werden, daß 90% der Spaltungen (1 und 2) über einen 
angeregten Zustand von Be® verlaufen müssen: 


la: By, T)Be®®  Bes* > Het + Het, 
lb: Boy, D)Bes*  Bes* > Het + Het. 


() O. Rocuar, P. Srorr: Helv. Phys. Acta, 25, 451 (1952). 
(?) H. GLÄTILI, O. SEIPPEL e P. STOLL: Helv. Phys. Acta, 25, 491 (1952). 
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a) Entwicklung von Kodax-NTi-Platten mit Boreinlagerungen. 


Wegen der alkalischen Wirkung von Borax sind diese Platten viele empfind- 
licher als gewöhnliche Emulsionen. Die Temperaturentwicklung mit Amidol bei 
einem pH-Wert von 7,6 gestattet die Unterscheidung der Spuren der Alphateilchen 
von denen, die einfach geladen sind. Eine Verwechslung mit 


den häufig auftretenden 
Kohlenstoffsternen ist deshalb ausgeschlossen. 


b) Unterscheidung der Reaktionen am Bt! und B!®. 


Für die Trennung der möglichen Reaktionen wird das Impulskriterium verwendet. 
Für jeden Stern wird der Fehlervektor A bestimmt. 


? 


3 
Au L. 


P;: Impuls der Alpha oder einfach geladenen Teilchen (Triton, Deuton); 


P.: Impuls vom Gammaquant. 


Die A-Verteilung des gesamten Materials gibt uns ein Maß für die Auslese der 
Sterne. Andere 3-Teilchen Reaktionen können mit Sicherheit ausgeschlossen werden, 
wenn als Grenze für den zuläßigen A-Wert | A| = 0,7 angenommen wird. (Impuls- 
einheit 1 — Impuls eines Alphateilchens von 1 MeV). 

Die praktische Auswertung wird mit einer eigens dazu konstruierten Rechnungs- 
maschine, besorgt, die in einer anderen Arbeit (*) beschrieben wird. 


c) Wirkungsquersehnittsverlauf der Reaktion B!!(y, T) 2He*. 


No.5 Mev 
10 


14 16 18 20 22 24 26 28 MeV 


Fig. 1. — Energieverteilung von 113 Sternen Der Reaktion B* (y, PB), 2. He% 


Als Eichpunkte können die beiden Wirkungsquerschnitte bei E, = 17,6 und 
14,8 MeV (Lithiumgammastrahlung (*)) benützt werden: 


By, T)2 Het: — oay—az,emew) = (1,75 + 0,53)-10-% em? 


(*) Helv. Phys. Acta (in Publikation). 


1234 R. MÜLLER und P. STOLL 


Fig. 2. - Wirkungsquerschnittsverlauf B!! (y, T) 2 He 
(eingezeichneter Fehler: mittlerer statistischer Fehler). 


Das Verhältnis der Wirkungsquerschnitte bei den entsprechenden Quanten- 
energien beträgt: 


5 (17,6) 


== — 3,800; 
o (14,8) 


oe) 


Dieses Verhältnis stimmt innerhalb der Fehlergrenze gut mit den Betatron- 
messungen (Fig. 1. und 2.) überein. 

Das Betatronspektrum wurde nach SCHIFF (5) berechnet. Als weitere Korrektur 
bei höheren Quantenenergien kommt die Austrittswahrscheinlichkeit aus der empfind- 
lichen photographischen Schicht dazu. 

Der maximale Wirkungsquerschnitt von oy, ~ 2,5:10-8 cm? wird bei 
E, = 23 MeV erreicht. Auffallend ist der konstante Wirkungsquerschnittverlauf 
von EH, = 16 MeV bis E, = 24 MeV (>2-102° cm’). Über den steilen Wieder- 
anstieg zu einem resonanzähnlichen Maximum bei 27,5 MeV können wegen der 
schwachen Statistik keine weiteren Angaben gemacht werden. Andeutungen für 
eine Feinstruktur der Kurve über den Wirkungsquerschnittsverlauf sind vorhanden. 


d) Wirkungsquersehnittsverlauf der Reaktion B'° (y, D) 2 He“, 


Der Wirkungsquerschnittsverlauf weist bei 10,8 MeV ein ausgesprochenes 
Maximum (Gaz = 1.5:10-27 em?) auf. Am B1°-Kern ist aus anderen Messungen ein 


(3) L. I. SCHIFF: Phys Rev., 70, 87 (1946). 
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breites Niveau bei 10,83 MeV (4) bekannt. Das Absinken des Wirkungsquerschnittes 
bei E,> 13 MeV kann mit dem Anstieg der Konkurrenzreaktion: B'(y, np) 2 He 


No,5 Mev 
10 
5 
0 LE | art 
10 12 14 16 18 20 22 MeV 


Fig, 3. — Energieverteilung von 110 Sternen der Reaktion B!!(y,D) 2 Het. 


erklärt werden. Bei Æ, = 14,8 MeV (Lithiumgammastrahlung (1)) beträgt c11,8 = 
— 5-108 em?. Eine zusätzliche Kontrolle ergibt sich aus der Tatsache, daß 


u 


mit Hilfe des Isotopenverhältnisses B1/B1° und dem bekannten Bor-Gehalt der 


Fig. 4. — Wirkungsquerschnittsverlauf B°° (y, D) 2 He* 
(eingezeichneter Fehler: mittlerer statistischer Fehler). 


NTI-Platten aus der Anzahl der B1-Sterne bei Æ, = 18,4 MeV der Betatronplatte die 
Anzahl der B!!-Sterne der gleichen Platte bei Æ, = 17,6 MeV berechnet werden kann. 


e) Die Anregungszustände (#*-Werte) des Be®. 
Für das Histogramm werden die Sterne der Reaktion B!!(y, T) 2 He! benützt: 


E* = Errotal 7 11/8 Erpriton ’ 


Emotal: Gesamtenergie des Sternes, 
Emriton: Energie des Tritons. 


(4) F. AJZENBERG, T. LAURITSEN: Energy Levels of Light Nuclei IV (in Publikation.) 
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Die Verteilung der #*-Werte ändert sich mit der Quantenenergie; so werden 
bei #, > 22 MeV vorwiegend die Bes-Niveaus von 4,9; 6,8 MeV und höhere ange- 


Fig. 5. — E*: Verteilung’ der Reaktion B!! (y, 7) Be°* 
Be‘* — He!+ He‘ 
E*: Anregungsenergie des Be°-Kernes 
schwarz: E, >21 MeV; strichliniert: 18 << Hy < 21 MeV; weiss: Ey < 18 MeV. 


regt (schwarz markiert); bei H. > 22 MeV, entsprechend den Ergebnissen (1) mit 
der Lithiumgammastrahlung, Niveaus von 2,9 und 3,4 MeV. Einige Grundzustände 
sind beobachtet worden; bei 17,6 MeV beträgt das Verzweigungsverhältnis : 


o (Gr. Zustand) 
o (Total) 


= 16+ 


Die Anderung der Ubergangswahrscheinlichkeit fiir die verschiedenen angeregten 
Zustinde des Be$ bei entsprechenden Quantenenergien deutet offenbar darauf hin, 
daß neben der wahrscheinlichen elektrischen Quadrupol und magnetischen Dipol- 
Absorption der einfallenden Gammaquanten, bei größeren Energien auch die ele- 
ktrische Dipol-Absorption eine wesentliche Rolle zu spielen beginnt. Eine Winkel- 
verteilungsmessung der Tritonen (1. Teilchen) gegenüber den einfallenden Gamma- 
quanten würde uns zusätzliches Material über den Absorptionscharakter der Gamma- 
strahlung liefern. 
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Spettro della radiazione y emessa nella reazione 
10B,(n, «). tit, it Li +y 


A. Braccı, U. FACCHINI e A. MALVICINI 


Laboratori CISE - Milano 


(ricevuto il 23 Ottobre 1952) 


Gli spettrografi a scintillazione per- 
mettono di effettuare la misura dell’ener- 
gia di radiazioni yin casi in cui Vintensita 
di queste radiazioni sia molto debole. In 
particolare & accessibile a questa tecnica 
lo studio dei y prodotti in reazioni di 
cattura di neutroni anche utilizzando sor- 
genti di neutroni molto modeste; a scopo 
preliminare abbiamo studiato 1 y emessi 
nella reazione di cattura di neutroni lenti 
da parte del !°B, che sono gid stati stu- 
diati con altra tecnica. 

Lo spettrografo per raggi y che ab- 
biamo realizzato, utilizza un fototubo 
RCA 5819 ed un cristallo di Nal atti- 
vato con TI (22%). Il cristallo cubico, di 
1,25 em di lato, viene accuratamente pu- 
lito e quindi posto in un piccolo reci- 
piente metallico, chiuso con un disco 
di plexiglas, a tenuta d’aria allo scopo di 
isolarlo dall’umiditä atmosferica (fig. 1). 
La parete interna del recipiente è resa 
diffondente per ottenere una buona col- 
lezione di luce sul fotocatodo ed il con- 


@) R. B., HOFSTADTER € J. A. McINTYRE: 
Nucleonies, 7, 32 (1950). 

(2) V. O. ERIKSEN, S. JENSEN e A. SUNDE: 
Physica, 18, 91 (1952). 

(®) S. JoHansson: Ark. jor Fysik, 3, 533 
.(1952). 
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tatto ottico fra le superfici del cristallo, 
del plexiglas e del fototubo è assicurato 
con olio di vaselina. Gli impulsi all’uscita 
del fotomoltiplicatore vengono amplifi- 


Fig. 1. — Disposizione sperimentale: 1) eristallo 
di Nal: 2) involucro del cristallo; 3) base di 
plexiglas dell’involucro; 4) boro compresso, spes- 
sore 0,3 g em”?; 5) fotomoltiplicatore; 6) piombo; 


7) sorgente di polonio-berillio; 8) 
9) scatola a tenuta di luce. 


paraffina ; 


cati con un sistema lineare e inviati ad 
un analizzatore di ampiezza a 99 canali. 

La fig. 2 mostra lo spettro delle ra- 
diazioni y ed X emesse nel decadimento 
del 1#Au e più precisamente la riga y 
di 411 keV e la riga X di 70 keV del 


1238 
198Ho, come ottenuti dal nostro spettro- 


grafo. 
La sorgente & in questo caso costi- 


0,411 MeV 


Fig. 2. — Spettro delle radiazioni y ed X emesse 

nel decadimento del Au, La riga a 0,411 MeV 

© emessa dal nucleo eccitato di 2e, mentre 

la riga X & successiva alla conversione interna 
dei y di 0,411 MeV. 


tuita da una lamina di Au attivato con 
neutroni, posta direttamente sopra l’in- 
volucro del cristallo. Lo spessore di questo 
è sufficiente ad assorbire la radiazione BT 
emessa dal 19%Au. 

La fig. 3 mostra la taratura dello 
spettrografo riguardo alla 
energia. 


misura di 


A. BRACCI, U. FACCHINI €@ A. MALVICINI 


Questa misura & stata fatta utiliz- 
zando le righe di decadimento del 19% Au 
e del %%Co che sono molto ben cono- 
seiute (4). La linearita € buona entro 
1220 

Il fondo dell’insieme cristallo-foto- 
tubo è fatto di impulsi di ampiezza cor- 
rispondente a circa 20 keV. 

Il dispositivo usato per la misura dei 
y emessi nel processo di cattura di neu- 
troni lenti dal boro & illustrato nella 
fig. 1. 

La sorgente di neutroni & costituita 
da 15 mC di polonio mescolati con berillio 
in polvere. Questo tipo di sorgente ha 
il vantaggio di presentare una poco in- 
tensa radiazione y emessa direttamente 
dalla sorgente. Il cristallo & protetto da 
questa radiazione (y del polonio e y della 
reazione « + Be) per mezzo di uno spes- 
sore di piombo di 5 em. Il cristallo é 
pure protetto con piombo dai y emessi 
dai processi di cattura di neutroni nella 
paraffina. 

La disposizione a cavitä in parafiına 


© 
N 
N 
~~ 
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pan 
MeV 
n —— 4 a 
04n 1 Tass 
Fig. 3. — Curva di taratura dello spettrografo. 


La riga 0,411 MeV proviene dal decadimento 
del Au e quelle a 1,17 e 1,33 MeV dal decadi- 
mento del Co. La linearitä ¢ buona entro il 2%. 


(4) D. LIND, A. HEDGRAN: Ark. för 
5, 29 (1952). 
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permette di ottenere un buon flusso di 
neutroni nella zona vicino al cristallo. 

Uno strato di boro di 0,3 cm? av- 
viluppa linvolucro del cristallo. In que- 
sto modo i neutroni catturati nel boro 
producono la reazione in studio e nello 
stesso tempo, a causa del grande assor- 
bimento dello strato di boro, non pene- 
trano nel cristallo dove, catturati dallo 
iodio, produrrebbero radiazioni BT e y 
di disturbo. 

Uno spettro ottenuto con il disposi- 
tivo illustrato @ mostrato nella fig. 4. 

La riga corrispondente a 0,477 MeV 
è dovuta al decadimento dei nuclei di 7Li 
eccitati formati, come é noto, secondo la 
reazione 


VB + n> hit + « 
NUIT RAS 


Nel 93% dei casi di cattura di neu- 
trone si ha formazione di nuclei di “Li 
eccitati, mentre nei restanti si hanno nu- 
clei finali nello stato fondamentale (5). 
La media su 6 determinazioni ha dato 
per l’energia dei y il valore 0,477 MeV 
con il 2% di errore, in buen accordo con 
i valori noti della posizione del livello 
del 7Li e con l’unica misura diretta del- 
l’energia dei y fatta utilizzando una pila 
ed uno spettrografo magnetico (°°). 


(5) W. F. HORNYAK, T. LAURITSEN, P. MOR- 
RISON e W. A. FOWLER: Rev. Mod. Phys., 22, 
291 (1950). 

(5) L. S. ELLIOT e R. E. 
74, 1869 (1948). 


BELL: Phys. Rev., 
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La riga a circa 70 keV corrisponde a 
radiazioni X emesse dal piombo in se- 
guito ad effetti fotoelettrici provocati in 
esso da y di varia origine. 
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Fig, 4. — 4) Spettro ottenuto con la disposizione 
di fig. 1. La riga a 0,477 MeV proviene dai nuclei 
di 7Li formatisi a seguito della cattura di neutroni 
nel boro. La riga X a ~70 keV proviene da 
atomi di piombo che subiscono ionizzazione nelle 
orbite interne da parte dei raggi y provenienti 
dalla sorgente o prodotti dai neutroni in vari 
processi. B) Spettro del “Au per taratura. 


Ringraziamo il prof. G. Bora per le 
stimolanti discussioni. 
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Equivalenza fisica dell’equazione di Dirac-Corben 
con l’equazione di Dirac. 


E. BELLOMO e A. LOINGER 


Istituto di Fisica dell Universita - Pavia 


(ricevuto il 28 Ottobre 1952) 


Uno di noi (1) ha mostrato recentemente che le soluzioni dell’equazione (?) 


(1) (€0,0;7%;, + e Ag + eom,c? + o3m c?)yp = ih(0y/0t) 
2 goo )Y y 


che era stata detta di Dirac-Corben, si possono ottenere delle soluzioni dell’equa- 
zione di Dirac 


(2) (01057: + A5 + 03me?)y, = ih(Cy,/et) , 


(dove m & la massa della particella, legata ad m,, m, dalla m = (m? + m?)/?) me- 
diante la formula 


m + mo\/2 „een, eS 
, = = = 1 + 4 ———— 2) yp, = Up 1 
: | 2m | m+ Mm; 1 de 


con y e y, normalizzate ad 1. 
Vogliamo ora far vedere un pit stretto legame tra la (1) e la (2). 
La (1) si puö scrivere 
1 9 . 
C010 x7, + CAy + — (02M + 03M) me?| y = A(0y/0t) 
m 


Poichè il valore fisico degli operatori 04, o, dipende esclusivamente dalla loro 
posizione in un’equazione del tipo (2) e dalle loro proprietà di commutazione, ne 
viene necessariamente che g,— (1/m)(gym, + 3m) & l’equivalente fisico di = Pß. 
Ne consegue che nella teoria di Dirac-Corben ad ogni osservabile si deve far cor- 
rispondere un operatore, la cui espressione è la stessa di quella della teoria di Dirac 


() A. LOINGER: Nuovo Cimento, 9, 855 (1952). 
(@) H. C. CORBEN: Nuovo Cimento, 9, 591 (1952). 
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se in quest’ultima esso non contiene ß; se invece contiene ß deve farsi la sostitu- 
zione B-> 9,. Ad esempio, l’operatore del momento intrinseco di dipolo elettro- 
magnetico che in notazioni non covarianti ha la forma i«;%,B, %j)8 deve essere mu- 
tato in %a,;0;,05, «03 (*). 

Si osservi che eseguendo sugli operatori Ö e sulla y, della (2) la trasformazione 
unitaria Ö'= UGU-, y, = Uy, la B si muta in g, e si ottiene la (1) con y = y» (4). 

Pertanto la teoria basata sulla (1) pud considerarsi una rappresentazione, un 
po’ diversa dall’usuale, della teoria di Dirac: tutte le conseguenze fisiche della (1) coin- 
eidono dunque con quelle della (2) (*). 


(8) Il non fare questa sostituzione porterebbe in questo caso all’esistenza anche di un momento 
elettrico di dipolo per una particella libera a impulso cinetico nullo. L’aggiunta nell’equazione (1) 
di un termine di Pauli (come, ad es., è stato fatto da P. CALDIROLA e P. GULMANELLI in Nuovo 
Cimento, 9, 834 (1952)) soddisfa ai requisiti dell’invarianza relativistica ed è, per quanto detto, non 
equivalente all’aggiungere lo stesso termine all’equazione di Dirac (2). Tale termine puod infatti 
seriversi: 


— 1o;( Ex + 03By) = — 49% 1/m { Mol Ex + 0, By) + Mel 0,1% — 0! By)} 5 


il termine in m,/m è legato all’esistenza di un momento elettrico anomalo a riposo. 

(#) Si noti che U commuta con tutti i 0, o, ad eccezione che con 03, 03. 

(*) Il prof. CORBEN ci ha fatto osservare che, da quanto sopra, «segue che si puö separare 
arbitrariamente la massa della particella in due parti (mo; m,) scalare € pseudoscalare rispettiva- 
mente », richiamando una particolare attenzione sull’esistenza di una certa libertà nella formula- 
zione della teoria di Dirac. 


Ringraziamo il prof. CORBEN della cortesia e dell’interessamento che ci ha dimostrato. 
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Un caso d’invecchiamento di seintillatori liquidi. 


M. AGENO, G. CORTELLESSA e R. QUERZOLI 


Istituto Superiore di Sanita, Laboratorio di Fisica - Roma. 


(ricevuto il 5 Novembre 1952) 


E cosa ben nota che non ¢’ molto 
accordo tra i dati dei diversi Autori 
relativi all’andamento del potere fluore- 
scente con la concentrazione, nel caso 
di soluzioni utilizzabili come scintillatori 
liquidi. Questa mancanza di accordo viene 
generalmente attribuita al diverso grado 
di purezza dei solventi e dei soluti ado- 
perati, dato che, come ha dimostrato 
KALLMANN, basta l’agoiunta di piccolis- 
sime quantita di un secondo soluto oppor- 
tunamente scelto per modificare comple- 
tamente la curva di fluorescenza. 

Particolarmente notevole & tuttavia 
il disaccordo, nel caso delle soluzioni di 
naftalina in xilolo, studiate da noi (1) 
nel 1949, ristudiate da KALLMANN (2) nel 
1950-51 e dalla Cıccone (*) nel 1952. 
Mentre infatti la nostra curva (dopo la 
rapida salita iniziale) continua costante- 
mente a crescere linearmente fino alla 
concentrazione massima di 170 g/l, le 
curve di Kallmann decrescono dopo aver 
raggiunto un massimo e quella della Cic- 
cone mostra un potere fluorescente co- 
stante, indipendente dalla concentrazione. 

Per chiarire la causa di questo disac- 


() M. AGENO, M. CHIOZZOTTO e R. QUER- 
ZOLI: Rend. Acc. Lincei, 6, 626 (1949); Phys. 
Rev., 79, 720 (1950). 

(?) H. KALLMANN e M. FÜRST: Phys. Rev., 
79, 857 (1950); 81, 853 (1951). 

(5) A. CICCONE: Nuovo Cimento, 9, 36 (1952). 


cordo qualitativo, abbiamo preparato una 
serie di soluzioni di natfalina in xilolo, 
con naftalina bicristallizzata da solu- 
zione alcoolica (p.f. 80 °C) e xilolo puro 
per analisi della ditta C. Erba di Milano. 
Abbiamo eonservato tali soluzioni in bot- 


EMI - 473 


Fig. 1. 


tiglie di vetro con tappo di sughero e in 
normale ambiente di laboratorio, per 
circa dieci mesi, misurandone a intervalli 
di tempo il potere fluorescente. Le mi- 
sure sono state fatte nelle condizioni geo- 
metriche della fig. 1, usando per ecci- 
tare le scintillazioni una sorgente di 
raggi y costituita da un preparato di 
MsTh di circa 10 mC, e misurando la 
fluorescenza con un fotomoltiplicatore 
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EMI 473 a temperatura ambiente, con- 
nesso ad un elettrometro con resistenza 
di fuga. 

Come é mostrato nella fig. 2, la curva 


$ Soluzioni di naftalina in xilolo 


50 —o— misure del 26-12-51 
--o— misure del-20-10-52 


50 100 150 


Fig. 2. 


di fluorescenza ottenuta subito dopo la 
preparazione delle soluzioni coincide pra- 
ticamente con quella da noi ottenuta 
nel 1949, pur essendo il dispositivo spe- 
rimentale completamente diverso. La 
curva ottenuta dopo dieei mesi ha invece 
un andamento assai simile a quello ot- 
tenuto da KaALLMANN nel suo primo 
lavoro. Ë da notare che le irregolaritä 
nell’andamento di questa seconda curva 
sono reali e non dovute ad errori nella 


della fluorescenza, come se il 
processo di invecchiamento avesse pro- 
ceduto in modo abbastanza irregolare e 
diverso da soluzione a soluzione. Resul- 
tati del tutto analoghi si sono ottenuti 
con altra serie di soluzioni preparate con 
naftalina commerciale. 

E da notare che le soluzioni a con- 
centrazione maggiore hanno assunto col 
tempo una leggera colorazione giallina 
e che, svaporata parte della soluzione à 
150 g/l, si & ottenuto un residuo di p-f. 
74°C. Si tratta quindi certamente di 
un’alterazione di natura chimica, che puo 
essere dovuta o alla formazione di pro- 
dotti di ossidazione o a lenta reazione 
della naftalina con impurezze origina- 
riamente contenute nello xilolo, forse 
sotto l’azione della luce. 

Crediamo comunque non inutile segna- 
lare il fatto, che pud, a nostro avviso, 
spiegare l’origine del disaccordo di cui 
abbiam parlato all’inizio, anche per Vin- 
teresse che pud presentare nei riguardi 
dei contatori a scintillazione a liquidi, 
non essendo escluso che fatti analoghi si 
possano verificare con altre soluzioni 
scintillanti. 


misura 
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Slow Heavy Mesons from Cosmic Ray Stars. 


R. Levi SETTI 


Istituto di Scienze Fisiche dell Universita - 


Milano 


Istituto Nazionale di Fisica Nucleare - Sezione di Milano 


G. TOMASINI 


Istituto di Fisica dell Universita - Genova 


(ricevuto 18 Novembre 1952) 


Two stars have been observed, in 
plates exposed at a maximum altitu- 
de of about 90000 feet, showing the 
ejection of a slow heavy meson which, 
after being brought to rest in the emuls- 
ion, disintegrates with the emission of 
one charged secondary particle. 

The plates were Ilford G5, of 1200 u 
and 600 u thickness, respectively; they 
have been exposed in the stratosphere 
by means of balloon launchings from 
Cagliari, in the course of the scientific 
expedition which took place there during 
the summer 1952, as a collaboration 
between the University of Bristol, the 
Centre de Physique Nucléaire of Brussels, 
the Italian Universities, other European 
laboratories and the Italian Navy and 
Air Force. 

Both stars occurred in heavy nuclei 
of the emulsion. The events are partially 
shown in the photomicrographs of fig. 1 
and 2; in the following they will be 
called as event N. 1 and N. 2. 


— The star has 
of which 


Event N. 1 (fig. 1). 


38 visible prongs, 13 are 


black, 10 grey and 15 at minimum ioni- 
zation. The star primary is most likely 
an «partie. The heavy meson is 
ejected in the backward direction, at an 
angle of about 80° to the primary, it 
stops in the emulsion after a path of 
5260 u and decays into a fast secondary 
particle 450 y long. An accurate direct 
measurement of the mass of the particle 
has been possible. The mass has been 
determined by scattering-range and by 
gap-range measurements. 

The method employed for scattering- 
range determination of that 
described by GOLDSCHMIDT et al. (1): the 
track has been divided into segments, 
the scattering measured with constant 
cell length for each segment, beginning 
with short cells in the segments nearest 
the slow end of the track and increasing 
the cell length going toward the fast end. 


mass 18 


@) Y. GOLDSCHMIDT-CLERMONT, G. W. KING, 
H. MUIRHEAD and D. M. Ritson: Proc. Phys. 
Soc. 61, 183 (1948); Y. GOLDSCHMIDT-CLER- 
MONT: Nuovo Cimento, 7, 331 (1950). 
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Fig. 1. — A heavy meson of mass 1049 + 9) m, is ejected by a star and decays 
into a charged secondary particle. (Observed by R. LEvI SETTI). 


R. Levi Serti and G. ToMASINI 


100 u 


Fig. 2.- A heavy meson of mass 1380 + 210 m, is ejected 
a charged secondary particle. (Observed by G 


by a star and decays into 
. TOMASINI). 
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The mass was determined for each seg- 
ment by taking into account the va- 
riation of energy with residual range and 
finally a weighted mean was taken. The 
cell length has been varied in a way 
such as to keep about constant the si- 
gnal to noise ratio and chosen to be 
the optimum cell. 

Angular and coordinate measure- 
ments were made, the value of the mass 
has been found to be 1030 + 120 m, as 
a mean of several determinations. 

The mass has also been determined 
by gap-length measurements using a 
method due to one of us (R. L. S), the 
«x method». The total gap length per 
unit length is determined by means of 
an ocular movable 
thread on a minimum ionization x par- 
ticle. An « particle, primary of a high 
energy star, satisfies this requirement. 
This furnishes a calibration of the plate 
at a ratio //I,= 4, on a track of con- 
stant gap-density. Such a measure will 
define a 45° straigth line in doublelog 
paper if the total gap length is plotted 
against range. The same gap measur- 
ement is made on the unknown track 
and the result plotted. The intersection 
of the straight line given by the latter 
measurement with the 45° line defines 
the range where the unknown particle 
has an ionization equal to 4 times the 
minimum value. The mass is then 
simply calculated. This method yielded 
for the particle in question a mass of 
1050 +140 m,. The error on the mass 
can be reduced by taking the mean of 
the two independent determinations and 
the result is: 1040 +90 m,. 

The secondary particle has a grain 
density indistinguishable from the mi- 
nimum, G/G, = 1.15 + 0.15. We con- 
sidered that, owing to the dip and 
insufficient length available, the scatter- 
ing measurements we have made and 
which give a lower limit of 80 MeV/c 
for the pp, are not trustworthy. How- 
ever the mass of the primary particle 
and the observed specific ionization 


micrometer with 
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of the secondary, require that this latter 
has a mass inferior to 400 m4. 

The energy of emission of the de- 
scribed heavy meson, in the laboratory 
system, is about 28 MeV; the lower limit 
of its lifetime can be determined and 
results to be about 107s. The star 
shows also, may be as a coincidence, 
the emission of a slow ~~ meson of 
2.4 MeV energy, ending with a blob, at 
an angle to the primary simmetrical to 
that of emission of the heavy meson. 

Event N. 2 (fig. 2). — The star 
has a total of 40 branches, respectively 
25 of evaporation, 3 grey and 11 sho- 
wer particles; the primary is most pro- 
bably a charge 1 particle, indicated as 
P in fig. 2. The ejection of the heavy 
meson occurs also in this case in the 
backward direction, at an angle of 16° 
to the probable primary; the particle 
ends in the emulsion after 1 885 y and 
decays into a secondary, 840 y long. 
Scattering-range determination of the 
mass, as described above, together with 
gap measurements give à 
13804210 m,. The 
length of the secondary particle of this 
meson, in comparison with the secon- 
dary of event N. 1, allows a tentative pf 
determination to be made by scattering. 
Several measurements give a value of 
100 +15 MeV/e. We consider that the 
error on this figure can be greater than 
the statistical one. Also this secondary 
is at minimum ionization, @/@, = 1.03 + 
+- 0.5. Its mass can be assumed to be 
probably less than 300 me. The energy 
of emission of the heavy meson is about 
18 MeV, the lifetime is greater than 
5-10-11s. Some considerations on the 
energy involved in the two parent stars 
follow. 

Further measurements are in pro- 
eress on some of the shower particles 
ejected by the two stars, the energy of 
the incident particles cannot be, however, 
directly measured owing to the insufficient 
length of the tracks. A rough estimate 


value of 
increased track 
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of the lower limit of this energy can be 
attempted from considerations in the 
C-system; an energy of about 2 GeV is 
at least required for each shower meson 
produced, giving a value of the order 
of 30 GeV for the « particle primary 
of the event N. 1, about 25 GeV for the 
primary of event N. 2. An estimate of 
the primary energy from the angular 
distribution of the shower particles has 
not been attempted, since the events are 
complex, giving evidence for secondary 
processes inside the nucleus; in fact 
there is a large spread in the angles of 
emission of the shower particles, some 
of them being ejected in the backward 
direction, and a large number of heavily 
ionizing particles is also present. By con- 
sidering the relation between the aver age 
multiplieity of shower particles and the 
primary energy, as shown by CAMERINT 
et al. (?), one gets a probable value of 
the order of 50 GeV for the primary 
of star N. 1 and of about 45 GeV for 
N. 2. Since the primary of star 
N. 1 is an « particle, it follows that both 
mesons have been produced in nucleon- 
nucleus collisions at an energy inferior 
than 50 GeV. This conclusion however 
must be regarded as only tentative. 

Data on the production of heavy me- 
sons in high energy nuclear interactions 
have been reported by DANIEL et al. (°). 
with regard to fast particles identificat- 
ion. From these observations it results 
that in stars accompanied by the emis- 
sion of large number of heavy prongs, 
the production of heavy mesons must 
be very low at energies inferior to about 
50 GeV. 

The results of the measurements on 
the two events described here are col- 
lected in the following table I. 


star 


(?) U. CAMERINI, J. H. DAVIES, C. FRANZI- 
NETTI, W. O. Lock, D. H. PERKINS and G. 
YEKUTIELI: Phil. Mag., 42, 1261 (1951). 

6) R. R. DANIEL, J. H. DAVIES, J. H. MUL- 
VEY and D. H. PERKINS: Phil. Mag., 43, 753 
(1952). 
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TABLE I. 


Heavy mesons. 


Event N. 1 | Event N. 2 | 

| | 

Range | 5 260 u 1885 u | 

mass | 1040-+90m, | 1380+210m, | 

lifetime | > 10710 3 >5-10 48 | 

Secondary particles. 

| | 

Event N. 1 | Event N. 2: | 

Length 4504 | 840 U. | 
G/G, C1642 015 1.03 £0.05 

pp ? 100-15 MeV/e 

mass < 400 m, < 300 m, 


Observations on the decay at rest of 
heavy charged mesons coming from out- 
side the emulsion, into a single visible 
secondary, in the photographic plate, 
have been described by O’CEALLAIGH (¢), 
more recently by PoweLL (5). The re- 
sults of these authors show the pro- 
bable existence of two types of heavy 
particles, the x meson, of mass around 
1050 m,, decaying into a u meson and 
more than one neutral particle, and the 
xy meson, of around 1450 me, 
decaying into a x meson and probably 
a single neutral particle. 

Production and subsequent decay at 
rest of slow heavy mesons (S particle) 
has been up to now observed in some 
cloud chamber photographs by BRIDGE 


mass 


(*) C. O’CEALLAIGH: Phil. Mag., 42, 1032 
(1951). 

6) C. F. POwELL: Köbenhavn Conference, 
1952 
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and Annis (%) and by ANNIS et al. (7); 
the behaviour of charged heavy mesons, 
called V+ particles, has been described 
by ARMENTEROS et al. (5). 

Our measurements show that heavy 


() EH. S. BRIDGE and M. ANNIS: Phys. Rev., 
82, 445 (1951). 

() M. ANNIS, H. BRIDGE, H. COURANT, S. 
OLBERT and B. Rossi: Nuovo Cimento, 9, 624 
(1952). 

(6) R. ARMENTEROS, K. H. BARKER, C. C. 
BUTLER, A. CACHON and A. H. CHAPMAN: Na- 
ture, 167, 501 (1951); R. ARMENTEROS, K. H. 
BARKER, C. C. BUTIER, A, CACHON and C. M. 
York, Phil. Mag., 48, 597 (1952). 
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mesons can be directly produced in stars; 
the evidence given by the mass mea- 
surements on the V like phenomenon 
N. 1 allows us to consider this event as 
belonging to a (Sx)-class; whereas event 
N. 2 seems to fall into a class (Sy). 
However, this attribution is less certain 
that in the event N. 1, owing to the 
shorter track. If a decay into two par- 
ticles is assumed, the second event, al- 
though with uncertainty, seems to indi- 
cate that the neutral particle ejected in 
the disintegration ought to be a heavy 
neutral meson V2 of mass around 800 me. 
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Anomalies of the absorption curve of Cosmic radiation in lead. 


M. App EL-WAHAB KHALIL 


Physics Department, Faculty of Science University - Alexandria, (Egypt). 


(ricevuto il 9 Novembre 1952) 


Many authors, investigating the ab- 
sorption of the vertical component of 
cosmic radiation in lead by a Geiger- 
Müller counter telescope, found small 
anomalies in the absorption curve (1.2.3.4), 
while others denied these changes of 
absorption (5). 

It was said that there is a relation 
between these anomalies in the absorpt- 
ion curve and the maxima observed in 
the transition curve of Rossi. 

First GEORGE and APPAPILLAI (1) 
found in 1945 at sea-level (London and 
Colombo) that the absorption curve 
showed a plateau at about 10 em lead, 
and a small anomaly between 15 and 
20 cm. 

KELLERMANN and WESTERMAN (?) 
using à method which avoided largely 
the influence of scattering have investig- 
ated the change of absorption found by 
them in the region of 12 em Pb. A 


(!) E. P. GEORGE and V. APPAPILLAI: Na- 
ture, 155, 726 (1945). 
(?) E. W. KELLERMANN and K. WESTER- 


MAN: Proc. Phys. Soc., A 62, 354 (1949). 

(?) E. FENYVES and O. HAIMAN: 
165, 244 (1950). 

() E. MazzoLLı1 De Marnov: Nature, 167, 
192 (1951). 

(°) G. R. HEYLAND and W. E. DUNCANSON : 
Nature, 167, S95 (1951). 


Nature, 


specifie change of the absorption coef- 
ficient was found also in the 22 cm Pb 
region. They interpreted the first change 
as due solely to the overlapping of the 
high momenta end of the electron 
spectrum and of the meson spectrum. 

Repeating the experiment of GEORGE 
and APPAPILLAI in Budapest, and elimin- 
ating barometric and temperature effects, 
FENYVES and Harman (?) found a second 
anomaly at 26.6 cm lead beside that 
mentioned by GEORGE and APPAPILLAI. 
The first anomaly was smaller than the 
standard deviation, while the 
was 2.5 times greater. 

De Marnoy (?) in 1951 has found a 
maximum of intensity under 13 cm 
lead + 30 em concrete roof. 

In the same year HEYLAND and Dun- 
CANSON (5) reported that the intensity of 
cosmic radiation decreased continuosly 
with increasing thickness of the absorber 
over a range of 30 cm lead. 

Repeating the experiment in Alexan- 
dria - Egypt (latitude 31°,2 N, iongitude 
290,9, geomagnetic latitude 280,5 N), I 
have found a third anomaly at 37.5 em Pb 
beside those at 15 cm and 27.5 cm Pb. 

Using a counter telescope of four tu- 
bes, the top and the bottom tubes being 
90.5 em apart, each counter tube of 
effective length 60 cm, diameter 4 em, 


second 
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a fourfold coincidence was established 
with a Rossi’s circuit of resolving power 
of 10-5 s. The curve obtained is shown 
in fig 1. From this curve we noticed that: 

1) The 2 first anomalies at 15 
and 27.5 cm lead are just within the 
standard deviation of our measurement 
untill today. 

2) The absorption curve for the 
soft component of cosmic radiation is 
in coincidence with that made by 
FEnyvss and HAIMANN. 

3) Making the experiment on the 
roof under 3 mm aluminium and down- 
stairs under a thickness of 40 em con- 
crete, it was noticed that these anomalies 
are more distinguished downstairs, the 
same as in the case of maxima of Rossi's 
transition curve. 

It should be 
this work is 


mentioned here that 
an extension to our (°) 
work in (1951) on the maxima found at 
about 15.25 and 37.5 cm Pb of Rossi’s 
transition curve. 

‘The second and the third maxima of 
Rossi’s curve were also found by many 
authors like Cray (7), Borue and 
Turn (8), CHAUDHURY (?) and CAMEDA- 
Miura (2°). 

This paper confirms the relationship 
between the anomalies of the absorption 
curve and the maxima found in the 
transition curve. 

It also confirms the existence of a 
fourth maximum at about 37.5 em Pb 
which till now is only found by us. 


(°) L. ALLEGRETTI e M. ABD EL-WAHAB 
KHALIL: Nuovo Cimento, 8, 1011 (1951). 

(7) J. Chay and W. L. SCHEEN: Physica, 
14, 489 (1948). 
(83) W. Borne and H. THURN: Phys. Rev., 
545 (1951). 
(HAE EL. 
274 (1951). 
(4%) T. KAMEDA and 1. Miura: Journ. Scien. 
Res. Inst. Tokio, 45, 65 (1951). 


79, 


SEN CHAUDHARY: Phys. Rev., 


81, 
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An interpretation for these observed 
anomalies, and for the relation men- 
tioned above can be shown from the 
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two following cases concerning showers 
created inside lead. 

a) Shower created inside lead by 
a charged particle incident from a 
direction inclined at a small angle to 
the vertical. This shower has a com- 
ponent which can strike the 2 lower 
counters at the same time the incident 
particle strikes the 2 above counters. 

b) Small angle shower created inside 
lead by an incident component of cosmic 
particles which passes only through the 
2 above counters and is stopped inside 
lead. The secondary particles strike 
the 2 lower counters at the same time. 

Thus this arrangement for measur- 
ing absorption behaves as a special case 
for the arrangement of Rossi for detecting 
transition showers. 


I thank Prof. Dr. L. ALLEGRETTI 
for his useful suggestions, and also Prof. 
Dr. M. Ex-Suprerni the director of the 
department of physics. 


LIBRI REC ENUPLI i> RE CENSTONGE 


C. MOLLER — The Theory of Relativity. 
Oxford, At the Clarendon Press, 
1952, pp. x11+ 386. 


Il presente volume, basato su un 
corso svolto ormai vent’anni fa, dall’ Au- 
tore all’Universita di Copenaghen, colma 
nel miglore dei modi una lacuna della 
letteratura fisica della prima meta del 900; 
è difficile infatti trovare un volume sulla 
teoria della relatività che contenga côn- 
temporaneamente la teoria speciale e 
generale e che accumuli il pregio di essere 
scritto con intento didattico con quello 
dell’estrema completezza. 

Il primo capitolo del volume & dedi- 
cato delle difficoltä che si 
presentavano nella teoria dell’etere. 

La risposta negativa alla domanda: 
«Permette una qualche esperienza di 
determinare il moto della terra rispetto 
all’etere? », conduce ad ammettere che 
il principio di relatività abbia validita 
generale; e se si ammette 
validità delle equazioni di Maxwell si & 
condotti al principio della costanza della 
velocità della luce e di qui (Cap. II) 
alla critica del concetto di simultaneita 
e alla deduzione delle trasformazioni di 
Lorentz. 

Le utili formule della cine matica rela- 
tivistica vengono poi stabilite, in parti- 
colare la formula di Thomas. 

Nel cap. III sono trattate le equa- 
zioni del moto della dinamica relativi- 
stica e le formule di trasformazione per 
le grandezze dinamiche pit importanti e 


all’esame 


ancora la 


si accenna a quelle che sono le verifiche 
sperimentali delle formule stesse : l’angolo 
formate dalle direzioni del moto di due 
elettroni urtantisi (CHAMPON), la defles- 
sione degli elettroni in un tubo a raggi 
catodici, l’energia di legame dei nuclei. 

Il cap. IV tratta la formulazione qua- 
drimensionale della teoria della relativita ; 
la covarianza delle leggi della natura 
rispetto alla trasformazione di Lorentz 
si esprime in effetti nel modo più sem- 


plice attraverso equazioni del tipo 
tru, = 0, dove t.,, © un quadritensore. 


Come l’Autore avverte nella prefazione, 
nel presente volume il formalismo qua- 
dridimensionale è tuttavia usato meno 
di quanto si faccia di solito. Sebbene 
ovviamente tale formalismo sia lo stru- 
mento pit indicato per trattare questioni 
di carattere generale, lo spazio ed il 
tempo, quali noi li concepiamo, sono enti 
fisicamente diversi e in una formulazione 
didattica della teoria è conveniente far 
risaltare questa differenza. Del resto, lo 
spazio che l’Autore dedica alla formula- 
zione quadrimensionale della teoria è più 
che sufficiente a fare risaltare l’utilità 
del formalismo in questione. 

In tutto il volume le due necessita, 
quella di separare nelle formule la parte 
spaziale e quella temporale e quella di 
far uso invece del formalismo a quattro 
dimensioni, si conciliano perfettamente. 

L’invarianza della carica elettrica, 
conseguenza della validità in ogni sistema 
di riferimento dell’equazione di conti- 
nuità, porta con sé che i potenziali sca- 
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lare e vettore costituiscono un quadri- 
vettore ed apre la strada alla scrittura 
delle equazioni di Maxwell in forma co- 
variante a vista. 

Queste equazioni vengono trattate nel 
cap. V (Elettrodinamica del vuoto). Il 
formalismo quadrimensionale si rivela 
qui, come & noto, quasi indispensabile, 
permettendo una rapidissima deduzione 
dei potenziali di Lienard e Wiechert e 
di tutte le questioni connesse. 

La trattazione dei sistemi chiusi, cioè 
di quei sistemi per i quali valgono i 
classici teoremi di conservazione, per- 
mette all’Autore di trattare la questione 
di definizione di baricentro di un sistema 
relativistico; come caso particolare di un 
sistema chiuso viene poi considerato un 
mezzo elastico e le equazioni del moto 
di tale mezzo sono stabilite. 

Un esempio, infine, di sistema con- 
tinuo relativistico è offerto da un qual- 
siasi campo quando la Lagrangiana ad 
esso associata è relativisticamente inva- 
riante. Questa questione è trattata nello 
stesso capitolo VI ed il tensore energia- 
impulso del campo viene costruito a 
partire dalla Lagrangiana. Segue poi nel 
capitolo seguente la trattazione di un 
insieme di questioni di solito poco di- 
seusse: l’elettrodinamica nei dielettrici e 
nei mezzi paramagnetici e la termodica- 
mica relativistica; particolarmente inte- 
ressante lo studio delle proprietà di tra- 
sformazione delle varie grandezze termo- 
dinamiche. 

Col cap. VIII si apre la trattazione 
della teoria generale della relativita che 
occupa i successivi quattro capitoli. Nel 
capitolo in questione vengono, come è 
naturale, trattati i campi gravitazionali 
non permanenti, cioè quei campi che si 
presentano quando si è fuori dalla pre- 
senza immediata di masse gravitazionali. 
In tal caso, se ci si mette in un sistema 
inerziale, la geometria dello spazio-tempo 
è la geometria pseudo-euclidea, onde lo 
studio dei campi gravitazionali non per- 
manenti non implica l’introduzione di 
una geometria di Riemann. Una tale 
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geometria è invece quella che si pre- 
senta quando si è in presenza di masse 
oravitazionali a distanza finita (campi 
eravitazionali permanenti). In tal caso 
non & più possibile introdurre un sistema 
di coordinate che renda le componenti 
del tensore g,, indipendenti dal punto; 
ma & sempre possibile far si che ciù sia 
vero in tutto un intorno di ciascun punto, 
cioè introdurre i sistemi locali di inerzia. 
Cid rende possibile, ammesso il principio 
di equivalenza, ricavare anche in questo 
caso le leggi fisiche a partire da quelle 
note dalla teoria speciale, pur di svi- 
luppare alquanto il ealcolo_ tensoriale; è 
cid che si fa nel cap. IX e i risultati ivi 
stabiliti servono (nel cap. X) per mo- 
strare come, ammesso di conoscere la 
metrica, le equazioni fondamentali della 
meccanica e dell’elettromagnetismo pos- 
sano essere immediatamente seritte. 
Resta il problema di determinare la me- 
trica; questo viene trattato nel cap. XI 
seguente, nel quale si mostra come, sup- 
posto di conoscere il tensore energia-im- 
pulso del sistema, i g,, Siano determinati; 
essi (ossia il ds?) vengono poi determinati 
in vari casi interessanti ed in partico- 
lare, seguendo SCHWARZSCHILD, nel caso 
di una distribuzione di massa sferica. 
Quanto stabilito permette di calcolare 
lo spostamento del perielio di Mercurio 
che, assieme allo spostamento verso il 
rosso e alla deflessione gravitazionale 
della luce, costituiscono le tre esperienze 
di prova per la teoria generale della 
relativita; problemi che vengono diseussi 
nell’ultimo capitolo. Il volume si chiude 
con una discussione sui modelli cosmo- 
logici. 

La completezza dell’opera, che spe- 
riamo sia risultata evidente da quanto 
abbiamo qui scritto, rende augurabile 
che essa incontri la pit larga diffusione; 
essa & veramente utile a chiunque si 
oceupi di Fisica: sia egli ancora studente, 
o gid dedito alla ricerca sperimentale, o 
teorico. 


G. MORPURGO 


H. S. W. Massey e E. H. S. BURHOP 
— Electronic and Ionic Impact Phe- 
nomena, pp. XVIH+669. At the 
Clarendon Press, Oxford, 1952. 


Questo nuovo grosso volume di The 
International Series of Monographs on 
Physics sara sicuramente accolto con 
grande interesse (stavamo per scrivere: 
con gioia) da tutti coloro che si trovano a 
dover combattere con le difficolta e coi 
problemi di una tecnica in cui inter- 
urti tra particelle di bassa 
energia: cioé dalla quasi totalita dei 
ricercatori. E infatti ben noto che i 
fenomeni d’urto di elettroni e di ioni 
tra loro, contro atomi e molecole neutre 
© contro superficie cristalline interven- 


vengano 


gono in modo diretto o indiretto e spesso 
in modo essenziale in quasi tutti i dispo- 
sitivi oggi in uso per la rivelazione e lo 
studio dei fenomeni elementari, nonchè 
in un numero sempre crescente di altri 
dispositivi, destinati a studi e ricerche di 
ogni genere. L’argomento trattato ha 
d’altra parte anche un notevolissimo inte- 
resse intrinseco, per le sue relazioni con 
la stiuttura dell’atomo e delle molecole 
e con le basi sperimentali della teoria 
dei quanti. Una limpida ed ordinata 
esposizione dello stato presente delle 
nostre conoscenze in materia, che pre- 
senti in un giusto equilibrio tutti i pit 
importanti resultati sperimentali da un 
lato e daltro lato (quando siano possi- 
bili) le impostazioni teoriche delle sin- 
gole questioni e i relativi sviluppi, por- 
tera quindi certamente un assai notevole 
contributo al progresso della ricerca e 
al chiarimento dei molti punti ancora 
oscuri od incerti. 

Il volume di cui ei stiamo occupando 
risulta diviso in due parti, di ampiezza 
poco diversa, la prima delle quali com- 
prende i primi sei capitoli (361 p.) e 
tratta per estesc i fenomeni d’urto in 
cui intervengono elettroni. Gli ultimi 
quattro capitoli, costituenti la seconda 
parte (285 p.), sono relativi ai fenomeni 
d’urto tra atomi ed ioni nei vari casi e 
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condizioni possibili e comprendono la 
trattazione dei fenomeni di ricombina- 
zione. 

La prima parte inizia con un capitolo 
relativo alla definizione e alla misura 
della sezione d’urto totale di elettroni 
contro atomi. Seguono due capitoli in 
eui l’urto di elettroni contro atomi viene 
diffusamente studiato, prima dal punto 
di vista sperimentale (II) e poi dal punto 
di vista teorico (III). Vengono descritti 
dapprima i metodi per la misura delle 
sezioni d’urto per ionizzazione e per ecci- 
tazione. Seguono lo studio della polariz- 
zazione della radiazione emessa per urto 
elettronico e della distribuzione angolare 
degli elettroni diffusi, sia elasticamente 
sia anelasticamente. La parte teorica 
inizia con la trattazione quantistica della 
diffusione elastica, col calcolo della 
relativa sezione d’urto e della distribu- 
zione angolare degli elettroni diffusi. Ven- 
gono quindi trattavi nell’approssimazione 
di Born gli urti anelastici e si discute il 
confronto tra le sezioni d’urto calcolate 
e i dati sperimentali disponibili. Tra gli 
argomenti vari, trattati alla fine di questo 
»apıtolo, segnaliamo un paragrafo sul 
limite teorico di potere risolutivo per 
un microscopio elettronico. 

Il capitolo IV riguarda l’urto di elet- 
troni contro molecole. Viene trattata dap- 
prima la diffrazione di elettroni, veloci 
e semiveloci, da parte di molecole e poi 
la diffusione di elettroni lenti. Tutto il 
resto del capitolo & dedicato alla ecci- 
tazione e dissociazione di molecole per 
urto elettronico: i dati disponibili rela- 
tivi ai vari gas e vapori sono qui esposti 
e discussi diffusamente. 

Il capitolo successivo (V) & dedicato 
all’emissione di elettroni secondari da 
parte di superficie sottoposte a bombar- 
damento elettronico. Sono esposti i me- 
todi di misura e le teorie dell’emissione 
secondaria e quindi discussi i dati spe- 
rimentali. Viene infine un capitolo (VI) 
dedicato agli urti elettronici con emis- 
sione di radiazione: si calcolano e discu- 
tono le sezioni d’urto e infine si tratta, 
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sia dal punto di vista teorico sia da 
quello sperimentale, la radiazione di fre- 
namento. 

Col capitolo VI ha inizio quella che 
abbiamo chiamato la seconda parte del 
volume. Esso tratta ampiamente gli urti 
tra atomi e molecole in condizione di 
disordine elementare. Le interazioni tra 
molecole sono studiate prima a partire 
dai dati macroscopiei quali la viscosità 
e poi dalle misure dirette di sezioni d’urto 
totali. Segue un paragrafo sulla mobilitä 
degli ioni positivi e quindi sono ampia- 
mente trattate le collisioni anelastiche 
accompagnate da transizioni elettroniche 
o da variazioni dello stato di vibrazione 
e rotazione delle molecole. 

I due capitoli successivi riguardano 
il passaggio di un fascetto omogeneo di 
raggi positivi o neutri attraverso gas 
(VIII) e l’urto di un tale fascetto contro 
una superfieie (IX). Discussi i metodi 
per la produzione dei raggi in questione, 
si espongono i metodi sperimentali usati 
e i resultati ottenuti per la misura di 
sezioni d’urto elastiche e anelastiche. Si 
tratta quindi elettroni 
secondari per bombardamento ionico € i 
vari fenomeni di riflessione e variazione 
di carica di ioni 
superficie. Un paragrafo à dedicato alla 
interazione tra ioni positivi e superficie 
catodica nei contatori di Geiger. Viene 
infine trattata l’evaporazione catodica e 
l’urto di molecole gassose contro le pa- 
reti del recipiente. 

L’ultimo capitolo (X) &, come gia si 
© detto, dedicato ai fenomeni di ricom- 
binazione tra ioni positivi e negativi e 
tra ioni positivi ed elettroni. Coneludono 
il volume aleune appendici, di eui una 
relativa a resultati moltc recenti, oltre 
ai soliti indici. 


Vemissione di 


nell’urto contro una 


Come si vede, il materiale raccolto e 
ordirato in questo volume & imponente 
e cid fa si che esso risulti un manuale 
di consultazione e di orientamento sulle 
singole questioni veramente prezioso. La 
chiarezza dell’esposizione, la ricchezza di 
dati numerici selezionati, di grafici e 
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tabelle lo rende inoltre di uso facile e 
di piacevole lettura. Per cid che riguarda 
la sostanza, il recensore ha invano cer- 
cato qualche svista da segnalare: se pur 
ye, essa è al di là delle sue cognizioni 
e capacità critiche. Non resta quindi 
altro che lamentare l’esclusione da questo 
bellissimo volume di alcuni argomenti, 
tra i quali primo quello degli urti tra 
molecole nei liquidi, argomento che va 
attirando oggi sempre più l’attenzione 
dei ricercatori. 

M. AGENO 


A. Mercier — Leçons sur les prin- 
cipes de l'électrodynamique  clas- 
sique. Pagg. 74. Édition du Griffon, 
Neuchâtel - Suisse 1952. 


Il programma che si & proposto l’Au- 
tore nello scrivere questo libro & chia- 
ramente delineato nella prefazione, do- 
vuta a L. DE BROGLIE. In essa & detto 
che l’Autore ritrova e giustifica le for- 
mule e le equazioni fondamentali delle 
teorie di Maxwell e Lorentz partendo 
da ragionamenti molto generali e preci- 
sando la parte che, in buona logica, è 
necessario riservare alla costante dielet- 
trica e alla permeabilità magnetica del 
vuoto. Un tale programma è svolto par- 
tendo dal convincimento che «una teoria 
ben fatta non deve prendere un fatto 
sperimentale per postulato: essa deve po- 
stulare proposizioni veramente prime, di 
per se stesse non. dimostrabili e diretta- 
mente non verificabili dall’esperienza; 
quindi dedurne relazioni tra grandezze 
misurabili, relazioni queste direttamente 
verificabili dall’esperienza ». 

Questa strana idea, dedotta forse da 
una vecchia classificazione delle teorie 
indescrittive ed esplicative, insostenibile, 
perché i «fatti» verificati o supposti, 
che stanno a base della teoria sono in 
realtà sempre ipotesi, almeno di genera- 
lizzazione, svela la posizione mentale del- 
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l’Autore che coi termini del problema 
degli Universali si chiamerebbe «rea- 
lismo ». 

Cid si conferma dovunque l’Autore 
parla di dimensioni delle grandezze fisiche 
e di costanti d’interazione. Qui egli 
prende continuamente per necessità lo- 
gica fatti di pura convenienza. Egli non 
dice per es. che una «costante d’intera- 
zione » c’è anche nella legge della dina- 
mica: F =k d(mv)/dt, e che per pura con- 
venienza noi la poniamo uguale ad uno, 
cosa che per es. non facciamo nella legge 
di gravitazione universale: F=G m,m,/r?, 
rinunciando per comodita a ridurre a 
due sole le grandezze indipendenti del 
nostro sistema di unita e sostituendo a 
un campione di massa difficilnente ripro- 
ducibile un campione completamente 
arbitrario. 

Anche la dove l’Autore, per dimo- 
strare l’assurdità logica dei due sistemi 
elettrostatico ed elettromagnetico, af- 
ferma che «una carica non pud essere 
uguale al prodotto di una carica per una 
velocità », c’& evidentemente la 
grossolana confusione d’idee.Si dimentica 
ancora una volta che le grandezze fisiche 
sono formazioni concettuali che derivano 
da un gruppo di operazioni effettivamente 
eseguibili e che cambiando il gruppo di 
operazioni, cambia il concetto corrispon- 
dente, anche se per avventura non se 
ne cambia il simbolo. Consigliamo all’ Au- 
tore di leggere almeno « La logica della 
fisica moderna» di BRIDGMANN. 


stessa 
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Con questo bagagho d’idee, l’Autore 
va alla ricerca delle proposizioni « vera- 
mente prime», «indimostrabili », neces- 
sarie per rendere l’elettromagnetismo una 
«teoria ben fatta ». Tali proposizioni sa- 
rebbero di due categorie. Una prima 
categoria comprende alcune « condizioni » 
generali ricavate in parte dalla necessità 
di introdurre per ogni nuova branca 
della fisica, almeno due (?) grandezze; ed 
in parte da considerazioni dimensionali. 
Una seconda categoria è data dalle condi- 
zioni particolari dell’elettricità che sono: 
la specificazione delle grandezze fisiche 
nuove proprie dell’elettromagnetismo e 
la constatazione della conservazione della 
carica e della propagazione dei campi 
con velocita finita. 

La costruzione, a partire dalle propo- 
sizioni prime, delle equazioni di Maxwell 
e Lorentz procede dalle equazioni di 
Coulomb-Gauss, ricavate facendo uso di 
condizioni di invarianza, che si suppon- 
gono necessariamente verificate dalle 
legei della fisica. 

Una volta giunti, con la trattazione 
contenuta nei primi capitoli, alle equa- 
zioni che legano i campi nel vuoto, 
l’estensione al caso della materia è fatto 
seguendo le linee classiche. 

Nettamente distaccato dal resto del 
libro € ultimo capitolo, che contiene un 
tentativo di assiomatica condotto facendo 
uso dei numeri di Clifford. 
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